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1.0  INTRODUCTION 


Over  the  years,  airborne  threats  to  the  United  States  and  its  forces  in  the 
field  have  evolved  to  include  low  RCS  technology  which  makes  these  threats 
much  more  difficult  to  detect  by  radar  systems.  The  small  RCS  detection 
problem,  per  se,  is  being  pursued  by  other  government  and  industry  programs 
and  was  not  the  primary  requirement  addresed  during  this  program.  However, 
when  such  low  RCS  threats  are  in  close  proximity  to  larger  targets,  other  factors 
are  present  which  exacerbate  the  detection  problem  as  well  as  the  related 
problems  of  track  acquisition  and  maintaining  track  over  time.  It  is  the 
characterization  and,  perhaps,  exploitation  of  this  "co-target  interference" 
phenomena  which  was  the  purpose  for  the  LDR  program. 

Cruise  missiles  may  be  launched  from  carriers  which  are  much  larger  in 
RCS.  Developing  feasible  processes  for  detecting  the  launch  by  the 
appearance  of  two  targets  where  only  one  originally  existed  is  vital  for  the 
surveillance  mission.  After  launch,  the  smaller  RCS  threats  travel  long 
distances  via  low  level  routes  in  the  presence  of  clutter  and  noise  sources,  are 
often  in  the  same  area  with  high  traffic  density,  and  their  paths  often  cross  with 
those  of  much  larger  targets,  all  of  which  may  cause  track  maintenance 
difficulties.  Additionally,  the  presence  of  jamming  and  other  Electronic  Counter¬ 
measures  (ECM)  will  cause  further  complications.  This  demanding 
environment  will  be  present  in  all  future  tactical  mission  scenarios. 
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2.0  REFERENCED  DOCUMENTS 


This  Large  Dynamic  Range  Final  Report  makes  references  to  the 
Software  User’s  Manual  and  the  Software  Test  Description, 

Also  the  following  references  are  cited  in  the  Report: 

1.  Kapogianis,  K.,  and  Macfarlane,  J.  F.,  Description  of  a  Multiple-Reflector 
Aircraft  Target  Model,  Hughes  Aircraft  Co.  internal  memo.,  Jan.  26, 1979, 

Ref.  2312.11/01. 

2.  Ruck,  G.  T.,  Barrick,  D.  E.,  Stuart,  W.  D.,  and  Krichbaum,  C.  K.,  Radar  Cross 
Section  Handbook,  Plenum  Press,  New  York  (1970). 

3.  Blackman,  S.  S.,  Multiple-Target  Tracking  with  Radar  Applications,  Artech 
House,  Norwood,  Mass.  (1986). 
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3.0  TECHNICAL  DISCUSSIONS 


3.1  PHENOMENOLOGY 

3.1.1  Introduction 

Two  targets  flying  in  close  proximity  present  a  challenge  for  radar  signal 
processing,  especially  if  one  of  the  targets  has  low  PCS.  When  the  two  are 
illuminated  by  an  electromagnetic  wavefront,  the  superposition  of  the  direct 
return  paths,  each  proportional  to  their  corresponding  PCS,  may  not  be 
appreciable  if  the  targets  are  in  close  proximity.  Mutual  target  coupling  can  now 
dominate  and  introduce  interesting  signal  behavior.  Figure  3.1. 1-1  shows  a 
possible  target  geometry  with  large  and  low  PCS  targets  with  wavefront 
multipath  present.  Because  the  large  target  will  have  a  complex  structure  with  a 
set  of  dominant  scattering  points,  many  bounce  paths  may  illuminate  the  low 
PCS  target  from  various  aspect  and  depression  angles.  This  can  in  turn 
enhance  the  small  monostatic  target  PCS  by  a  path  illuminating  a  favorable 
bistatic  angle. 


Rguro  3.1. 1-1.  Co-Target  Interference  May  Enhance  Small  Target  PCS 


To  understand  the  nature  of  the  mutual  target  interaction,  it  is  necessary 
to  understand  how  the  wavefront  bounces  off  each  target  at  different  bistatic  and 
monostatic  angles.  Thus  it  is  necessary  to  have  an  adequate  target  PCS  model 
that  allows  for  easy  and  reliable  predictions  of  target  PCS,  both  monostatic  and 
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bistatic.  Once  the  RCS  of  different  target  components  are  obtained  through  the 
RCS  model,  an  interaction  model  is  needed  to  describe  how  the  wavefront 
bouncing  from  one  target  is  received  by  the  other  target.  The  two  models  are 
combined  to  describe  and  quantify  the  target  interaction  signals  that  are 
eventually  detected  at  the  receiver. 


3.1.2  Target  RCS  Model 

As  a  part  of  the  interaction  prediction  model,  we  predict  target  RCS  using 
a  hybrid  model  that  is  a  combination  of  geometrical  optics  and  physical  optics. 
These  are  considered  to  be  high-frequency  techniques,  applicable  when  target 
dimensions  are  large  compared  with  wavelength.  In  the  LDR  simulation,  the 
radar  wavelengths  used  are  on  the  order  of  a  foot  or  shorter  so  most  airborne 
targets  will  be  at  least  10  wavelengths  long.  Therefore,  the  high-frequency 
modeling  is  justified.  At  X-band  (0.1  ft),  the  targets  will  be  hundreds  of 
wavelengths  long,  and  high  frequency  techniques  are  even  more  useful. 

It  must  be  pointed  out  that  the  electrical  size  requirement  actually  applies 
to  individual  scattering  features  and  not  the  overall  target  length.  This  was 
assumed  in  the  LDR  interaction  simulation  and  typically  assumed  in  practice 
because  high-frequency  techniques  can  only  be  applied  to  relatively  simple 
shapes  that  are  easily  described  in  mathematical  terms.  Therefore,  we  must 
break  up  the  actual  target,  whether  it  is  an  aircraft  or  a  missile,  into  a  collection 
of  approximate  simple  geometrical  shapes.  The  high-frequency  size 
requirement  then  applies  to  these  shapes,  and  not  necessarily  to  the  overall 
target.  Even  so,  most  target  features  are  still  within  the  high  frequency 
scattering  region.  They  should  be  at  least  five  wavelengths  in  size,  although 
reasonable  accurate  results  may  be  obtained  for  some  bodies  even  smaller 
than  this.^ 

At  the  lower  end  of  the  high-frequency  region  and  in  the  resonant 
scattering  region,  where  the  target  features  are  less  than  5  wavelengths,  high- 
frequency  techniques  lose  their  accuracies.  Other  prediction  techniques  such 
as  the  Method  of  Moments  (MOM)  should  be  used,  and  this  would  certainly 
enhance  the  capability  of  LDR  software  to  predict  target  RCS.  On  the  other 
hand,  it  must  be  kept  in  mind  that  the  scope  of  the  LDR  project  is  to  investigate 
2-target  bistatic  interaction,  it  is  not  a  comprehensive  target  RCS  modeling 
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project.  Refined  high-frequency  prediction  techniques  such  as  geometric  theory 
of  diffraction  (GTD),  method  of  equivalent  currents  (MEC),  and  physical  theory  of 
diffraction  (PTD)  have  to  be  left  out  of  the  LDR  model.  These  techniques 
demand  a  great  deal  of  mathematical  and  software  complexity  that  can  not  be 
provided  within  the  scope  of  this  study. 

In  the  high-frequency  region,  collective  interactions  are  very  small,  so 
that  a  body  can  be  treated  as  a  collection  of  independent  scattering  centers. 
Each  part  of  the  body  scatters  energy  essentially  independently  of  all  other 
parts.  The  fields  induced  on  a  portion  of  the  target  are  only  due  to  the  incident 
wave  and  not  the  energy  scattered  by  other  parts.2  The  simplest  high  frequency 
method  is  geometric  optics  (GO).  The  RCS  is  given  by  a  simple  formula  that 
involves  only  the  local  radii  of  curvature  at  the  specular  point.  When  the  radii  of 
curvature  becomes  infinite,  as  in  the  case  of  a  cylinder  of  fiat  plate,  then  GO  fails 
and  we  would  then  have  to  use  physical  optics  (PO)  instead.  Physical  optics 
give  good  results  if  the  surface  is  not  too  small  and  if  the  scattering  direction 
does  not  swing  too  far  from  the  specular  direction.  At  wide  angles  from 
specular,  however,  physical  optics  fails.  This  is  because  the  contributions  from 
the  edges  are  ignored,  which  can  be  accounted  for  by  the  GTD  method.^ 

The  LDR  interaction  model  uses  geometric  optics  to  predict  the  RCS  of 
the  generic  aircraft  and  physical  optics  for  the  cruise  missile.  This  is  because 
the  aircraft  can  be  modeled  as  a  collection  of  ellipsoids  that  are  several 
wavelengths  in  size.  With  an  ellipsoid,  the  easiest  way  to  compute  the  RCS  is 
to  use  geometric  optics  formulas  that  are  easily  available.  The  cruise  missile, 
being  smaller  than  the  aircraft  and  being  more  "flat",  is  modeled  as  a  collection 
of  other  geometric  shapes.  The  main  component,  the  fuselage,  is  modeled  as  a 
hemispherically-capped  cylinder.  The  rest  of  the  missile  is  modeled  using  flat 
plates  and  wires.  For  these  shapes,  the  RCS  formulas  are  derived  using 
physical  optics  theory. 

3.1. 2.1  Aircraft  RCS  Modfll 

When  the  aircraft  Is  illuminated  by  electromagnetic  energy,  it  has  been 
observed  experimentally^  that  most  of  the  energy  is  reflected  from  a  discrete 
number  of  locations.  Each  location  can  be  modeled  as  a  single  scattering  point, 
and  thus  a  finite  number  of  scatterers  can  be  used  to  model  the  entire  aircraft. 
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In  practice,  the  major  sections  of  the  aircraft  are  modeled  as  basic  geometrical 
shapes  such  as  ellipsoids.  Each  ellipsoid  has  one  reflection  point,  referred  to 
as  the  "specular  point",  which  is  determined  by  geometric  optics.  A  good 
approximation  of  the  aircraft's  radar  cross  section  can  then  be  estimated  as  a 
function  of  aspect  angle. 

The  LDR  interaction  simulation  employs  this  model  as  developed  in 
Reference  1.  A  Cessna  A-37B  aircraft  was  used  although  any  aircraft  can  be 
modeled  when  its  actual  physical  dimensionings  are  specified. 


3.1.2.1.1  Coordinate  System 

An  illustration  of  the  ellipsoids  is  shown  in  Figure  3.1. 2.1. 1-1.  The 
coordinate  system  X.J.,  Y^,  and  Zj  is  the  reference  coordinate  system  for  the 

aircraft  RCS  computation.  The  origin  is  the  center  of  the  fuselage  ellipsoid  with 
the  Xj  axis  always  pointing  in  the  nose  direction  of  the  aircraft.  Each 


component  ellipsoid  has  its  own  local  coordinate  system  generated  by  rotating 
the  reference  coordinate  system  by  the  angles  a,  p,  and  y,  and  then  translated 
by  X,  Y,  and  Z,  A  local  coordinate  system  is  shown  in  Figure  3.1. 2.1. 1-2.  The 
ellipsoids  and  their  related  parameters  for  the  Cessna  are  listed  in  Table 
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Figure  3. 1.2. 1.1-1.  Ellipsoid  Components  of  Aircraft 
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Rgure  3.1 .2.1. 1-2.  Local  coordinate  system  of  an  ellipsoid 

To  transform  from  the  reference  coordinate  system  to  a  local  coordinate 
system  we  use  the  transformation  matrix  T: 


Xl 

cosa  cosp 

sina  cosp 

-sinp 

Yl 

cosa  sinpsiny-sina  cosy 

sina  sinpsinyfcosa  cosy 

sinycosp 

• 

Zl 

cosa  sinpcosyfsina  siny 

m 

sina  sinpcosy-cosa  siny 

cosycosp 

or  more  compactly. 

Xl=  t(xt-xo)^ 

where  Xq  is  the  displacement  vector  composed  of  X,  Y,  and  Z. 

The  radar  location  is  specified  in  terms  of  the  azimuth  and  elevation 
angles  as  measured  in  the  reference  coordinate  system.  An  illustration  is 
shown  in  Rgure  3. 1.2. 1.1 -3.  Thus,  nose-on  aspect  is  0“  azimuth  and  0° 
elevation  while  tail  aspect  is  ISO**  azimuth  and  O'*  elevation.  The  line  of  sight  to 
the  radar  is  a  vector  with  reference  coordinates  given  by; 


■Rx 

cos  az  cos  el 

Ry 

sin  3z  cos  ei 

Rz 

Sind 
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Figure  3.1 .2.1 .1-3.  Reference  coordinate  system 


3.1. 2.1. 2  Overview  of  Geometric  Optics 

Some  basic  principles  of  geometric  optics  are  presented  here  for 
reference.  This  technique  is  a  ray-tracing  procedure  whereby  the  energy 
incident  upon  the  surface  of  the  object  travels  in  ray  bundles.  The  reflection 
coefficient  for  a  perfectly  conducting  surface  is  -1,  implying  no  reduction  of  the 
intensity  of  the  reflected  wave.  This  is  only  true  at  the  point  of  reflection  (the 
specular  point),  and  as  is  noted  above,  the  ellipsoids  are  modeled  in  terms  of 
their  specular  points.  Because  of  this,  the  decay  in  intensity  of  the  reflected 
wave  does  not  come  into  play  when  computing  the  ellipsoid  RCS. 

The  polarization  of  the  incident  wave  is  also  neglected.  The  RCS  of  the 
ellipsoid  is  the  same  regardless  of  the  type  of  incident  polarization. 

Geometric  optics  is  not  valid  in  the  neighborhood  of  sharp  points  or 
boundaries,  such  as  the  tip  of  a  cone,  the  edge  of  a  wedge,  etc.  For  these 
cases,  other  techniques  must  be  used,  such  as  GTD  or  PTD.  As  mentioned 
previously,  the  Large  Dynamic  Range  project  does  not  use  these  more 
sophisticated  theories. 

The  results  of  geometric  optics  give  the  RCS  in  terms  of  the  radii  of 
curvature  of  the  body  at  the  specular  point.  The  formula  is  given  by: 

a  =  ^1^2, 
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TABLE  3.1 .2.1 .1-1 .  PARAMETERS  FOR  THE  AIRCRAFT  ELLIPSOIDS 
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where  ai  and  aa  are  the  principal  radii  of  curvature  of  the  body  at  the  specular 
point.  From  this  expression,  it  can  be  seen  that  the  RCS  is  not  dependent  upon 
the  frequency  of  the  incident  wave.  Also,  this  formula  is  valid  for  perfectly 
conducting  bodies.  The  difficulty  with  this  formula  is  that  for  flat  or  singly  curved 
surfaces  (not  quadric),  one  or  both  radii  of  curvature  are  infinite,  giving  infinite 
RCS.  Physical  optics  is  needed  to  overcome  this  dilemma. 

The  geometric  optics  approach  is  valid  for  wavelengths  that  are  small 
compared  to  the  dimensions  of  the  objects.  In  the  LDR  simulation,  we  use 
wavelengths  that  are  1  foot  or  less.  Thus  the  model  is  more  accurate  at  the 
higher  ends  of  the  spectrum,  such  as  above  5  GHz  (0.2  ft)  or  so.  Some  of  the 
smaller  ellipsoids  used  in  the  model  of  the  aircraft  have  dimensions  of  this  order 
(0.2  ft)  so  their  RCS  values  are  not  as  accurate  as  the  larger  ellipsoids.  Ideally, 
the  objects  should  be  at  least  5  wavelengths  in  size  although  as  mentioned 
above,  reasonably  accurate  results  can  still  be  obtained  for  smaller  objects.  It  is 
also  necessary  to  have  objects  that  are  sufficiently  smooth.  The  ellipsoids  in  the 
LDR  simulation  certainly  satisfy  these  conditions. 

For  longer  wavelengths  more  sophisticated  RCS  computational  methods 
must  be  used  such  as  method  of  moments,  integral  equations,  etc.  However, 
the  LDR  simulation  is  not  expected  to  simulate  targets  at  below  L-band  (1  GHz 
or  1  ft). 


3.1. 2.1. 3  RCS  Computations  for  Ellipsoids 

The  RCS  of  an  ellipsoid  as  given  by  geometric  optics  is:® 


a  =  7C 


2.  2  2 
a  b  c 


where 

p  =  a^x  +  b^y  +  c^z 

a,  b,  c  are  the  ellipsoid  semi-axes,  and 

Dx,  Dy,  and  D^  are  the  components  of  the  local  radar  line  of  sight. 
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Geometric  optics  assumes  that  each  ellipsoid  is  represented  by  only  one  point, 
the  specular  point,  where  most  of  the  electromagnetic  energy  is  reflected  back. 
The  local  coordinates  of  this  point  is  given  by  the  vector  q  with  components: 


Qx 


Vp  . 


An  illustration  of  the  specular  point  geometry  and  the  radar  scenario  can 
be  seen  in  Figure  3. 1.2. 1.1-1. 


3.1. 2.1. 4  Monostatic  RCS  of  Aircraft 

When  computing  the  scattered  power  from  a  complex  object  in  the  high- 
frequency  region,  the  fields  scattered  by  the  many  different  components  must  be 
calculated  and  then  added  together  before  squaring  to  obtain  the  scattered 
power.  Interference  effects  may  thus  be  represented  because  the  phase 
relationship  between  the  various  scatterers  are  preserved.  Since  the  scattered 
power  is  proportional  to  the  RCS,  this  implies  that  the  square  root  of  the  RCS  of 
the  various  aircraft  components  are  algebraically  summed  with  appropriate 
phases  to  obtain  the  total  RCS. 

The  total  monostatic  RCS  of  the  aircraft  can  be  computed  by  summing 
the  square  root  of  the  RCS  of  ail  ellipsoids,  it  is  given  by: 


Ot  = 


,2 


N  1/2 

X(aj)  exp(i<|)j) 
i*i 


where  Oj  is  the  RCS  of  the  ]-th  ellipsoid  and  <t»j  is  the  relative  phase  angle  of  the 

j-th  ellipsoid.  The  magnitude  of  <t)j  is  determined  by  selecting  a  reference  plane 
normal  to  the  direction  of  incidence  and  measuring  the  distance  from  this  plane 
to  the  specular  point  of  the  ellipsoid.  The  reference  plane  is  chosen  to  pass 
through  the  reference  coordinate  origin.  This  plane  and  the  direction  of 
incidence  as  the  X-axis  form  what  is  called  the  reference  line-of-sight  (LOS) 
coordinate  system.  To  go  from  the  reference  coordinate  system  to  this  system 
we  rotate  by  a  transformation  matrix  S: 
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‘los 


=  SXt 


Therefore,  the  specular  point  q  in  the  reference  LOS  coordinate  system  is 
given  by; 


‘los 


=  s(t^q+; 


with  being  the  transpose  of  T.  The  first  component  of  Q|o,  is  the  distance 
from  the  reference  plane  to  the  specular  point  of  an  ellipsoid.  Call  this  dj  for 
the  j-th  ellipsoid.  Then  the  phase  is  given  by 


«  2di 
— ! 

X  . 


The  factor  2  is  to  account  for  the  incident  and  reflecting  path. 

3.1. 2.1. 5  Shadowing  Effects 

Because  the  geometric  optics  model  assumes  that  for  each  ellipsoid 
most  of  the  electromagnetic  energy  is  reflected  from  a  single  specular  point,  it 
follows  that  there  will  be  instances  where  the  specular  point  is  obstructed  by 
another  ellipsoid.  This  is  called  "shadowing".  An  ellipsoid  is  assumed  to  not 
contribute  to  the  total  aircraft  RCS  when  it  is  shadowed.  Large  sections  of  the 
ellipsoid  may  be  visible,  but  its  RCS  will  be  set  to  zero  if  its  specular  point  is 
shadowed.  In  reality,  a  small  amount  of  energy  is  scattered  from  these  sections 
because  they  are  accessible  to  the  incident  plane  wave.  However,  for  the 
purpose  of  the  LDR  simulation,  this  energy  can  be  neglected  since  it  contributes 
little  to  the  total  scattered  energy  from  the  ellipsoid.  Most  of  the  scattered 
energy  would  have  come  from  the  unobstructed  specular  point,  which  in  this 
case  is  blocked. 

Another  type  of  shadowing,  "layered  shadowing"  occurs  when  part  of  the 
ellipsoid  model  is  nonexistent.  Sometimes  it  is  necessary  to  model  certain 
target  components  as  partial  ellipsoids  to  represent  the  true  target  as  closely  as 
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possible.  For  example,  the  wings  of  the  aircraft  are  modeled  as  ellipsoids  but 
they  are  actually  half-ellipsoids.  Thus  if  the  specular  point  of  a  wing,  when  seen 
from  a  particular  aspect  angle,  happens  to  fall  in  the  non-existent  portion  of  the 
ellipsoid  then  the  wing  is  shadowed,  and  its  RCS  is  set  to  zero.  Another 
example  is  the  tailplane  where  the  central  portion  of  the  ellipsoid  does  not  exist, 
but  is  actually  considered  a  part  of  the  tailfin.  Other  regions  where  layered 
shadowing  occurs  are  listed  in  Table  3.1. 2.1. 5-1. 

The  components  that  are  partial  ellipsoids  have  reduced  RCS.  If  the 
specular  point  of  such  an  ellipsoid  fails  in  an  actual  portion  of  the  ellipsoid  (i.e. 
layered  shadowing  does  not  occur),  and  the  first  type  of  shadowing  does  not 
occur,  then  the  RCS  of  the  ellipsoid  is  computed  from  the  ellipsoid  RCS  formula 
above  using  dimensions  of  the  full  ellipsoid.  This  value,  however,  is  then 
reduced  in  the  simulation  by  the  corresponding  fraction  of  the  ellipsoid  that  is 
real.  For  example,  the  RCS  of  the  wings  of  the  airplane  are  half  of  what  would 
be  computed  using  the  formula. 

The  actual  implementation  of  the  shadowing  algorithm  begins  by  sorting 
the  specular  points  in  order  of  their  distances  from  the  LCS  reference  plane.  In 
other  words,  the  ellipsoids  with  larger  dj  are  considered  to  be  in  front  of  other 
ellipsoids  with  respect  to  the  oncoming  wavefront.  Suppose  that  the  i-th 
ellipsoid  is  in  front  of  the  k-th  ellipsoid,  then  the  k-th  ellipsoid  is  shadowed  if  its 
specular  point  satisfies  the  equation: 

2  .2  2  ” 
a  b  c  , 

for  some  ^>0,  where 

dxt  3rtd  Pj  =  the  components  of  the  k-th  specular  point  in  the  i-th 
locai  coordinate  system, 

Lj,,  Ly,  and  =  the  components  of  the  radar  LCS  in  the  i-th  local 
coordinate  system,  and 

a.  b,  and  c  are  the  semi-axes  of  the  i-th  ellipsoid. 

An  illustration  of  this  equation  is  shown  in  Figure  3.1. 2.1. 5-1. 
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- - -  Wave  front 

‘Spec'jiaf  pcint 

^  ^ z) 


.  ■  7~fr  -  . 


V.'avefrcr.t 


Figure  3. 1.2. 1.5-1.  Shaccwing  Geometry 

Any  ellipsoidal  component  which  is  shadcwed  will  have  its  RCS  {a 
to  zero  in  the  equation  for  a-j  above. 


TABLE  3.1. 2.1. 5-1.  LAYERED  SHADOWING  REGIONS 


Fuselage 


Filets  Helmet 


Tailpiane 


Left  Wing 


Rigtit  'Mrg 


(SO-  <  4  s  270*  and  9  <  145=] 
and 

-0.267  <0111  <0.257  ft 


1.  Xy(3)  <  0.0 

2.  3.20  <  Xt,;31  <  3.34  1 


Xji2)  <  2.S9  ft 


Xn2)  >  -2.39  ft 


Left  Seat 


Hi^t  Seat 


Outer  Lett  Wing 
Ho  Tank 


Cuter  Right  Wing 
Tio  Tank 


C<i),  02).  and  Gf3)  am  (Tie  ioc3i  oorrsonent  coorcinates  of  soecjta/'  potnt 
XT(n.XTt21,Xr(3)  are  in  Oh®  aircrart  .'efarenca  coorcinates. 


3.1. 2.1. 6  Bistatic  RCS  of  Aircraft 

In  the  high-frequency  region  it  can  be  shown  that  most  of  the  scattered 

field  comes  from  the  immediate  neighborhood  of  the  specular  point  even  in 

bistatic  scattering®.  Figure  3.1. 2.1. 6-1  illustrates  the  geometry  of  bistatic 

^  ^  ^  ^ 

scattering.  For  incident  direction  k|  and  scattering  direction  ks.kjjiks,  the 
specular  point  is  that  point  whose  normal  vector  n  bisects  the  angle  between 
kj  and  k*  and  lies  in  the  plane  formed  by  these  vectors.  The  specular  point  can 
be  found  by  assuming  that  an  equivalent  monostatic  radar  line-of-sight  is  along 
the  vector  n  and  applying  the  formulas  above.  This  Is  also  the  basis  for 
computing  bIstatic  RCS  from  monostatic  RCS  using  what  is  know  as  the 
monostatic/bistatic  equivalence  theorem.  This  theorem  states  that  as  the 

wavelength  tends  to  zero,  the  bistatic  RCS  for  transmitter  (incident)  direction  kj 
and  receiver  (scattering)  direction  k*.  kj;6ks ,  is  equal  to  the  monostatic  RCS  for 

the  transmitter-receiver  direction  kj  +  k,  for  bodies  that  are  sufficiently  smooth. 
Therefore,  we  can  determine  the  bistatic  RCS  of  the  aircraft  by  applying  the 
monostatic  results  in  the  direction  kj  +  kj. 


Rgure  3.1 .2.1 .6-1 .  Bistatic  scattering  geometry 

In  computing  the  total  bistatic  RCS  of  the  aircraft,  the  relative  phase  of 
each  ellipsoid  is  computed  in  a  slightly  different  way  than  that  of  the  monostatic 
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case.  We  compute  the  distance  of  each  specular  point  as  measured  from  two 
reference  planes.  One  plane  is  normal  to  the  transmitter  direction  kj  while  the 

^  (t)  (r) 

other  is  normal  to  the  receiver  direction  kg.  Call  these  distances  dj  and  dj 
for  the  j-th  ellipsoid.  Then  the  phase  is  given  by 


<l>i  =  2k 


The  monostatic/bistatic  equivalence  theorem  has  certain  limitations  that 
limit  its  usefulness.  First  of  all,  it  is  not  accurate  for  bistatic  angles  greater  than 

about  135°.  The  bistatic  angle  is  the  angle  between  kjand  kg .  Furthermore,  its 
accuracy  decreases  as  the  radar  frequency  is  lowered.  Despite  these 
weaknesses,  the  theorem  is  used  for  bistatic  scattering  from  the  airplane 
because  it  is  the  most  accessible  method  for  computing  bistatic  PCS  with 
ellipsoids  and  geometric  optics.  However,  the  theorem  is  not  used  to  compute 
bistatic  scattering  from  the  missile. 


3.1. 2.1. 7  Bistatic  Shadowing 

Shadowing  in  the  bistatic  case  is  similar  to  the  monostatic  case  but  is 
more  complex.  Even  though  the  direction  kj  +  k,  is  used  to  compute  the 
equivalent  bistatic  PCS,  the  actual  transmitting  and  receiving  directions  of  the 
electromagnetic  energy  must  be  used  to  check  for  shadowing.  If  the  specular 
point  of  an  ellipsoid  is  obstructed  when  seen  from  either  the  transmitting  or  the 
receiving  direction,  then  a  shadowed  condition  for  that  ellipsoid  must  exist. 
Thus  it  is  necessary  to  check  for  shadowing  with  respect  to  two  directions  in  the 
bistatic  case. 

Bistatic  shadowing  also  differs  from  monostatic  shadowing  in  some 
interesting  ways.  First,  a  specular  point  can  be  in  the  clear  with  respect  to  one 
LOS,  but  is  shadowed  with  respect  to  the  other  LOS.  In  this  case,  the  net  effect 
is  that  the  ellipsoid  is  shadowed.  Second,  in  the  bistatic  case,  when  the 
specular  point  of  an  ellipsoid  is  in  front  of  the  specular  point  of  another  ellipsoid, 
as  seen  from  a  given  direction,  it  can  still  be  shadowed  by  the  ellipsoid  that  is 
behind.  This  means  that  we  must  be  careful  to  check  if  "front"  components  are 
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shadowed  by  "rear"  components.  In  the  monostatic  case,  components  in  front 
are  never  shadowed  by  components  behind  them.  Figures  3.1.2.1.7-1  (a)  &  (b) 
contain  some  examples  of  bistatic  shadowing. 


XMfT  I 

Specular  Point  of  B 


B  shadowed  by  C  as  seen 
from  RCV  direction 


RCV 


Figure  3.1. 2.1. 7-1  (a).  Either  line-of-sight  can  shadow  ellipsoid 


Figure  3.1 .2. 1.7-1  (b).  Another  type  of  bistatic  shadowing 
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3.1. 2.1. 8  Computed  RCS  of  Airplane 

In  this  section  we  present  the  results  for  the  computed  RCS  of  the 
airplane.  Figures  3.1. 2.1. 8-1  (a)  &  (b)  show  the  monostatic  RCS  of  the  aircraft 
as  a  function  of  the  azimuth  angle  for  L-band  and  X-band.  Both  figures  are  for 
the  case  of  a  radar  in  the  0°  elevation  plane.  The  most  striking  features  of  the 
plots  are  the  peaks  at  broadside  (SO*’).  This  is  not  surprising  since  at  broadside 
the  fuselage  of  the  aircraft  presents  itself  most  completely  to  the  incident  wave. 
The  peaks  are  about  35  dBsm  in  magnitude  while  the  average  RCS  in  both 
cases  is  about  0  to  5  dBsm.  The  X-band  curve  fluctuates  much  more  than  the  L- 
band  curve,  and  this  is  due  to  the  much  smaller  wavelength,  thus  increasing  the 
phase  difference  from  aspect  to  aspect. 


AZIMUTH  ;C=-2., 

Figure  3. 1.2. 1.8-1  (a).  Monostatic  RCS  of  aircraft  @  0°  elevation 

L-band 
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Rgure  3. 1.2. 1.8-1  (b).  Monostatic  RCS  of  aircraft  @  0*  elevation 

X-band 

It  is  interesting  to  look  at  the  individual  component  RCS  for  this  geometry. 
Figures  3. 1.2. 1.8-2  (a)  &  (b)  present  the  RCS  of  the  ten  ellipsoids  making  up 
the  aircraft  as  functions  of  azimuth.  At  90®,  the  main  contributors  to  the  aircraft 
RCS  are  the  tailfin,  the  fuselage,  and  the  left  tank.  Surprisingly,  the  tailfin  is  the 
main  component  responsible  for  the  peak  at  90®.  The  discontinuities  in  the 
curves,  where  a  component  RCS  suddenly  drops  off  to  zero,  are  the  aspects 
where  that  component  is  shadowed.  For  example,  the  helmet  is  visible  all  the 
way  out  to  1 40®  azimuth,  then  it  is  shadowed. 

Next,  Rgures  3.1. 2.1. 8-3  (a)  &  (b)  show  the  aircraft  RCS  for  the  case 
where  the  radar  is  pointing  down  at  the  aircraft  with  an  elevation  angle  of  -45®. 
We  Immediately  see  that  the  strong  sharp  peaks  at  broadside  have  been 
replaced  by  weaker  broader  peaks.  A  look  at  the  component  RCS  plots  for  this 
elevation  plane  in  Figures  3. 1.2.1. 8-4  (a)  &  (b)  shows  why.  The  tailfin  is  no 
longer  a  big  contributor  to  the  aircraft  RCS;  its  RCS  has  now  fallen  way  down  to 
-30  dBsm  at  90®.  The  fuselage  is  now  the  biggest  contributor  but  at  90®  it  is 
shadowed,  as  shown  by  the  sudden  "gap"  in  RCS.  The  other  main  contributors 
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Figure  3.1 .2.1 .8-2  (a).  PCS  of  aircraft  components  @  0®  elevation 


Figure  3.1 .2.1 .8-2  (b).  PCS  of  aircraft  components  @  0°  elevation 
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Rgure  3.1 .2.1 .8-3  (b).  Monostatic  RCS  of  aircraft  @  45®  elevation  -  X  band. 
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Figure  3. 1,2. 1.8-4  (a).  RCS  of  aircraft  components  @  45^  elevation 
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Figure  3. 1.2. 1.8-4  (b).  RCS  of  airc.'aft  ccmpcnents  @  45’  elevation 
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Next,  we  show  the  monostatic  RCS  of  the  aircraft  as  a  function  of  varying 
elevation  angle  in  Figure  3. 1.2. 1.8-5.  The  azimuth  angle  of  the  radar  is  0°,  or 
nose  aspect.  Negative  elevation  angles  mean  that  the  radar  is  underneath  and 
locking  up  at  the  aircraft.  Positive  elevation  is  radar  looking  down.  The  figure 
shows  that  between  -10®  and  -80®,  the  aircraft  RCS  is  weak,  in  comparison  with 
10®  to  80®.  Between  0®  and  -10®  the  aircraft  RCS  is  essentially  constant  at 
around  5  dBsm.  This  is  useful  to  know  in  the  case  of  a  ground  radar  locking  up 
at  the  airplane  from  a  distance. 


E'_£VAT!CM  ;C£G.) 

Rgure  3. 1.2. 1.8-5.  Monostatic  RCS  of  aircraft  @  0®  azimuth 

L-band 

Finally,  the  bistatic  RCS  of  the  aircraft  is  shewn  in  Figure  3. 1.2.1. 8-6. 
The  transmitter  and  the  receiver  are  always  set  at  150®  apart  in  azimuth,  and  the 
angle  plotted  is  the  varying  azimuth  angle  of  the  transmitter.  In  comparing  this 
plct  with  the  plot  in  Figure  3. 1.2. 1.8-3  (a)  above,  it  is  interesting  to  note  that  the 
RCS  obtained  bistatically  in  this  way  is  generally  higher  than  the  monostatic 
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RCS.  In  certain  situations  this  would  be  very  helpful.  It  may  be  possible  to 
exploit  bistatic  scattering  when  we  discuss  the  target  interactions  later. 


L-BAND  (1.3  GHZ) 

XMTR/RCVR  150  DEG.  APART 


Figure  3. 1.2. 1.8-6.  Bistatic  RCS  of  aircraft  @  45®  elevation 

L-band 


3. 1.2.2  Cruise  Mi.s.Rile  RCS  ModPi 

Like  the  aircraft,  the  cruise  missile  is  mcdeled  as  a  collection  of  discrete 
components.  Unlike  the  aircraft,  however,  the  missile  includes  non-spherical 
shapes  such  as  a  cylinder,  flat  plates,  and  wires.  As  noted  before,  the  smaller 
and  "flat"  shape  of  the  cruise  missile  dees  not  lend  itself  well  to  ellipsoid 
modeling,  and  other  geometric  objects  must  be  used. 

The  missile  fuselage  is  the  main  contributor  to  missile  RCS.  It  is  modeled 
as  a  hemispherically-capped  cylinder.  This  means  that  the  central  portion  of 
the  fuselage  is  a  right  circular  cylinder  with  flat  ends  while  the  rounded  ends  are 
modeled  as  two  hemispheres  (half-spheres).  The  wings,  tailfin,  and  tailplane 
are  modeled  as  a  collection  of  triangular  flat  plates  while  the  straight  edges  of 
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these  parts  are  modeled  as  thin  wires.  Other  parts  of  the  missile  such  as  air 
ducts,  dams,  propellers,  etc.  are  not  modeled.  Different  views  of  the  generic 
cruise  missile  are  shown  in  Rgure  3.1. 2.2-1. 

The  RCS  of  the  non-spherical  shapes  are  computed  using  the  physical 
optics  technique.  Since  geometric  optics  fails  for  shapes  such  as  cylinder,  and 
flat  plates,  we  must  turn  to  physical  optics  instead  as  the  other  high-frequency 
technique  that  is  easily  implemented. 


PUN  VIEW 


SIDE  VIEW 


DERIVED  FROM  DETAILED  SURFACE 
DESCRIPTION  PROVIDED  BY  BOEING 
997  POINTS 
1690  FACETS 


Figure  3.1 .2.2-1 .  Views  of  generic  cruise  missile 


3.1. 2.2.1  Coordinate  System 

The  reference  coordinate  system  for  the  missile  is  oriented  in  the  same 

way  as  that  of  the  aircraft  (shown  in  Figure  3.1. 2.1. 1-3).  The  origin  is  the  center 
of  the  missile  fuselage  cylinder  with  the  X-axis  pointing  in  the  direction  of  the 

nose  of  the  missile.  Each  missile  component  has  its  own  local  coordinate 
system  generated  by  rotating  the  reference  coordinate  system  by  the  angles  a. 
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P,  and  Y,  and  then  translated  by  X,  Y,  and  Z.  For  each  type  of  geometric  shape 
the  local  coordinate  system  is  shown; 


j-^— b— ►! 


Yl 
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Rectangular  plate 


A 


H^b-H 


The  wire  does  not  have  a  local  coordinate  system.  We  simply  look 
at  its  dimension  along  the  missile  reference  coordinate  system. 

The  various  missile  components  and  their  related  parameters  are  listed  in  Table 

3.1.2.2.1-1, 

As  with  the  aircraft,  to  transform  from  the  reference  coordinate  system  to  a 
local  coordinate  system  we  use  the  transformation  matrix  T; 

t(xt-Xd)^ 


where  Xq  is  the  displacement  vector  of  the  missile  component  from  the  origin 
and  is  composed  of  X,  Y,  and  Z.  For  most  of  the  missile  components,  the  matrix 
T  is  exactly  the  same  matrix  as  that  given  in  the  aircraft  section  above. 
However,  the  triangle  plates  making  up  the  tailfin  of  the  missile  are  treated 
differently.  The  orientation  of  these  plates  cannot  be  adequately  represented  by 
the  three  rotation  angles  a,  p,  and  y.  More  angles  are  needed  which  change 
the  transformation  matrix  for  these  components.  This  alternative  matrix  is  given 
in  the  Appendix. 
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TABLE  3.1.2.2.1-1.  COMPONENTS  OF  THE  CRUISE  MISSILE 
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3.1. 2.2.2  Overview  of  Physical  Ootics 

Physical  optics  overcomes  the  difficulty  of  infinite  radii  of  curvature  of  flat 
or  singly  curved  surfaces  in  computing  RCS.  It  does  this  by  approximating  the 
induced  surface  fields  and  integrating  them  to  obtain  the  scattered  field.  In  the 
general  case  of  scattering  from  an  arbitrary  body,  the  scattered  fields  are  given 
by  the  Stratton-Chu  integral  equations.^  They  are  extremely  difficult  to  solve, 
however,  in  the  general  case,  and  many  simplifying  assumptions  can  be  made. 
In  solving  the  Stratton-Chu  integrals,  physical  optics  makes  these  assumptions: 

1 .  The  incident  field  is  assumed  to  be  time  harmonic  and  a  plane  wave 
over  the  entire  finite  scattering  object. 

2.  The  total  fields  at  the  surface  on  the  "shadowed  side"  of  the  object  are 
zero.  This  assumption  is  valid  only  when  the  wavelength  is  small  in 
comparison  with  the  dimensions  of  the  object. 

3.  The  observation  point  (or  receiving  antenna)  is  far  from  the  object  in 
terms  of  wavelength  and  object  dimensions:  r  >  2D2/X,  where  r  is  the 
distance  from  object  to  the  observation  point,  D  is  the  largest  object 
dimension,  and  X  is  the  wavelength. 

4.  Tannpnt-Planft  Aonroximation.  In  this  approximation,  the  surface  at  and 
near  any  point  is  assumed  to  be  an  infinite  plane  tangent  to  the  surface 
at  that  point.  This  leads  to  the  following  expressions  for  the  total 
surface  fields  for  a  perfectly  conducting  surface: 

^  T 
nxE  =0 

T  ^  i 

nxH  =2nxH 

I 

where  and  are  the  total  surface  E  and  H  fields,  respectively,  H ' 
is  the  incident  H  field,  and  n  is  the  normal  vector  to  the  surface. 

These  assumptions  simplify  the  Stratton-Chu  integrals  tremendously  and  lead 

s 

to  the  physical  optics  integral  for  the  scattered  H  field: 

H®(P)  =  exp(ikRo)  f  ^  ^  H*)  x  k^exp(-ikok  Tjds 

2  X  Rq  Js, 

where 
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P  is  the  far-field  observation  point, 

Rq  is  the  distance  from  the  origin  to  P, 
k  =  27:/X. 

$1  is  the  surface, 

k  is  the  unit  vector  in  the  scattered  field  direction  of  propagation,  and 
r  is  the  vector  from  the  origin  to  a  surface  patch  ds’. 


When  the  surface  such  as  a  flat  plate  lies  on  the  xy-plane,  the  integral  can  be 
written  in  terms  of  x  and  y  and  simplified  to: 


H  (P)  = 


2  jr  Rq 


fT 

k  X 

HqX  z 

[1  exp 

ik  (  k  -  k  )  •  r 

dxdy 
(Eqn.  3.1.2.2.2-1) 


where 


k  is  the  unit  vector  in  the  incident  field  direction  of  propagation, 

Hq  is  found  in  H*  =  Hoexp  ik'  t  ,  and  He  is  a  constant  vector 
perpendicular  to  the  direction  of  incidence, 
z  is  the  unit  vector  in  z. 


This  is  the  equation  we  later  use  to  derive  the  bistatic  scattering  formula  for  the 
triangular  plates. 


3. 1.2. 2.3  Cylinders 

While  only  ellipsoids  were  used  to  model  components  of  the  aircraft,  the 
central  portion  of  the  cruise  missile  is  modeled  as  a  hemispherically-capped 
cylinder.  Thus,  it  was  necessary  to  develop  an  RCS  model  for  the  cylinder. 

The  scattering  geometry  and  symbols  for  the  cylinder  are  shown  in 
Figure  3. 1.2. 2.3-1.  The  plane  of  incidence  is  defined  as  the  plane  with  the 
incident  wave  propagation  direction  and  the  cylinder  axis  (the  x-axis  here),  and 
is  the  x-2  plane  in  the  figure  (although  not  in  general).  The  incident  direction 

makes  an  angle  V,  as  shown  with  the  y-z  plane.  The  plane  of  scattering  is 
defined  as  the  plane  with  the  scattering  propagation  direction  and  the  cylinder 
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axis,  and  is  shown  rotated  by  azimuth  angle  (])'  from  the  z  axis.  The  scattering 
direction  makes  an  angle  ¥s  as  shown  with  the  y*z  plane.  As  shown,  the  two 

angles  ¥i  and  ¥s  have  the  same  sign.  There  are  two  polarization  modes  of 
interest.  The  TM  mode  is  when  the  E-field  is  in  the  plane  of  incidence  and 
scattering  (H-field  perpendicular).  The  TE  mode  is  when  the  E-field  is 
perpendicular  to  the  plane  of  incidence  and  scattering. 


X 


Figure  3.1 .2. 2.3-1 .  Scattering  geometry  for  circular  cylinder 


In  general,  a  scattered  field  from  a  cylinder  of  finite  length  propagates  in 
all  directions,  so  that  one  can  consider  the  scattered  field  at  arbitrary  angle  ¥s. 
However,  for  cylinders  whose  length  is  very  much  greater  than  the  radius,  the 
scattered  field  at  and  near  Vg  =  -¥i  is  much  larger  than  the  scattered  field  in 
other  directions.  This  is  specular  reflection.  For  the  cruise  missile,  the  length  is 
much  larger  than  the  radius  so  this  statement  is  true.  The  angle  <|)'  is  arbitrary 
so  the  propagation  vectors  of  the  scattered  field  from  the  cylinder  form  a  cone  in 

the  forward  direction  which  has  a  half-angle  of  k  /  2  -  \|/g  with  respect  to  the  x 
axis  as  shown. 
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Given  the  range  of  radar  frequency  we  will  be  dealing  with  (L  to  X  band) 
and  the  radius  of  the  missile  cylinder  (0.83  ft),  there  are  two  cases  of  interest 
when  computing  the  RCS  of  the  cylinder. 

Case  1.  Cylinder  radius  of  same  order  as  wavelength  (0.5  <  kga  <  1 0  ) 

In  this  case  the  scattered  field  becomes  more  concentrated  in  the 
specular  direction  \j/g  =  -vKj  because  the  cylinder  length  (20  ft)  is  large 
compared  to  the  radius.  Using  the  techniques  in  Reference  2,  pp.302-304,  the 
bistatic  scattering  cross  section  is  given  by 


4  k  h  2 

a  (Vi  ({)’ )  = - cos  Vs  (¥i.4>') 


K 


Sin 


k  (sinvi+  sinvc)  h 


-.2 


k  (sinVj  +  sinVg)  h 


C 

where  a  (Vj,  <t>' )  is  defined  as  the  bistatic  scattering  width  for  an  infinitely  long 
cylinder  at  oblique  incidence,  which  we  can  use  for  a  finite  cylinder.  It  is  defined 
for  the  TM  and  TE  polarizations  and  is  given  by: 


c  TM 

( <t>’ ) 


cTl 


kcos  V 


n_TM  ind' 


2-  (-1)  Cn  e 


kcos  V 


n_TE  ino' 


I(-1)  C^e 


with 


TM  Jn(kaoCOSV)  ^TM 

"*“7T)  = 

Hn  (kaocosv) 


and 


TE  Jo(kaoCosv)  ^TE 

(iT  “ 

Hn  (kaoCosv) 

where  is  the  Bessel  function  of  the  first  kind  and  Hp  ^  is  the  Hankel  function  of 
the  first  kind.  The  prime  on  the  J  and  H  denotes  first  derivative.  In  the  LDR 
software,  the  infinite  series  is  truncated  to  only  go  from  n=-10  to  10. 


32 


Case  2.  Cylinder  radius  and  length  much  greater  than  wavelength 
(ka>10,2h>X) 


With  this  case,  physical  optics  gives  an  expression  for  the  bistatic  RCS 
that  is  valid  for  directions  at  and  near  the  specular  direction  (to  within  about  45° 
on  either  side  of  the  specular  direction).  In  this  specular  region,  the  strong 
return  from  the  cylinder  sides  makes  the  contributions  from  the  ends  and  edges 
negligible.  For  the  perfectly  conducting  cylinder,  the  bistatic  RCS  for  either  TE 
or  TM  polarization  is  given  by: 


2  2 

4  k  ah  Vs 
=  - ^cos(i)72 

It  COSVj 


sin 


k  {sin\yj+  sinvj/  )  h 


k  (sin\yj+  siny.)  h 


where  2h  is  the  cylinder  length  and  a  is  the  radius. 

In  the  specular  direction,  Vg  =  ~Vi .  and  the  cross  section  and  scattered 
power  are  maximum,  or 

2 

a  (y j  ,-y j ,(!)')  =  4  k  ah  cos  y,  cos  (J)72 


In  the  backscattering  direction.  Vs  -  Vi  and  ())’  =  0. 
incidence  angles  near  broadside  becomes 

-.2 


The  RCS  for 


a  (yj  .y; ,  0 )  =  4  k  ah  cos  y , 


sin  (2k  h  sinyj) 


2k  h  siny. 


and  at  broadside  (Vi=0),  the  backscattering  RCS  is 
c  (0,0,0)  =  4kah^ 


We  now  present  results  of  the  bistatic  RCS  of  the  cylinder  using  the 
equations  in  this  Section.  The  cylinder  used  to  model  the  fuselage  of  the  cruise 
missile  is  19  ft  long,  and  has  a  radius  of  10  in.  We  will  be  using  the  geometry 
shown  below: 
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A 


Figure  3.1.2.2.3-2.  Bistatic  scattering  geometry  for  cylinder 


The  azimuth  angle  is  ^  and  the  elevation  angle  is  0.  Figures  3.1.2.2.3-3 
(a),  (b),  and  (c)  show  the  bistatic  PCS  of  the  cylinder  as  a  function  of  the 
scattering  elevation  angle  9.  Figure  3. 1.2. 2. 3-3  (a)  is  the  case  where  (})|  =  45° 
and  Gj  =  12°.  As  stated  above,  we  would  expect  scattering  in  the  specular 


AZIMUTH  I  •  45  DEO. 


Figure  3.1.2.2.3-3  (a).  Bistatic  PCS  of  cylinder  (§>  12°  incident  elevation 

X-band 


34 


direction  to  be  dominant  since  the  cylinder  used  in  the  missile  has  length  much 
greater  than  radius.  In  this  geometry,  the  specular  direction  is  where  0  =  Bj.  As 
shown  in  the  plot,  the  peak  RCS  occurs  for  0  around  12°  just  as  expected  in  the 
specular  direction.  The  highest  sidelobe  is  about  13  dB  down  from  the  peak. 
The  magnitude  of  the  peak  is  about  26  dBsm,  which  shows  that  in  the  specular 
direction  the  cylinder  scattering  is  considerable. 

Next,  Figures  3.1.2.2.3-3  (b)  &  (c)  both  show  further  that  specular 
reflection  is  the  dominant  scattering  contribution  from  the  cylinder.  Figure 
3.1.2.2.3-3  (b)  is  the  case  where  (j))  =  45°  and  Bj  =  45°  while  in  Figure  3.1 .2.2. 3-3 
(c)  is  the  case  where  =  0°  and  0;  =  75°.  We  see  that  in  both  cases  the  peak 
goes  as  high  as  32  dBsm  while  it  narrows  at  the  same  time. 


AZIMUTHJ  .  45  DEG. 

Rgure  3.1 .2.2.3-3  (b).  Bistatic  RCS  of  cylinder  @  45°  incident  elevation 

X-band 
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A21MUTHJ  -0  DEG. 


Rgure  3.1.2.2.3-3  (c).  Bistatic  RCS  of  cylinder  (S)  75°  incident  elevation 

X-band 


At  L-band  (1.3  GHz)  specular  reflection  is  also  strong,  as  shown  in 
Figures  3.1.2.2.3-4  (a)  &  (b).  Figure  3.1.2.2.3-4  (a)  is  the  case  where  bj  =  0°,  e,  = 
15°  and  TM  mode,  and  Figure  3.1.2.2.3-4  (b)  is  the  case  where  (fij  =  45°,  e,  =  60°, 
and  TE  mode.  The  peaks  in  these  two  Figures  are  14  dBsm  and  18  dBsm, 
which  are  smaller  than  in  the  X-band  case.  Also,  because  the  frequency  is 
lower,  the  mainlobes  are  much  bigger  in  width,  especially  in  Figure  3. 1.2. 2.3-4 
(a). 


Now  we  look  at  the  situation  when  the  scattering  vector  is  on  the  same 
side  of  the  z-axis  as  the  incident  vector  in  Figure  3. 1.2. 2. 3-2  above.  This  would 
include  the  monostatic  backscattering  case.  We  have  two  cases  to  illustrate  this 
situation.  Figure  3. 1.2. 2. 3-5  (a)  is  for  the  same  incident  vector  as  in  Figure 
3.1 .2.2.3-3  (b)  above  at  X-band.  As  can  be  seen,  the  scattering  is  very  weak  in 
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RCS  (DBSM) 


this  region  by  comparison.  Even  in  the  backscattering  direction  (0  =  45=’)  the 
RCS  is  only  -25  dBsm.  Figure  3. 1.2. 2. 3-5  (b)  is  for  the  same  geometry  at  L- 
band  as  in  Figure  3.1.2.2.3-4  (a)  above.  Once  again,  the  scattering  is  weak'in 
this  region  by  comparison,  about  10  to  20  dB  weaker. 
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AZIMUTH  I  *45  DEG. 


Figure  3. 1.2.2. 3-5  (a).  Bistatic  RCS  of  cylinder  in  backscatter  region 

@  45^  incident  elevation  -  X-band 


The  RCS  formulas  for  the  cylinder  were  not  derived  using  geometric 
optics,  and  therefore  a  specular  point  is  not  defined.  The  cylinder  forms  the 
central  portion  of  the  missile  body,  and  we  assume  that  this  portion  is  never 
shadowed  by  any  other  component  of  the  missile. 

The  two  ends  of  the  missile  are  modeled  as  hemispheres,  which  are  half¬ 
ellipsoids  that  are  easily  modeled  as  before.  We  can  define  specular  points  for 
these  two  components,  and  they  can  be  shadowed  by  the  cylinder  at  certain 
aspect  angles.  For  example,  at  head-on  incidence,  the  rear  hemisphere  is 
shadowed. 
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AZIMUTH  1.0  DEG. 


Rgure  3.1 .2.2.3-5  (b).  Bistatic  RCS  of  cylinder  in  backscatter  region 

@  15°  incident  elevation  -  L-band/TM  mode 


A  comparison  of  the  results  of  the  LDR  model  and  the  Lincoln  Lab  model 
of  the  capped  cylinder  is  shown  in  Figure  3.1 .2.2.3-6.  As  can  be  seen,  the  two 
models  agree  quite  well  with  each  other.  In  both  plots  the  peak  RCS  at 
broadside  is  about  24  dBsm.  The  RCS  at  small  aspect  angles  average  about  - 
10  dBsm  in  both  plots.  The  very  deep  nulls  in  the  LDR  model  is  due  to  the  TM 
mode,  the  TE  mode  does  not  have  them.  It  is  possible  that  the  Lincoln  Lab  plot 
is  due  to  a  combination  of  both  modes,  and  therefore,  it  has  the  shallow  nulls  at 
high  aspect  angles  plus  the  peak  at  broadside.  In  any  event,  the  similarity  of 
results  is  very  encouraging  to  the  further  development  of  the  LDR  missile  model. 
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re  3. 1.2. 2.3-6.  Lincoln  Lab  capped  cylinder  model  vs.  LDR  model 


3.1. 2.2.4  Triangular  plates 

Using  physical  optics  and  Equation  3.1.2.2.2*1  above,  we  can  derive  the 
RCS  of  the  triangular  plate  for  any  aspect  angle.  The  coordinate  system  used  in 
the  derivations  are  shown  in  Figure  3.1.2.2.4-1.  As  in  the  cylinder,  the  plane  of 
incidence  is  defined  as  the  plane  with  the  incident  wave  propagation  direction 
and  the  x  axis.  The  plane  of  scattering  is  defined  as  the  plane  with  the 
scattering  propagation  direction  and  the  x  axis.  It  is  useful  to  note  that  the 
incident  and  scattering  unit  vectors  are  given  respectively  by: 


k  =(cos0cos(p,  cos9sin(p,sin0) 
k  *  (cos0’  coscp',  cos0'sin(p’,sin0') 


Z  ^s 


Figure  3.1 .2.2.4- 1 .  Bistatic  geometry  for  triangular  plate 

We  need  to  derive  the  expressions  for  the  incident  field  H  q  above.  There 
are  2  polarizations  to  consider:  TM  and  TE. 

Case  1 .  TM  mode 

The  H  field  is  perpendicular  to  the  incidence  plane  while  the  E  field  is  in 
the  incidence  plane.  This  H  field  is  perpendicular  to  both  the  incidence  vector 
and  the  x  axis,  or: 


TM  ^ 

Ho  =  Ho(k  X  X  ) 

»  Ho(0,  sin0  ,  -cos0sin(p  ) 
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=  Ho(0.k3,-I<2) 

where 

X  is  the  unit  vector  along  x, 

and  Hg  is  a  constant  proportional  to  the  H-field  magnitude. 

Case  2.  TE  mode 

The  H  field  is  in  the  incidence  plane  while  the  E  field  is  perpendicular  to 
the  incidence  plane.  This  H  field  is  perpendicular  to  both  the  incidence  vector 

and  Ho  ,  or: 


TE  TM 

Ho  =  Ho(k  X  Ho  ) 

2  2  2  2 

=  Hq  ( -cos  9  sin  (p  -  sin  6,  cos  6  sincp  coscp,  sin0cos0  coscp) 

=  Ho(-k2*k3,kik2,kik3) 


Finally,  using  the  equation  for  the  scattered  field,  Equation  3.1.2.2.2-1, 
we  can  write  from  the  definition  of  cross-section: 


O  =  471 


lim 


H 


TE/TM 


K 


r®  TljTE/TM  -1 

k  x[Ho  xzj 

1  n* 

. 

.  .2 

-  exp 

ik  ( k  -  k  )  •  r 

dxdy 

L.  te/tm 

1  *. 

The  double  integral  is  integrated  over  the  triangular  surface  in  the  x'y’  plane. 
We  wil  define  the  variable  J  by: 


J  = 


k  x[h 


TE/TM  - 
0  XZ 


]|Vk 


TE/TM 


TE/TM 

Ho 
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TE/TM 
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For  the  monostatic  case,  the  cross  section  is  the  same  for  either 
polarization  and  is  given  by: 


where 


•iaZ  (ibZtana)/2sinbX 
e  e 


-  ^s\n  (bW/2) 


X  =  k  (coscp  cos9  tan  a  -  sincp  cos9)  ^ 

Z  =  2  k  cos(p  cos9^ 

W  =  2  k  sincp  cos9  _ 

and  a,  b,  and  a  are  the  dimensions  and  angle  of  the  triangular 
plate  as  in  the  Figure  in  Section  3.1. 2.2.1. 


If  Z  =  0  then  the  cross  section  is  given  by: 


(Eqn.  3.1.2.2.4-1) 


In  the  bistatic  case,  the  cross-section  is  different  for  each  polarization.  It 

is  given  by:  ^ 

4k^j  -iaZ  (ibZtana)/2  sin  (bQ/2)  sin  (bW/2) 

^  e  e  - 0 - - ^ - 

where 

W=k(sin(pcos9  +  sincp’ cos9’ )  (Eqn.  3.1.2.2.4-2) 

Z  =  k  ( coscp  cos9  +  coscp'  cos9’ ),  and 
Q  =  Z  tan  a  -  W 

The  only  term  that  changes  with  polarization  is  the  term  J.  When  Z  =  0,  the 
cross  section  is  given  by  Equation  3. 1.2. 2. 4-1  above,  except  that  the  bistatic 
definition  for  W  in  Equation  3.1.2.2.4-2  is  used. 
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As  stated  before  in  the  section  on  bistatic  scattering,  it  can  be  shown  with 
physical  optics  that  in  the  high  frequency  case,  most  of  the  scattered  field  comes 
from  the  immediate  neighborhood  of  the  specular  point,  as  long  as  k  is 
sufficiently  large®.  Thus,  we  would  expect  that  for  the  flat  plate  the  strongest 
scattering  is  in  the  specular  direction,  just  as  in  a  mirror. 


Figure  3.1. 2.2, 4-2  (a)  &  (b)  shows  the  bistatic  RCS  of  a  plate  for  the  case 
where  the  incident  vector  makes  a  45=*  elevation  angle  (Sj  =  45")  with  the 
surface  of  the  plate  as  shown  in  the  small  figure.  The  elevation  angle  is  defined 
as  the  angle  0  in  the  small  figure.  The  frequency  used  is  X-band  (10  GHz).  We 
would  expect  the  strongest  scattering  when  9  =  45",  and  this  plot  certainly  verify 
this  fact.  Even  though  the  mainicbe  is  quite  large,  the  sidelobes  are  about  30 
dB  down. 


Piste 


Figure  3.1 .2.2.4-2  (a)  &  (b).  Eistatic  RCS  of  triangular  plate 

@  45"  incident  elevation  -  X-band 
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Next  we  have  a  case  where  0j  =  30®.  Figure  3.1.2.2.4-2  (c)  shows  that 
the  strongest  scattering  is  in  the  neighborhood  of  G  =  30®.  The  mainlobe  is 
much  bigger  here,  probably  because  the  incident  vector  grazes  the  surface  of 
the  plate  more  closely.  Still,  if  we  go  further  out.  the  scattering  power  is  at  least 
20  dBsm  down  from  the  peak. 


ELSV  (DEG.) 

A21MUTHJ  .  30  DEG. 

Rgure  3.1 .2.2.4-2  (c).  Bistatic  RCS  of  triangular  plate  @  30®  incident  elevation 

X-band 


Next  we  have  the  case  where  the  incident  vector  is  almost  head-on  to  the 
plate,  making  an  elevation  angle  of  75®  with  the  plate.  Figure  3. 1.2.2. 4-2  (d) 
shows  that  the  strongest  scattering  does  indeed  occur  for  the  specular  reflection 
(9  =  75®).  The  mainlobe  is  even  narrower  than  the  45®  case  of  Figure  3.1 .2.2.4- 
2(b). 
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Rgure  3.1 .2.2.4-2  (d).  Bistatic  PCS  of  triangular  plate  @  75=  incident  elevation 

X-band 

At  L-band  (1.3  GHz),  however,  the  physical  optics  formulas  derived 
above  do  not  give  good  specular  reflection  results  for  the  bistatic  PCS  of  the 
plate.  This  is  shown  in  Figures  3.1.2.2.4-3  (a)  &  (b).  In  Figure  3.1 .2.2.4-3  (a)  the 
incident  vector  makes  an  elevation  angle  of  15°  with  the  surface  of  the  plate,  yet 
the  peak  PCS  is  somewhere  around  40°  elevation  angle.  Even  though  the 
PCS  at  0  =  15°  is  only  2  dB  or  so  down  from  the  peak,  specular  reflection  is  not 
as  definite  in  this  case  as  with  X-band.  Next,  Figure  3.1.2.2.4-3  (b)  suggests 
that  when  0j  is  larger,  then  the  peak  scattering  will  tend  to  occur  more  at  the 
angle  9  =  0j.  However,  the  mainicbe  is  still  much  too  large  compared  to  the  X- 
band  case. 
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AZIMUTH  J  -  30  DEG. 

Rgure  3.1.2.2.4-3  (a).  Bistatic  RCS  of  triangular  plate  @  15^  incident  elevation 

L-band/TM  mode 


0  10  20  30  40  SO  60  70  80  90 

ELEV  (DEG.) 

AZIMUTHJ  -  30  DEG. 

Figure  3.1 .2.2.4-3  (b).  Bistatic  RCS  of  triangular  plate  @  60®  incident  elevation 

L-band/TM  mode 
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The  backseat! e ring  region  for  the  same  triangular  plate  is  shown  below: 

\  I  ^ 

\  N  1 


^  '  I 


We  would  expect  the  scattering  to  be  weak  on  this  side  of  the  plate 
because  this  is  not  in  the  specular  region.  This  is  definitely  the  case,  as  shown 
in  Figures  3. 1.2. 2.4-4  (a)  &  (b).  Figure  3. 1.2. 2.4-4  (a)  shows  the  X-band  case 
where  Qj  =  30°.  Scattering  is  nonexistent  for  small  0;  compare  with  Figure 
3. 1.2.2. 4-2  (c)  above.  In  fact,  at  the  backscattering  angle  (0  =  30°),  the 
scattering  is  beneath  the  minimum  plotted.  The  scattering  climbs  up,  however, 
as  0  gets  closer  to  the  normal  angle.  This  is  reasonable  because  in  this  region 
we  move  toward  the  specular  region  on  the  other  side  of  the  normal  vector  z. 


'o  10  20  30  40  50  60  70  80  90 


ELEV  (DEG.) 

A2IMUTHJ  -  30  DEG. 

Figure  3.1.2.2.4-4  (a).  Bistatic  RCS  of  triangular  plate  in  backscatter  region 

@  30°  incident  elevation  -  X-band 
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For  the  L-band  case,  there  is  more  scattering  in  this  backscattering 
region  than  for  X-band.  Figure  3.1.2.2.4-4  (b)  shows  results  for  the  same 
geometry  but  at  L-band  (1 .3  GHz).  The  RCS  of  the  plate  in  this  region  is 
between  -30  dBsm  and  about  -10  dBsm.  Comparing  this  Figure  with  Figure 
3.1.2.2.4-3  (a),  we  see  that  the  RCS  is  10  to  20  dB  lower  in  this  region. 


Figure  3.1.2.2.4-4  (b).  Bistatic  RCS  of  triangular  plate  in  backscatter  region 

(§>  30®  incident  elevation  -  L-band/TE  mode 
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3.1. 2.2.5  Straight  Edges 

The  edges  of  the  wings,  tailfin,  and  tailplane  can  contribute  significant 
scattering.  We  model  these  elements  of  the  missile  by  assuming  that  their 
scattering  behavior  is  similar  to  that  of  thin  wires.  Thin  wires  can  be  thought  of 
as  cylinders  with  radius  much  smaller  than  the  wavelength.  When  this  is  the 
case  it  makes  essentially  no  difference  how  the  cylinder  ends  are  terminated  as 
far  as  the  scattered  fields  are  concerned.  The  cylinder  ends  may  be  flat, 
rounded,  or  pointed,  and  the  results  are  the  same.^ 


There  are  two  cases  of  interest.  The  first  case  is  when  the  length  of  the 
wire  is  of  the  same  order  as  the  wavelength.  This  is  same  situation  as  in 
antenna  theory,  the  dipole.  There  are  two  simple  expressions  for  the  bistatic 
cross  section  of  the  thin  wire  when  its  length  is  the  half-wave  dipole  and  the  full- 
wave  dipole. 


For  the  half-wave  dipole  the  bistatic  cross  section  is; 


Vs)  =  0-86 


2  2 
COS  YjCOS  Yg 


cos  (71/2  sin  \|/|) 

2 

cos  [k/2  sin  x/gj 

cos  \j/j 

cosvg 

where 

¥i  =  the  angle  between  the  direction  of  incidence  and  the  vector  that  is 
perpendicular  to  the  wire, 

¥s  =  the  angle  between  the  direction  of  scattering  and  the  vector  that  is 
perpendicular  to  the  wire, 

Yi  =  the  polarization  angle  of  the  E-field  measured  from  the  plane  of 
incidence.  Yi  =  0®  is  the  TM  mode  while  Yi  *  90°  is  the  TE  mode, 

Ys  ■  the  same  as  T  except  with  respect  to  the  plane  of  scattering. 

As  can  be  seen  from  the  expression  above,  the  cross  section  is  zero  for  the  TE 
mode  because  Yi  =  90°  and  Ys  =  90°  in  that  mode. 

For  the  full-wave  dipole,  the  bistatic  cross  section  is  given  by: 
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7i.  Jt 


(Vi.  Vs) 


2  2 

=  0.93  X  cos  Y| 


2 

cos  y. 


sin  (ti  sin  \|/|) 

2 

sin  (jt  sin  \yj 

2 

2 

COS  \\f■^ 

cos  \|/g 

When  the  radar  wavelength  in  the  LDR  simulation  is  about  between  0.5  ft 
and  1  ft,  some  of  the  wire  edges  will  have  length  on  the  order  of  a  wavelength. 
In  these  cases,  we  will  make  the  approximation  that  the  wire  is  either  a  half¬ 
wave  or  a  full-wave  dipole. 


The  other  case  of  interest  is  when  the  length  of  the  wire  is  more  than 
about  X.  Then  the  general  bistatic  scattering  cross  section  is  given  by: 


2  cos  Va  2  2 

o(Vi.Vs)  =  4jth  - ^cos  YjCOS  Ys 

cos  Yj 


1 


2  2 
In  (0.8905  k  a  cos  \j/i)  + 7c /4 


sin 

k|sin\|rj+sin\yg|h 

.  ''I 

sinVj+sinYjjh 

where 

a  =  radius  of  wire,  which  we  will  take  to  be  Ay85.^o 

Results  for  the  cross  sections  of  the  wires  are  shown  below.  Figures 
3.1.2.2.5-1  (a)  &  (b)  show  the  monostatic  RCS  of  different  wires  for  the  case  of 
0°  and  45"  azimuth.  The  wires  are  all  assumed  to  lie  along  the  x-axis.  The 
elevation  angle  is  defined  as  the  angle  0  in  the  small  figure. 
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Figure  3.1.2.2.5-1  (a).  Monostatic  RCS  of  wires  @  0®  azimuth  -  L  band/TM 

Except  for  the  full-wave  dipole  (L=  X),  all  the  wires  have  the  peak  RCS  at 
broadside  (0  =  90®).  When  the  elevation  angle  is  below  10®,  the  RCS  is  very 
small  or  zero,  except  for  the  full-wave  dipole.  Curiously  enough,  at  broadside, 
the  RCS  of  the  full-wave  dipole  is  zero.  Another  thing  to  note  is  that  the  longer 
wires  have  peak  RCS  between  -5  to  0  dBsm. 

The  case  of  45®  azimuth  is  shown  in  Rgure  3.1.2.2.5-1  (b).  It  is  similar  to 
the  0®  azimuth  case,  except  that  the  wires  which  are  not  the  full-wave  dipole 
have  higher  RCS  at  low  elevation  angles  than  before.  This  is  especially  true  of 
the  half-wave  dipole. 


wire 
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Figure  3.1 .2.2.5-1  (b).  Monostatic  RCS  of  wires  @  45®  azimuth  -  L  band/TM 


We  also  have  bistatic  RCS  results  for  the  wire.  The  geometry  shown  in 
Figure  3.1.2.2.5-2  (a).  The  incident  vector  is  at  0®  azimuth  but  has  a  45® 
elevation  angle  with  respect  to  the  wire.  The  scattering  vector  is  varied  as  a 
function  of  the  elevation  angle  8. 


Rgure  3.1.2.2.5-2  (a) 

The  bistatic  RCS  of  the  same  4  wires  are  shown  in  Figure  3.1 .2.2.5-2  (b).  It  can 
be  seen  that  with  this  geometry  the  full-wave  dipole  gives  the  highest  RCS 
value  through  most  of  the  aspects.  The  other  wires  have  similar  shapes  as  the 
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3.1. 2.2.6  Computed  RCS  of  Missile 


In  this  section  we  present  the  results  for  the  computed  RCS  of  the  cruise 
missile.  The  RCS  is  a  combination  of  all  the  individual  equations  we  presented 
above  for  the  different  types  of  shapes. 

Figure  3.1.2.2.6-1  (a)  is  a  monostatio  RCS  plot  of  the  missile  varying  with 
azimuth  angle.  The  radar  is  at  0®  elevation,  i.e.  in  the  same  xy-plane  as  the 
missile.  0®  azimuth  is  when  the  radar  is  nose-on  with  the  missile  while  180®  is 
the  radar  looking  at  the  tail  of  the  missile.  The  frequency  is  L-band  (1 .3  GHz). 


Figure  3.1.2.2.6-1  (a).  Monostatic  RCS  of  missile  @  0®  elevation  -  L  band 

As  the  plot  shows,  the  average  RCS  of  the  missile  for  the  TE  mode  over 
most  aspects  is  about  -10  dEsm,  except  fcr  the  region  near  broadside  (90®).  In 
that  region,  as  is  expected,  the  missile  has  a  peak  RCS  of  about  10  dEsm.  The 
average  RCS  is  about  5  to  10  dEsm  less  than  the  average  RCS  of  the  aircraft 
for  the  same  radar  geometry.  Recall  that  the  TE  mode  is  the  mode  where  the 
electric  vector  is  perpendicular  to  the  plane  of  the  missile  cylinder  axis  and  the 
direction  of  incidence.  Thus,  the  electric  vector  is  not  aligned  with  the  cylinder 
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axis.  In  contrast,  the  TM  mode  has  the  electric  vector  aligned  in  the  same  plane 
with  the  cylinder  axis.  This  is  a  more  favorable  alignment  and  as  expected  the 
average  and  peak  RCS  are  larger  than  in  the  TE  mode.  The  average  RCS  is 
about  0  dBsm  while  the  peak  at  broadside  jumps  up  to  25  dBsm.  Thus  the 
polarization  makes  a  big  difference  in  the  missile  RCS  at  L-band. 

In  Rgure  3.1.2.2.6-1  (b)  we  show  the  RCS  results  for  the  case  of  a  radar 
pointing  down  at  the  missile  with  an  elevation  angle  of  -45=.  Once  again  the  TM 
mode  gives  far  larger  RCS  values  than  the  TE  mode.  At  broadside  the  TE  mode 
does  not  even  have  a  peak  while  the  peak  for  TM  there  is  about  24  dBsm. 


Figure  3.1.2.2.6-1  (b)  Monostatic  RCS  of  missile  @  -45’  elevation  -  L  band 


At  X-band  (10  GHz),  however,  the  difference  between  the  two 
polarization  modes  are  not  as  pronounced.  This  is  shown  in  Figures  3.1.2.2.6-2 
(a)  &  (b).  The  mainlobe  region  around  broadside  has  about  the  same  width  and 
amplitude  for  both  polarizations.  This  can  be  explained  by  the  fact  that  the 
cylinder  contribution  dominates  in  this  region,  and  the  RCS  equation  for  the  cyl- 
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inder  when  the  frequency  is  high  (ka  >  10)  is  the  same  for  either  polarization.  In 
the  other  aspect  regions,  the  TM  mode  stiil  tends  to  give  higher  RCS  values 
than  the  TE  mode. 

In  Rgures  3.1.2.2.6-3  (a)  &  (b)  we  present  the  X-band  case  where  the 
radar  looks  down  with  an  elevation  angle  of  *45'^.  Once  again,  the  two  modes 
give  essentially  the  same  RCS  results  with  little  differences. 


iCOCT  91  S.I/PU 
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Figure  3.1.2.2.6-3  (a).  Monostatic  RCS  of  missile  @  -45°  elevation 

X-band  -  TE 
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Figure  3.1 .2.2.6-3  (b).  Monostatic  PCS  of  missile  @  *45®  elevation 

X-band  -  TM 


And  finally  in  Figure  3.1.2.2.6-4  we  have  the  RCS  of  the  missile  as  a 
function  of  the  elevation  angle  of  the  radar  with  respect  to  missile.  -90® 
elevation  in  the  plot  is  when  the  radar  is  on  the  ground  looking  straight  up  while 
90®  is  when  the  radar  is  looking  straight  down  at  the  missile.  The  radar  has  0® 
azimuth  angle,  i.e.  the  nose  aspect.  The  RCS  for  both  modes  is  fairly  symmetnc 
with  respect  to  0®  elevation.  As  expected  at  L-band,  the  TM  mode  gives  higher 
RCS  values,  especially  around  0®  elevation. 

These  results  indicate  that  polarization  plays  an  important  role  at  the 
lower  range  of  radar  frequencies.  The  TM  mode,  where  the  electric  vector  is  in 
the  same  plane  as  the  missile  main  axis,  results  in  greater  scattering  than  the 
TE  mode.  Missile  RCS  averages  around  -10  dBsm  at  most  aspects  except  near 
broadside  where  the  RCS  can  climb  up  to  as  high  as  30  dBsm.  At  higher  radar 
frequencies,  the  missile  RCS  tends  to  be  the  same  regardless  of  polarization. 
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Figure  3. 1.2.2. 6-4.  Monostatic  RCS  of  missiie  for  radar  @  0®  azimuth 

L  band 


In  Figure  3.1 .2.2.6-5  we  compare  the  missile  RCS  obtained  by  our  model 
with  that  of  the  Uncoln  Lab  model.  The  first  plot  is  the  same  RCS  curve  for  the 
TM  mode  that  is  show  in  Figure  3.1.2.2.6-1  (a)  above.  The  second  plot  is  the 
Lincoln  Lab  results  for  the  missile  RCS.  It  is  not  clear  if  the  Uncoln  Lab  results 
are  measured  RCS  or  computed  RCS.  In  any  case,  the  plots  are  quite  similar 
for  most  aspects,  except  that  the  LDR  results  show  a  stronger  peak  at  90®  than 
the  Lincoln  Lab  results,  about  7dBsm  more.  This  could  be  because  of  the 
different  types  of  polarization  used.  Nevertheless,  the  similarity  of  results 
reaffirms  the  validity  of  the  LDR  cruise  missiie  model  to  a  large  extent. 
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3.1.3  Target  Interaction  Model 

After  obtaining  the  RCS  models  for  the  aircraft  and  missile,  monostatic 
and  bistatic,  we  now  develop  the  interaction  model  between  the  two  targets. 
This  section  is  at  the  heart  of  the  phenomenology  that  is  explored  by  the  LDR 
program. 

We  begin  with  the  general  bistatic  radar  geometry  shown  in  Rgure  3.1.3- 
1.  The  aircraft  and  missile  are  at  respective  distances  Rj;^  and  R■Ij^^  from  the 
transmitter,  and  distances  Rpi^  and  Rrj^  from  the  receiver.  R^y^^  is  the  distance 
between  them.  The  radar  energy  follows  a  bistatic  path  as  it  bounces  off  the 
aircraft  at  a  bistatic  angle  p^,  travels  to  the  missile  and  bounces  back  to  the 
receiver  at  a  bistatic  angle  p2.  Restricting  ourselves  for  now  to  just  the  area 

around  the  aircraft,  we  know  that  the  power  density  of  the  incident  wave  as  seen 
by  the  aircraft  is  given  by: 

p  Q 

Power  density  = — i- 
4;r  R?a 

I 

where 

Pt  =  transmitter  power,  and 
Gt  =  transmitter  antenna  gain. 

Transmitter 

/f- 


Figure  3. 1.3-1.  Basic  interaction  path  (path  1) 


62 


Now  suppose  that  the  aircraft  is  modeled  as  just  one  component  and  has 

bi 

a  bistatic  RCS  at  angle  given  by  Then  the  radar  power  scattered  by  the 
aircraft  toward  the  missile  is  given  by: 

Pa  “ - 


Ak  R 


TA 


In  a  sense  the  aircraft  is  acting  as  an  antenna  radiating  energy  to  the  missile, 
and  similarly,  the  missile  is  an  antenna  capturing  energy  from  the  aircraft.  This 
being  the  case,  we  can  introduce  two  quantities:  and  A^.  A^  is  the  effective 

aperture  area  of  the  aircraft  acting  as  an  antenna  while  Am  is  the  same  for  the 
missile.  Not  all  the  energy  scattered  by  the  aircraft  toward  the  missile  will  be 
captured  by  the  latter.  If  the  missile  aperture  is  too  small  compared  to  the 
aircraft  aperture,  then  the  missile  can  only  capture  a  fraction  of  the  energy.  This 
fraction  of  energy  is  given  by 

Am 
Aa’ 

1  . 


if  Am  < 
if  Am  ^  A^ 


This  is  one  of  the  key  elements  in  the  LDR  target  interaction  model.  We  also 
implicitly  assume  that  Ram  «  Rrm- 


The  power  captured  by  the  missile  is  then  given  by; 


M 


P.G 


An  Rja 


bi 

min 


Am 

Aa 


.1 


Next,  this  power  is  reflected  back  to  the  radar  receiver.  The  distance  it  has  to 
travel  in  going  back  is  given  by  R^m  +  Prm-  Also,  assume  that  the  effective 
aperture  area  of  the  receiving  antenna  is  given  by  A,.  Then  the  total  power 
received  by  the  receiver  is  given  by: 

Ar 

(  Ram  +  Prm) 


Pr 


PtGt 

r,2 

Air  R.r> 


bi 

min 


^  1 
Aa 
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Now,  Af  can  be  expressed  in  terms  of  the  antenna  gain  though  the 
relationship: 


Ar  = 


G,X 
4  K 


Also,  the  effective  antenna  aperture  areas  of  the  aircraft  and  the  missile  can  be 
related  to  the  bistatic  RCS  of  each  through  the  expressions: 


bi 

= 


4jt  Aa 


,  and  a 


bi 

M 


4jcAm 


bi 

where  is  the  bistatic  cross  section  of  the  missile  at  bistatic  angle  p2> 


Therefore,  we  can  finally  express  the  interaction  power  received  by  the 
receiver  as: 


PlAfcfi  =  Pf  = 


PtG  X 


2  2  /  5i  BT 

V  <j»  Ow 


3 

47c  Rta(Pam+Prm) 


bi  bi 
if  aA>aM 


2  2  bi 
P.G^X  Ga 


471  Rta(PaM+Prm) 


bi  bi 
if  Ga^Gm 


where  we  have  also  assumed  that  Gt  =  Gf  =  G. 


This  equation  for  Pjamr  is  the  basic  interaction  equation  between  an 
aircraft  and  a  missile.  It  assumes  that  the  interaction  path  is  from  transmitter  to 
aircraft  to  missile  and  back  to  receiver.  There  is  another  signal  path  possible.  It 
is  from  the  transmitter  to  missile  to  aircraft  and  back  to  receiver.  This  path 
scenario  is  shown  in  Figure  3. 1.3-2.  The  expressions  for  this  path  is  almost 
identical: 
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TMAR  = 


P 


2  2  /  b?  bi 

PtG  X  Voa^m 
4jt  R-im(RaM+ ^Fw) 


bi  bi 
if 


Pt 


2  2  bi 

Cm 


3  2 

47c  Rtm(Ram+Pra) 


2 


bi  bi 
if  Cm^C;^ 


t 


where  now 

bi 

-  bistatic  RCS  of  missile  at  angle  P3, 

bi 

Ca  =  bistatic  RCS  of  aircraft  at  angle  P4, 

Transmitter 

/K 

\ 

\ 

\ 

\ 


Figure  3. 1.3-2.  Alternative  interaction  path  (path  2) 


The  expressions  for  Ptamr  Ptmar  starting  point  that  enables 

us  to  derive  the  full  interaction  equations  in  the  general  case  when  the  targets 
consist  of  many  components.  These  two  paths  are  called  the  primary 
interaction  paths.  There  are  two  more  paths  possible.  They  are  illustrated  in 
Figures  3. 1.3-3  and  3. 1.3-4.  In  path  3,  the  energy  is  scattered  from  aircraft  to 
missile  but  back  to  aircraft  then  to  receiver.  Thus,  unlike  paths  1  and  2,  there  is 
a  monostatic  RCS  of  the  missile  to  be  computed  here.  In  path  4,  the  energy  is 
scattered  from  missile  to  aircraft  and  back  to  missile  then  to  receiver.  So  there 
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is  a  monostatic  RCS  of  the  aircraft  to  be  computed  here.  These  two  alternate 
paths  are  called  the  secondary  interaction  paths. 


Transmitter 


Transmitter 


's 


V  R 


TA 


R 


RA 


Figure  3. 1.3-3.  Interaction  path  3 


Figure  3.1. 3-4.  Interaction  path  4 


Receiver 


Receiver 


The  interaction  equations  for  paths  3  and  4  are  derived  in  the  same  way 
as  before.  We  make  use  of  the  same  concept  of  effective  antenna  aperture  area 
of  each  target.  For  path  3,  the  interaction  power  received  by  the  radar  is  given 
by: 
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Pr 


bi, 

PtG 

4jc  Rja  (  2Ram+  Pra) 


mm 


mono 


r-,1 


mm 


bis 


mono 


,1 


GX 

4k 


where 


bii 

c^A  s  bistatic  RCS  of  aircraft  at  angle 

Am  s  effective  aperture  area  of  missile  acting  as  antenna, 

bi 

Aa'  =  effective  aperture  area  of  aircraft  acting  as  antenna  at  angle  pi, 

bij 

Aa  =  effective  aperture  area  of  aircraft  acting  as  antenna  at  angle  P4. 


Using  the  relationship  between  effective  aperture  area  and  RCS  that  was 
presented  before,  we  can  simplify  the  expression  for  path  3  interaction  to 
become: 


2,  2  /  bi,  bij 


TAMAR  = 


PtG^X  Voa  Oa 
Rta(2  Ram+ Pra) 


bi^  mono  mono  bi2 
ifOA>OM  ,Om  XTa 


2  2  /  bit  mono 

PtG  X  Voa 
4k  Rja  (2  Pam  +  Pra) 

2  2  bi,  /  bis 

PtG  X  Oa  Voa 


So  2  / 

4k  Rta(2  Pam+ Pra)  »  <^m 


bh  mono  mono  bi} 
ifOA>aM  .Om  ^Oa 


bii  mono  mono  bij 

ifOA^OM  •<^M  >^A 


2  2  bi, 

PtG  ^  Oa 


Pta  (2  Pam  +  Pra) 


bi,  mono  mono  bis 
ifOA^OM  .Cm  ^Ca 


where 

b'l 

Ca  =  bistatic  RCS  of  aircraft  at  angle  P4, 
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mono 


monostatic  RCS  of  missile. 


And  for  path  4  the  interaction  equations  are: 


TMAMR  = 


2  2  /  bi,  bij 

Pt  G  X  V 
Rtm  (2  Ram  +  I^rm) 


bi^  mono  mono  bi2 


if  Ou  >  a 


M  > 


2  2  j  BTJ 
P,G^  Vc^a 


mono 

A 


Rtm  (2  Ram  +  Rrm) 


bi^  mono  mono  bi^ 
if  >  Ga  .  Ga  ^  Gm 


PjG^X  o 


bi, 

M  VGm 


Rtm  (2  Ram  +  Rrm)  '^^a 


mono 


bi,  mono  mono  bi2 
ifGj^^OA  >®A  X^M 


2  2  bi, 

P.G^X 


2 

Rtm  (2  Ram  +  Rrm)^ 


bi,  mono  mono  bij 

ifOM^<^A  .Oa  ^^^M 


where 

bi, 

^M  -  bistatic  RCS  of  missile  at  angle  |33, 

bij 

^M  =  bistatic  RCS  of  missile  at  angle  32, 

mono 

=  monostatic  RCS  of  aircraft. 


The  interaction  equations  for  the  4  paths  presented  so  far  are  for  the 
ideal  case  of  targets  with  just  1  component.  Also,  the  equations  can  be  used  for 
rough  estimations  of  the  interaction  when  both  targets  are  given  by  a  single 
RCS  value,  i.e.  after  the  component  cross  sections  have  been  summed  up 
appropriately.  To  compute  the  interaction  signals  more  correctly,  however,  it  is 
necessary  to  compute  the  interaction  signal  with  respect  to  each  target 
component.  This  is  because  each  aircraft  component  interacts  individually  with 
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each  missile  component.  Each  of  these  individual  interaction  signals  has  a 
phase  that  is  proportional  to  its  interaction  path  length.  In  the  end,  the  signals 
can  be  summed  algebraically,  just  as  in  computing  the  RCS  of  the  entire  target. 
Each  interaction  path  will  have  a  signal  of  its  own.  We  now  derive  the  full 
interaction  equations  for  each  path. 

3.1. 3.1  Path  1  (TAMm 

The  interaction  signal  power  between  the  i-th  aircraft  component  and  the 
j-th  missile  component  as  seen  by  the  radar  is  given  by: 

(1) 

Pij  =  ‘<1<^i(A)^i(M)  ^i(A)  >  , 

®i(A)  ^i(A)  ^  ^j(M) 

where  =  bistatic  RCS  of  i-th  aircraft  component  at  angle 

=  bistatic  RCS  of  j-th  missile  component  at  angle  p2.  and 

2  ^ 

PtG 

'  ““"3l - 1 

RTA(PAM+PnM) 

When  the  Interaction  signals  from  different  components  are  summed  with 
appropriate  phases,  it  is  necessary  to  work  with  voltage  rather  than  power.  The 
interaction  signal  voltage  is  given  by: 

(1)  1/2 

Vij  =  if  . 

.1/2  ’^2  a  re 

•^1  "  ^i(A)  ^ 

Now  we  sum  over  all  aircraft  and  all  missile  components  to  find  the  interaction 
signal  for  path  1 : 


Vi 


2-^i(A)®  * 


'j(M) 


t€  S  ^ 


where  S ,  =  {i  s.t.  }. 
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S  2  =  {i  s.t.  ^  }, 

<|>jj  is  the  relative  phase  angle  of  the  ij-path.  This  path  is  from  the 
transmitter  to  the  i-th  aircraft  component  to  the  j-th  missile 
component  and  back  to  the  receiver. 


3.1. 3.2  Path  2  (TMARI 

The  interaction  signal  for  this  path  is  similarly  derived.  The  interaction 
signal  voltage  between  the  i-th  aircraft  component  and  the  j-th  missile 
component  as  seen  by  the  radar  is  given  by: 

(2)  1/2  1^^ 

Vjj  =  k2  <yi(A)<^j(M)  ^j(M)  >  ^i(A)  , 

,  1/2  '•^2 

k2  Oj^M)  if  ^j(M)^<^i(A) 

I 

where  =  bistatic  RCS  of  i-th  aircraft  component  at  angle  P4, 

cTj(|^  =  bistatic  RCS  of  j-th  missile  component  at  angle  Pa,  and 
2  2 

u  _  PtG  X 

K2  =  — j— - 

Rtki(Ram  + Rra)^. 

Summing  over  all  aircraft  and  all  missile  components  we  find  the  interaction 
signal  for  path  2; 


V2  = 


I 

Allj 


1/4 

^i(M) 


icS, 


1/4 

^i(A) 


1/2 


e  +  G:, 


j(M) 


i^i| 


where  S,  -  {i  s.t. 

S2  =  {i  s.t.  ^<^i(A)}- 


3.1. 3.3  Path  3  (TAMARI 

The  interaction  signal  voltage  between  the  i-th  aircraft  component  and 
the  j-th  missile  component  as  seen  by  the  radar  is  given  by: 
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.  1/2 

1/4 

1/4 

ks 

^i(A), 

®i(A)2 

if 

<^i(M)><^i(A),  ^ 

1/2 

1/4 

1/4 

if 

ka 

^i(A), 

^j{M) 

^i(A),>®i(M)’ 

®j(M)^<^i(A), 

1/2 

1/2 

1/4 

-1/4 

ka 

^i(A), 

^i(A)2< 

^j(M) 

"  ^i(A),^^KM)»®i(M)><^i(A)j^ 

.  1/2 

1/2 

ka 

^i(A), 

"  1^i(A),^<^j(M)»<^j{M)^<^i(A), 

where  ^^ka),  =  bistatic  RCS  of  i-th  aircraft  component  at  angle  pi, 

®i(A),  =  bistatic  RCS  of  i-th  aircraft  component  at  angle  |34, 

Oj  =  monostatic  RCS  of  j-th  missile  component ,  and 
2  2 

.3-  3.^'^-^ - J 

An  Rta(2Ram+ Rra)^. 

As  in  paths  1  and  2,  the  interaction  signal  for  path  3  can  be  written  as: 

Vc  -  I  Ivf  e'*’ 

Attj  Alii 


3.1. 3.4  Path  4  (TMAMR) 

The  interaction  signal  voltage  between  the  i-th  aircraft  component  and 
the  j-th  missile  component  as  seen  by  the  radar  is  given  by: 


1/2  1/A  1/A 

1/2  1/4  1/4 

*^A  <^KM),^i(A)  ^j(M),>®i(A)»  ®i(A)^®i(M), 

1/2  1/2  1/4  -1/4 

*^A  ^i(M),<^j(M)j<^i(A)  ^j(M),^<^i(A)«  ®i(A)><^KM)2^ 

.1/2  1/2 
•^A 


where  ®km)i  =  bistatic  RCS  of  j-th  missile  component  at  angle  P3, 

s  bistatic  RCS  of  j-th  missile  component  at  angle  P2, 
Oi(A)  =  monostatic  RCS  of  i-th  aircraft  component ,  and 
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2  2 

.  PtG 

K4  =  — - - 

Ptvi(2Ram+ Prm) 

As  in  path  3,  the  interaction  signal  for  path  4  can  be  written  as: 

V4  - 1 

Atlj  Aili 
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Radar 


Figure  3.1 .3.5-1 .  Air-to-air  monostatic  radar  scenario 


The  interaction  signals  between  the  aircraft  and  the  missile  are  computed  as  a 
function  of  time  for  a  monostatic  radar  looking  down  at  the  targets  with  an  initial 
elevation  angle  of  -45*^.  The  aircraft  and  the  missile  fly  along  a  straight 
trajectory  in  parallel  at  a  constant  Mach  1  speed.  The  initial  range  from  the 
radar  to  the  targets  is  50  nm  while  the  interaction  angle  between  the  targets 
stays  constant  at  30®.  The  radar  polarization  used  is  the  TE  mode. 

The  skin  return  signals  for  this  scenario  is  shown  in  Figure  3.1. 3.5-2  (a) 
while  the  interaction  results  are  shown  with  paths  1  and  2  in  Rgure  3.1. 3.5-2  (b) 
and  paths  3  and  4  in  Figure  3.1, 3.5-2  (c).  The  skin  return  signals  are  the  PCS 
signals  received  directly  from  each  target  by  the  radar  without  any  cotarget 
interference. 
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3.1. 3.5  Computed  Tar^pf  Interaction 

This  section  presents  the  results  of  the  interaction  model  developed  thus 
far.  Before  continuing,  a  word  needs  to  be  said  about  the  appropriate  units  for 
the  interaction  signals.  Looking  back  at  the  interaction  equations,  it  will  be  seen 
that  the  interaction  signal  power  received  by  the  radar  is  always  proportional  to 
the  cross  section  a: 


P,  =  ko.  or 

1/2  1/2 

k  Og  Ob  or 
1/2  -1/2 
kOgOb  Oc 

I 

where  k  is  always  in  the  form: 

2  2 

3  2  2 
(4JI)  Rj(  Ry 

Now,  from  the  basic  radar  equation  we  know  that  a  skin  return  signal  received 
by  the  radar  from  a  target  is  given  by: 

2  2 

„  PtG"x  o 

'  “  r~r 

(Ak)  r 

Obviously,  the  received  power  for  the  target  skin  return  is  also  in  the  form  ko, 
and  therefore  it  is  proportional  to  the  cross-section.  When  specifying  the  power 
of  this  signal,  we  do  so  in  terms  of  the  size  of  the  cross-section  (m2).  Thus,  it 
seems  natural  to  specify  the  interaction  signal  strength  in  terms  of  an  equivalent 
cross-section.  The  equivalent  cross  section  has  units  of  m2  and  it  is  given  by 
the  different  forms  shown  above,  depending  on  the  particular  interaction  path: 

1/2  1/2  1/2  1/2  -1/2 
<^0=  or  or  O3  Ob  Oc 

The  interaction  signals  will  be  plotted  in  units  of  dBsm. 

The  first  case  is  for  an  air-to-air  monostatic  radar  at  L-band.  The  scenario 
is  shown  in  Figure  3.1 .3.5-1 . 
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INTERACTION  SIGNAL  (OBSM) 


L-band/TE  mode 


Figure  3.1. 3.5-2  (b).  Interaction  paths  1  and  2  for  mono,  radar  pointing 

@  -45“  elevation.  L-band/TE  mode 


Figure  3.1 .3.5-2  (c).  Interaction  paths  3  and  4  for  mono,  radar  pointing 

@  -45®  elevation.  L-band/TE  mode 

At  about  T=218  sec.,  the  targets  pass  directly  beneath  the  radar  for  a  broadside 
aspect.  Hence,  as  shown  in  Figure  3.1. 3.5-2  (a),  the  RCS  of  the  aircraft  peaks 
up  to  about  30  dBsm  and  the  RCS  of  the  missile  is  about  12  dBsm.  Figures 
3.1. 3.5-2  (b)  &  (c)  show  that  the  Interaction  signals  of  path  1  and  path  3  are  the 
strongest  with  peaks  of  about  20  dBsm,  Path  2  is  slightly  weaker  than  path  1 . 
Figure  3.1 .3.5-2  (c)  shows  that  path  3  dominates  after  about  T=120  sec.  where  it 
is  consistently  20  dBsm.  Path  4  is  about  10  to  20  dB  weaker  In  strength. 

Most  of  the  time  in  this  run,  the  RCS  of  the  missile  is  between  -10  and  10 
dBsm  as  seen  In  Rgure  3.1. 3.5-2  (a).  The  other  figures  here  show  that  2  out  of 
the  4  interaction  paths  have  signals  consistently  stronger  than  10  dBsm. 
Hence,  the  results  suggest  that  the  interaction  signals  provide  us  with  much 
more  information  about  the  low  RCS  target  (the  missile)  than  the  skin  returns  by 
themselves  can.  The  presence  of  these  strong  interaction  paths  alert  us  to  the 
presence  of  another  target  that  may  be  masked  by  the  larger  target.  It  is 
interesting  to  note  that  in  some  instances,  an  interaction  signal  can  even  be 
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stronger  than  the  skin  return  from  the  aircraft.  For  example,  this  happens  with 
paths  1  and  2  at  T=50  sec. 

Next  we  have  the  case  of  a  bistatic  radar  looking  at  essentially  the  same 
geometry.  The  transmitter  and  the  receiver  will  be  offset  150®  in  azimuth 
although  the  initial  elevation  angle  for  both  will  still  be  -45®  (look-down).  The 
aircraft  and  the  missile  will  be  flying  along  a  straight  trajectory  in  parallel  as 
before.  The  initial  range  to  both  the  transmitter  and  the  receiver  is  50  nm.  The 
geometry  is  shown  below.  The  receiver  position  shown  is  not  entirely  accurate, 
and  should  be  slightly  out  of  the  page. 


Transmitter 


Receiver 

y\ 


Rgure  3.1 .3.5-3.  Air-to-air  bistatic  radar  scenario 


The  skin  return  signals  are  shown  in  Figure  3.1. 3.5-4  (a).  For  both 
targets,  the  signals  decrease  slowly  with  time  as  this  run  progresses.  The 
aircraft  has  peaks  of  about  23  dBsm  while  the  missile  has  peaks  at  about  17 
dBsm. 


76 


INTERACTION  SIGNAL 


Q  LOfllSQiNfEMACr  7  JO*r 


TIME  (SEC.) 

Figure  3.1 .3.5*4  (c)  .  Interaction  paths  3  and  4  for  bistatic  radar  pointing 

@  -45^  elevation.  L-band/TE  mode 

A  very  interesting  case  is  the  case  of  a  monostatic  ground  radar  looking 
up  at  the  targets.  This  scenario  is  illustrated  in  Figure  3.1. 3.5-5  below: 

Airplane 


Radar  ^ 

Figure  3. 1.3. 5-5.  Ground-to-ground  monostatic  radar  scenario 
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The  skin  return  signals  are  shown  in  Figure  3.1. 3.5-6  (a).  The  signals  are 
essentially  constant  until  about  T=:245  sec.  This  can  be  explained  as  followed. 
The  initial  range  to  the  targets  is  50  nm  while  the  height  of  the  targets  is  10  kft. 
Hence,  the  initial  elevation  angle  the  radar  points  up  is  1.89®  (radar  at  -1.89® 
elevation  with  respect  to  targets),  so  the  radar  is  looking  up  at  the  underside  of 
the  targets.  Even  at  T=245  sec.,  the  elevation  angle  is  still  only  9°.  As  Figure 
3.1. 2.1. 8-5  shows,  the  RCS  of  the  aircraft  for  small  negative  elevation  angles 
(between  0®  and  about  -10®)  is  mostly  constant  at  5  dBsm.  Similarly,  Figure 
3. 1.2. 2. 6-4  shows  that  with  the  TE  mode  the  RCS  of  the  missile  for  small 
negative  elevation  angles  is  also  fairly  constant. 


TIME  (SEC.) 


Figure  3.1. 3.5-6  (a).  Skin  returns  for  mono,  radar  pointing  @  1.89®  elevation 

L-band/TE  mode 


Then  in  Figures  3.1. 3.5-6  (b)  &  (c)  we  have  the  interaction  signals  for  this 
case.  Only  path  4  is  larger  than  the  skin  returns,  but  not  until  after  T=240  sec. 
The  interaction  signals  are  mostly  about  -10  dBsm,  which  is  the  same  as  the 
missile  skin  return.  This  is  a  case  where  during  most  of  the  run,  the  interaction 
signals  do  not  give  us  more  information  about  co-target  interference.  One 
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reason  why  the  interaction  signals  are  weak  for  this  case  is  that  since  the  radar 
is  on  the  ground,  some  of  the  paths  involve  forward  scattering  toward  the 
shadowed  sides  of  certain  missile  components,  such  as  the  triangular  plates. 
Recall  that  one  of  the  assumptions  of  physical  optics  is  that  there  is  no  surface 
current  on  the  shadowed  side  of  an  object.  If  the  incident  wave  is  on  one  side  of 
the  plate,  the  RCS  model  assumes  that  there  is  no  scattering  toward  the  other 
side  of  the  plate.  Hence,  many  of  the  missile  components  have  zero  RCS  on 
the  interaction  paths.  This  will  reduce  the  interaction  signals. 
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Figure  3.1. 3.5-6  (c).  Interaction  paths  3  and  4  for  mono,  radar  pointing 

@1.89®  elevation.  L-band/TE  mode 


Finally,  we  present  the  case  of  a  bistatic  ground  radar  where  the  receiver 
is  offset  in  azimuth  from  the  transmitter  by  150°.  The  initial  elevation  angle  for 
both  units  is  1 .89°.  The  geometry  is  shown  below: 


Figure  3.1. 3.5-7.  Ground-to-ground  bistatic  radar  scenario 


Receive 
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This  is  another  case  where  as  shown  in  Figures  3. 1.3. 5-8  (a),  (b),  and  (c), 
the  interaction  path  signals  are  mostly  weaker  than  the  skin  returns.  Path  2  is 
the  only  path  that  is  comparable  to  the  aircraft  skin  return.  The  rest  of  the  paths 
are  -10  dBsm  or  less,  which  is  even  less  than  the  missile  skin  return.  Once 
again,  with  a  ground  radar,  many  of  the  missile  components  have  zero  RCS  on 
the  interaction  paths,  and  this  leads  to  weaker  interaction  signals. 


TIME  (SEC.) 


Figure  3.1. 3.5-8  (a).  Skin  returns  for  bistatic  radar  pointing  @  1.89=  elevation 

L-band/TE  mode 
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INTERACTION  SIGNAL  (DBSM) 


TIME  (SEC.) 


Figure  3.1. 3.5-8  (b).  Interaction  paths  1  and  2  for  bistatic  radar  pointing 

@1 .89“  elevation.  L-band/TE  mode 


TIME  (SEC.) 

Figure  3.1 .3.5-8  (c).  Interaction  paths  3  and  4  for  bistatic  radar  pointing 

@1.89“  elevation.  L-band/TE  mode 


Figures  3.1. 3. 5-9  (a),  (b),  and  (c)  show  results  for  the  same  geometry  as 
in  Figure  3.1. 3.5-1  above,  except  with  the  TM  polarization  instead.  For  the  most 
part,  the  TM  case  is  very  similar  to  the  TE  case,  except  that  we  note  the  missile 
skin  return  has  a  peak  at  broadside  of  about  26  dBsm.  This  does  not  occur  for 
the  TE  mode.  We  also  note  that  path  3  is  not  as  smooth  as  in  the  TE  case. 


Figure  3.1 .3.5-9  (a).  Skin  returns  for  mono,  radar  pointing  @  -45°  elevation 

L-band/TM  mode 
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INTERACTION  SIGNAL  (DBSM) 


Figure  3.1. 3.5-9  (b).  Interaction  paths  1  and  2  for  mono,  radar  pointing 


@-45®  elevation.  L-band/TM  mode 


Figure  3.1 .3.5-9  (c).  Interaction  paths  3  and  4  for  mono,  radar  pointing 

@-45®  elevation.  L-band/TM  mode 
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We  also-  now  present  the  results  for  a  case  at  X-band  (10  GHz).  This 
case  is  for  an  air-to-air  bistatic  radar  with  TM  mode.  The  geometry  is  exactly  the 
same  as  in  Figure  3. 1.3. 5-3  above.  We  see  from  Figure  3.1.3.5-10  (a)  that  the 
skin  returns  at  X-band  are  similar  in  magnitude  with  that  of  Figure  3. 1.3. 5-4  (a) 
(L-band  case),  except  that  there  are  much  more  fluctuations  here.  Also,  the 
missile  skin  return  here  has  a  peak  of  32  dBsm.  Paths  1,  2,  and  3  in  Figures 

3.1.3.5-10  (a)  &  (b)  are  similar  in  shape  to  the  L-band  case.  But  path  4  is  much 
weaker  than  the  L-band  case. 


Figure  3.1.3.5-10  (a).  Skin  returns  for  bistatic  radar  pointing  @  -45=  elevation 

X-band/TM  mode 
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Figure  3.1.3.5-10  (b).  Interaction  paths  1  and  2  for  bistatic  radar 

@-45°  elevation.  X-band/TM  mode 
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Figure  3.1 .3.5-1 0  (c).  Interaction  paths  3  and  4  for  bistatic  radar  f 

@-45°  elevation.  X-band/TM  mode 


PATH  1 
PATH  2 


XMTR@  0  OEG.  azimuth 
RCVfl@  1 50  DEG.  AZIMUTH 
XMTR/RCVR  HEIGHT  -2.24e05  FT. 
AIRCRAFT  HEIGHT  -  I  , OE04  FT. 
DlST^MA  -  100  FT. 

INTERACTION  ANGLE-30  DEG. 


3.2 


Given  that  the  large  /small  target  interaction  exists,  a  logical  question  that 
follows  is  wether  any  target  parameters  can  be  extracted  from  these 
signals.That  is,  can  one  form  any  measurement  discriminants.  Figure  3.2-1 
shows  an  interaction  geometry  and  the  appropriate  terms  associated  with  each 
path.  Table  3.2-1  lists  all  of  the  signals,  the  corresponding  path  length  and 
observed  angles  which  form  the  basis  in  forming  any  discriminants.  The 
variables  R^,  R2>  L,  62  are  the  unknown.  With  a  monopuise  radar,  if  the  three 
interaction  signals  are  observed,  then  all  of  the  target  parameters  can  be 
determined. 


Figure  3.2-1  Information  content  embeded  in  the  interaction  signal  allows 
prediction  of  small  target  skin  return. 
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TABLE  3.2-1  Tabulation  of  the  signal  paths  and  observations 


Path 

Signal 

Length 

Observed  Angle 

R  to  T1  to  R 

Yi(t) 

2Ri(t) 

01  (t) 

R  to  T2  to  R 

2  RgO) 

02^) 

R  to  T1  to  T2  to  R 

Plus  vice-versa 

Y3(t) 

Rl(t)  +  L(t)  +  R2(t) 

e,(t)+  02(1) 

2 

R  to  T1  to  T2  to  T1  to  R 

Y4(t) 

2(Ri(t)  +  Li(t)) 

01  (t) 

R  to  T2  to  T1  to  T2  to  R 

Ygit) 

2  (R2(t)  +  L(t) ) 

egfl) 

Notation  :R=Radar  ,T  1  =Target  1 ,  T2=Target2 


With  the  geometry  illustrated  in  the  figure  and  applying  the  law  of 
cosines,  one  obtains 

L  (t)  =  R?(t)  +  R2(t)  -2  R,  (»)  R2(t)  cos  (e,  -  Bj) 

This  equation  can  be  used  to  check  for  consistency  if  one  directly  estimates  the 
target  parameters  from  the  analytically  tractable  discriminants..Consider  a  linear 
least  squares  estimation  that  directly  solves  for  the  target  parameters.  From  the 
table,  one  can  write  the  set  of  vector  equations  in  the  following  two  tables. 
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TABLE  3.2-2.  Direct  parameter  estimation  vector  values  for  the  example 

computation. 


•  Assume  Y2(t)  not  detectable 


Range  and  Doppler  Vectors 


y1  s{2  0  0)B  +  n^ 


y3  si(1  1  1 


y4  =  ( 2  2  0 )  B  +  n^ 


Angle  Vectors 


y1  =  (  1  0)fi  +  n.| 


y5  =  ( 0  2  2 )  B  +  ns 


R  =  (R^  L  R2)^or(Ri  LRg)^ 


and  T  denotes  the  transpose  of  the  row  vector 


Given  obsen/ations 


+  n  :  n-1.2 . N 

Q,  is  the  parameter  vector 
Find  e  that  minimizes  the  least-square  error 


y5  »  (  0  1  )  ft  +  ns 


e  =  (9,  62)''' 


2 

n 

n 

or  Aj-E(e-8)(9-0)*  =  R’ 


where 


Figure  3.2-2.  Least-Squares  Estimate:  Direct  approach  Mathematics 


Figure3.2"3.  An  example  of  angle  estimation  error. 


An  example  of  the  predicted  accuracy  follows.  The  required  mathematics  are 
given  in  Figure  3.2-2  which  includes  the  estimation  error  prediction.  A  sample 
computation  is  given  in  Figure  3.2-3.  The  computation  in  the  figure  assumes 
that  small  target  skin  return  is  not  detected  and  uses  the  weights  from  the  Table 
3.2-1.  The  results  indicate  that  target  parameter  estimation  error  is  proportional 
to  the  discriminant  errors.  Because  enhanced  interaction  paths  provide  small 
target  parameter  estimates  even  with  an  undetectable  skin  return  ,  one  can 
question  the  need  for  co-target  interference  rejection  techniques.  This 
intriguing  result  has  the  following  implications;. 

1.  Standard  resolution  is  sufficient  for  small  target  detection  and 
estimation. 

2.  Dynamic  range  requirement  is  reduced  by  the  interaction  signals. 

3.  Power-Aperture  products  for  small  target  detection  is  reduced  to 
values  as  required  by  the  interaction  signals. 

4.  The  masking  tactic  of  small  target  hiding  near  the  large  target  risks 
detection  via  these  secondary  signals. 
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3.3  Radar  Resoltition  Trade  Study  and  Workload 


In  order  to  property  conduct  a  radar  resolution  trade  study  with  respect  to 
radar  type,  frequency,  waveform,  pulse-repetition  frequency,  pulse  compression 
codes,  dwell  time,  antenna  agility,  etc.,  one  must  quantitatively  define  a 
resolution  "cell”  as  it  applies  to  surveillance  radars.  The  most  tractable 
definition  is  basically  a  four  dimensional  hypercube  p,  defined  as 


where 


(PaZ  ^  two-dimensional  antenna  beamwidth  proportional  to 

the  covered  angular  volume  by  a  beam 

Ar  is  the  radar  range  resolution  (instantaneous  bandwidth) 

Av  is  the  radar  velocity  resolution  (proportional  to  dwell  time 

Rq  a  specified  range. 

it  should  be  noted  that  two  conditions  must  be  satisfied  in  order  for  two 
targets  to  be  resolved  with  conventional  radar  processing.  First,  the  mainlobe 
response  of  the  4D  cell  must  be  narrower  than  the  two  target  separation. 
Second,  the  sidelobe  response  does  not  mask  the  neighboring  targets.  That  is, 
the  sidelobe  response  supports  the  target  dynamic  range.  The  small  target 
masking  levels  establish  the  cancellation  or  suppression  requirements  or 
equivalently  the  sidelobe  relative  to  the  mainlobe  level  required.  To  derive  this 
requirement,  one  must  establish  the  small  and  large  target  RCS  model.  Our 
frequency  resolution  trade  approach  attempts  to  depend  minimally  on  specific 
system  parameters  and  address  the  two  conditions  parametrically  with  respect 
to  wavelength  ( frequency). 

The  first  class  of  systems  examined  are  surveillance  radars  with 
mechanically  scanned  antennas.  The  surveillance  systems  operating  at  VHF, 
L,  S,  C  and  X  are  configured  to  satisfy  the  surveillance  and  track  workloads. 
This  implies  the  same  surveillance  volume,  frame  time,  number  of  target  types, 
track  update  times,  RCS  versus  frequency,  etc.,  which  establish  detection 
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performance  requirements  for  a  given  surveillance  system  when  targets  are 
non-interfering.  A  reasonable  model  for  large  complex  targets  is  that  the 
backscatter  is  independent  of  frequency.  Small  targets,  modelled  as  cone 
spheroids,  have  RCS  that  vary  proportionally  to  a  fractional  power  of 
wavelength  until  a  base  level  limit  (eg.  -30dBsm  @  3GHz)  which  then  remains 
constant  at  that  level  for  all  higher  frequencies.  For  such  a  system  satisfying  the 
surveillance  equation,  the  velocity  resolution  AV  satisfies 

Av=(Xy2)^V8/pTf 

where  ^  is  the  beam  overlap  factor,  Ve  is  the  search  volume,  ^  is  the  two 
dimensional  beamwidth  and  Tf  is  the  frame  time.  Therefore,  the  four 
dimensional  cell  is 

p=  K  Ro2Ar(0  V«/2  Tf)X 

Now  Ar  will  be  evaluated  for  two  design  constraints.  For  a  system 
without  the  processing  constraints  and  equal  RF  complexity,  range  resolution  is 
limited  by  the  system  percent  bandwidth  and  therefore,  the  4D  cell  is 
proportional  to  wavelength  squared.  For  a  system  that  is  processor  limited, 
driven  by  range  bin  processing,  range  resolution  is  independent  of  wavelength 
and  therefore,  the  4D  cell  is  only  proportional  to  wavelength.  Thus,  depending 
on  which  constraint  is  active,  the  4D  resolution  ceil  p  will  be  proportional  to  X  or 
X2. 

TABLE  3.3-1.  System  Configurations  with  HPRF  waveform.  Velocity  resolution 
is  reduced  by  the  range  cell  resolution  factor. 
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TABLE  3.3-2.  System  Configurations  with  LPRF  waveform.  Range  resolution  is 
reduced  by  the  velocity  cell  resolution  factor. 


FREQUENCY 

REQVEL 
COVERAGE  (kHz) 

NO.  OF  FILTERS 

MIN  NUMBER 
PRFS 

CELLFACTOR 

UHF 

2.0 

2 

2 

L 

7.5 

2 

2 

S 

15.0 

4 

3 

3 

C 

30.0 

X 

60.0 

Waveform  ambiguities  reduce  the  resolution  dependence  on  frequency. 
High  PRF  requires  range  resolving  and  depends  on  radar  instanteneous 
bandwidth.  Table3.3-1  list  a  HPRF  waveform  configuration  for  each  operating 
frequency  where  the  range  resolution  is  100ft.  and  velocity  coverage  includes 
target  velocities  up  to  3000  ft/sec.  The  number  of  PRFs  required  to  resolve  to 
200  n.m.  while  satisfying  the  surveillance  constraint,  velocity  resolution  is 
reduced  by  a  value  defined  as  the  range  cell  factor  in  the  table.  Therefore, 
HPRF  velocity  resolution  improves  as  -log  (X). 
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Figure  3.3-1 .  Normalized  Time  on  Target  versus  Operating  Frequency 


Low  PRF  requires  velocity  ambiguity  resolution  which  depends  on 
doppler  filter  resolution  or  equivalently  the  time  on  target.  This  implies  that  a 
radar  workload  must  be  defined  to  determine  dwell  time.  Figure  3.3-1  shows  a 
plot  of  normalized  time  on  target  versus  operating  frequency  for  a  fixed  power 
aperture  product.  The  small  target  RCS  frequency  dependence  is  the  same  as 
used  in  the  surveillance  constraint.  Fixing  range  resolution,  and  using  300Hz 
PRF,  the  cell  factor  is  almost  constant  with  wavelength  as  shown  in  Table  3.3-2. 

Based  on  the  above  discussion,  the  following  can  be  summarized  for  a 
surveillance  system  with  mechanically  scanned  antenna.  The  cell  size  is 
proportional  to  a  value  between 

X  and  X2  LPRF  Waveform,  and 

X  log{M\)  and  X  VogCMX)  HPRF  Waveform. 
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From  a  system  design  point  of  view,  an  electronically  scanned  antenna 
system  decouples  the  surveillance,  track  acquisition  and  track  maintenance 
modes.  With  this  level  of  beam  agility,  resolution  criteria  may  be  more  relevant 
to  track  acquisition  and  track  dwells  rather  than  surveillance.  Consider  an 
Alert/Confirm  Policy  utilized  with  an  electronically  scanned  antenna  radar  in 
which  the  target  resolution  is  only  attempted  during  the  confirm  stage.  This 
allows  time-energy  management  and  adaptive  dwells  as  they  are  required.  For 
example,  it  is  feasible  to  increase  track  dwell  times  for  attaining  resolution  and 
minimize  smearing.  The  predictive  nature  of  tracking  fitters,  as  large  and  small 
RCS  targets  cross  each  other's  path,  allows  one  to  allocate  the  appropriate 
dwell  time  for  resolution  as  well  as  for  detection.  Then  the  40  cell  is  strictly  a 
function  of  the  two  dimensional  beamwidth  and  range  resolution  established 
from  the  point  search  requirement.  This  implies  that  p  is  proportional  to 

and  X*  HPRF  and  LPRF  waveform 

TABLE  3.3-3.  Dynamic  Range  Requirements  versus  Frequency.  Requirements 
are  less  stringent  for  the  low  operating  frequency. 


FREQUENCY 

DYNAMIC  RANGE  REQUIREMENT 

0.25 

22  to  28  dB 

1.0 

34  to  40  dB 

2.0 

40  to  46  dB 

3.0 

40  to  46  dB 

0 


5. 


40  to  46  dB 


TABLE  3.3-4.  Dynamic  Range  Error  Limit.  Optimal  designs  favor  the  highest 

frequency  choice. 


RESOLUTION  DOMAIN 

DYNAMIC  RANGE  ERROR  LIMIT 

Angle 

70  to  90  dB 

Range 

40  to  50  dB 

Doppler 

90dB 

The  dynamic  range  versus  frequency  is  given  in  Table3.3-3.  It  is 
assumed  that  the  large  target  RCS  is  10  dBsm  and  is  independent  of  frequency. 
The  dynamic  range  follows  from  the  dependence  of  the  small  target  RCS  with 
frequency  (described  in  the  above  paragraphs).  The  only  limiting  factor  is  the 
capability  of  designing  the  hardware  so  that  sufficient  dynamic  range  exists  to 
maintain  an  unmasked  small  target.  Table  3.3-4  shows  typical  dynamic  range 
limits  due  to  hardware  tolerances,  it  Is  clear  that  range  resolution  dynamic 
range  tends  to  be  the  most  difficult  as  it  is  limited  40  to  50  dB. 
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3.4  Statistical  Characterization 


The  RCS  model  developed  in  Section  3.1,  the  Phenomenology  Section, 
gives  the  capability  to  predict  target  RCS  as  function  of  target  aspect  angle.  The 
interaction  model  then  quantifies  the  interaction  between  the  targets.  To  apply 
those  models  in  further  statistical  and  detection  studies,  however,  it  is  necessary 
to  account  for  the  cross-section  and  interaction  signal  fluctuations  by 
characterizing  their  probability  distribution  functions. 


A  common  model  of  target  fluctuation  is  the  Swelling  model.  In  the 
Swelling  1  and  2  cases,  the  probability-density  function  for  the  cross  section  is 
given  by 


p(o)  =— exp 
O'av 


a 

a 


\ 

av/ 


o^O 


(Eqn.  3.4-1) 


where  Oav  is  the  average  cross  section  over  ail  target  fluctuations.  In  the 
Swelling  3  and  4  cases,  the  probability-density  function  is  given  by: 


p(o)* 


o 


av 


The  Swelling  cases  1  and  2  have  the  exponential,  or  Rayleigh-power 
distribution.  Analysis  of  measurements  on  actual  aircraft  have  shown  that  the 
cross-section  fluctuations  are  well  fitted  by  this  distribution,  except  at 
broadside^^  We  now  examine  the  cross  section  profile  of  the  aircraft  as 
computed  by  the  LDR  model  in  Section  3. 1.2.1,  and  determine  if  the 
exponential  distribution  fits  the  predicted  profile. 

First  of  all,  a  histogram  of  the  aircraft  cross  section  was  made.  Using 
Figure  3.1. 2.1. 8-1  (a),  the  case  at  L-band  and  0*  radar  elevation  angle,  we 
exclude  the  broadside  peak  from  the  computations.  The  region  with  azimuth 
between  TS*”  and  105*’  was  removed  because  the  peak  may  bias  the  statistical 
characterization.  Then  we  made  a  histogram  of  the  frequency  of  the  cross 
section,  and  this  is  shown  in  Figure  3.4-1.  The  histogram  represents  the 
empirical  probability  distribution  function  of  the  aircraft  cross  section.  The  figure 
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also  shows  the  result  of  curve-fitting  Sweriing  cases  1  and  3  to  the  histogram 
using  the  least-mean  square  criterion.  In  other  words,  we  minimize  the  error 
given  by 


N 

n-1 


2 


Where  is  the  histogram  at  each  RCS  point  x^.  p  is  the  probability  density 
function  being  tested  at  each  RCS  point  x^,  and  N  =2000. 


Figure  3.4*1 .  Histogram  and  Sweriing  models  for  aircraft  RCS 
@0®  elevation.  L-band 


It  can  be  seen  that  the  Sweriing  1  case,  the  exponential  distribution,  fits 
quite  well,  and  that  it  is  a  much  better  fit  than  the  Sweriing  3  case.  The  average 
RCS,  represented  by  the  variable  in  Equation  3.4=1  is  1.2  m2.  The  average 

error  between  the  histogram  and  the  Sweriing  1  density  is  also  quite  good  at 
0.011. 


In  Figure  3.4-2  we  have  the  X-band  case  also  at  0®  elevation.  Once 
again  the  region  out  to  15®  on  either  side  of  the  peak  was  removed.  Only  the 
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Swerling  case  1  is  shown  here.  The  exponential  distribution  fits  about  as  well 
with  this  case  as  with  Rgure  3.4-1 . 


Rgure  3.4-2.  Histogram  and  exponetial  distribution  for  aircraft  RCS 
@0®  elevation.  X-band 


Next,  we  look  at  the  missile  cross  section.  We  use  the  monostatic  RCS  of 
the  missile  for  L-band,  TE  mode,  and  0®  elevation  in  Rgure  3.1.2.2.6-1  (a). 
Rgure  3.4-3  shows  the  histogram  and  exponential  distribution  fit  for  this  case.  It 
is  apparent  that  the  fit  is  not  as  gccd  as  with  the  aircraft  case.  The  average  error 
per  point  is  now  0.04.  Only  the  Swerling  1  case  is  shown;  the  Swerling  3  case 
is  even  worse. 


101 


fvi  U5nuioncs.krt3.oar_' 


Figure  3.4-3.  Histogram  and  exponential  distribution  for  missile  RCS 
@0®  elevation.  L-band  -  TE  mode 


In  searching  for  another  probability  distribution  for  the  missile  RCS,  the 
log-normal  distribution  has  been  suggested  There  is  no  theoretical  model  of 
target  scattering  that  leads  to  the  log-normal  distribution,  although  it  has  been 
suggested  that  echoes  from  some  satellite  bodies,  cylinders,  plates,  and  arrays 
can  be  approximated  by  a  log-normal  probability  distribution^^^  This  distribution 
(density)  can  be  expressed  as 


p(o)  = 


a>  0 


where  s^  =  standard  deviation  of  In  (c/am),  and  Cm  =  median  of  c.  This  form  of 
the  log-normal  distribution  can  easily  be  written  in  the  standard  form  of  the 
normal  distribution  as 
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where  ^  =  In  a,  and  ^  =  In  Cm  -  s<j  =  In  Cp. 

Figure  3.4-4  shows  the  result  of  applying  the  log-normal  distribution  to 
the  missile  PCS.  The  average  error  per  point  is  now  only  0.018,  much  better 
than  the  exponential  distribution.  Even  though  it  seems  there  are  large  areas  of 
error  between  the  histogram  and  the  distribution,  we  must  keep  in  mind  that  the 
x-axis  is  plotted  here  in  terms  of  dBsm,  or  10  logio  Por  comparison,  most  of 
the  error  in  Figure  3.4-3  is  in  the  region  of  0.5  <  o  <  6  m^,  which  in  dBsm 
corresponds  to  -3  <  o  <  7.8  dBsm.  Figure  3.4-4  shows  that  in  this  region  the 
log-normal  distribution  is  much  closer  to  the  histogram  than  the  exponential 
distribution.  The  large  error  region  in  the  left  tail  of  the  distribution  is  not 
significant  enough  to  give  a  large  average  error.  We  note  also  that  the  mean 
PCS  for  the  missile  In  this  case  is  -8.1  dBsm,  which  is  what  we  expect  for  the 
missile. 


Rgure  3.4-4.  Histogram  and  log-normal  distribution  for  missile  PCS 
@0®  elevation.  L-band  -  TE  mode 
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In  Figure  3! 4-5  we  fit  the  log-normal  distribution  to  the  missile  RCS  at  X- 
band,  TM  mode,  and  0®  elevation  of  Figure  3.1.2.2.6-2  (b).  The  log-normal  fits 
the  histogram  very  well,  and  the  average  error  per  point  is  only  0.005.  The 
mean  RCS  Is  now  -6.8  dBsm.  The  good  fit  of  the  log-normal  distribution  makes 
it  a  good  probability  distribution  to  use  for  the  missile. 


Figure  3.4-5.  Histogram  and  log-normal  distribution  for  missile  RCS 
@0®  elevation.  X-band  -  TM  mode 


Naturally,  the  log-normal  distribution  is  also  applied  to  the  interaction 
signals.  We  present  two  cases  in  Rgures  3.4-6  (a)  &  (b).  The  first  case  is  the 
path  1  interaction  signal  of  the  air-to-air  monostatic  scenario  of  Figure  3. 1.3. 5-2 
(b),  and  the  second  case  is  the  path  3  interaction  signal  of  the  air-to-air  bistatic 
scenario  of  Rgure  3.1.3.5-10  (b).  In  both  cases  the  log-normal  distribution  does 
a  good  job  of  fitting  the  histogram  with  the  average  error  per  point  about  0.025. 
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Rgure  3.4-6  (a).  Histogram  and  log-normal  distribution  for  interaction  path  1 
Mono,  radar  pointing  @-45®  elevation.  L-band/TE  mode 


LonuvMMitn«er_o«r>o«t 


INTERACTION  SIGNAL  (OeSM) 

Rgure  3.4-6  (b).  Histogram  and  log-normal  distribution  for  interaction  path  3 
Bistatic  radar  pointing  @-45®  elevation.  X-band/TM  mode 
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3.5  ALGORITHMS 


3.5.1  Co-Taroet  Interference  Reduction  Formulation:  Unified  Theory 

Mathematical  formulation  provide  algorithmic  understanding  and 
extentions  that  otherwise  would  be  difficult  to  deduce.This  section  derives  the 
optimal  techniques  for  cotarget  interference  rejection.  The  following 
paragraphs  show  that  the  optimal  detection  statistics  lead  towards  suppression 
or  cancellation  schemes.  Thresholds  for  Constant  False  Alarm  Control  are  also 
be  defined  as  byproducts.lt  should  be  noted  that  there  are  two  distinct 
characteristics  that  make  co-target  interference  rejection  different  from  adaptive 
noise  interference  suppression.  These  are: 

1 .  Insufficient  number  of  samples  limit  any  hope  of  direct  adaptivity  and, hence,  a 
signal  model  must  eventually  be  assumed. 

2.  Sample  cross-covariance  of  sensor  degree  of  freedom  yield  high  correlation 
between  large  and  small  target.  Thus,  adaptive  cancellation,  as  in  Sample 
Matrix  Inversion  Techniques,  do  not  apply  because  both  large  and  small 
target  will  be  suppressed. 

Because  of  these  two  reasons.  Adaptive  Processing  Theory  is  not  directly 
applied  to  reject  large  targets.  However,,  adaptive  techniques  are  still  required 
for  ECM  ,  clutter,  and  thermal  noise  suppression  which  is  shown  to  be  part  of 
the  Whitening  Process  of  the  co-target  interference  rejection  algorithm. 

Our  approach  is  to  use  statistical  decision  theory  and  derive  the  optimal 
solutions  based  on  likelihood  functions.  The  model  (1)  assumes  combined 
simultaneous  space,  time  (doppler),  and  range  observations,  (2)  includes  ECM 
and  noise  processes,  (3)  the  presence  of  clutter ,  and  (4)  multiple  interracting 
targets.  A  maximum  likelihood  test  statistic  is  derived  for  all  Swerling  target 
models.  This  general  approach  establishes  the  optimal  performance  and  the 
algorithm  structures  which  can  be  modified  to  account  for  design  constraints 
compromises  Figure  3.5.1 -1  shows  mathematical  formulation  for  small  target 
detection.  This  formulation  assumes  that  in  a  small  range  doppler  region,  two 
targets  can  be  present  and  possibly  interfering  with  each  other.  Range,  doppler 
and  angle  are  unknown  parameters.  Without  loss  of  generality,  the  formulation 
can  be  extended  to  more  than  two  targets.  Our  problem  focuses  on  co-target 
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interference  in  a  radar  with  adequate  energy-aperture  products  to  detect  the  low 
RCS  target.  The  challenge  is  to  decide  between  hypothesis  H.|  and  Hg  and 

determining  performance  of  such  a  decision  rule. 


Hypothesis  Set 

H( 

H- 

H/ 

interference 
where 


y(t)=l(t) 

interference  only 

y(t)  =  2iSi(t)  +  l{t) 

signal  and  interference 

y(t)  =  2i  Si  (t)  +  22S2(^)  +  i(t)  • 

two  signals  and 

l(t)  = 

n(t) 

fi(t) 

JK(t) 

fln 

21.22 

Si(t) 

S2(t) 

S(t) 


model 


•kW 

an 

I 

Rk 


n(t)+£(t)+ X  1  j.(0 
k=1  ^ 

is  thermal  noise  vector  with  independent  components 
is  the  clutter  process  as  observed  from  each  spatial  channel 

is  the  jammer  process 

gain  and  phase  vector  associated  with  jammer  process 

observations  from  each  channel 
signal  1  and  2  gains  from  each  of  the  spatial  channels 

aiS(t-ti{t)) 

a2S(t-t2(t)) 

is  a  replica  of  the  transmitted  signal 

“  (Rj^  +  Rj^t  + ...) 

deterministic  or  random  variable  depending  on  target  Swelling 

is  the  k^^  target  range 
is  the  k^^  target  doppler 


Figure  3.5.1 -1 .  Hypothesis  Testing  Formulation  for  Small  Target  Detection  with 
the  Presence  of  Large  Interfering  Target,  Clutter,  Noise  and  Jammers. 


To  simplify  our  notation,  one  can  equivalently  consider  a  super  vector 
with  component  indices  being  samples  of  space  range,  and  doppler  domain. 
Thus,  one  defines  the  vector  SjLas  the  large  signal  vector  and  £2  as  the  small 
target  vector.The  first  derivation  that  follows  assumes  that  no  apriori  information 
is  used  .  One  can  basically  eliminate  the  noise  hypothesis  since  large  target 
presence  competing  with  noise  can  be  declared  with  high  probability. For 
Swerling  I  through  IV  target  model,the  log  likelihood  function  can  be  written  in 
terms  of  a  quadratic  functional  given  below 


-In  L  (y)  = 


y 


•1  -1 

(At  -A2)y  +  ln 


> 

< 


T 


where  Ag  _  is  the  noise,  clutter,  and  ECM  interference  covariance  matrix  and 


Ai  —  Ag  +  A2  —  Ai  +  S2S2 


The  ratio  of  the  two  determinants  can  be  shown  as 


A2 


1  +  SjAi^S 


The  use  of  the  Matrix  Inversion  Lemma 
the  two  covariance  inverses 


=  ^~0 


+  A-1  ^ 


\  + 


§1  Aq 


I-hS^Ao’St 


yields  the  following  relationship  among 


Define  (a,b)  =  a"^  Ag^  b 

as  an  inner  product  between  vector  a  and  ]2_,  then  the  likelihood  terms  can 
be  rewritten  as 
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SjA7’S2=|S2f-ii^ 

1+  §1 


H  +  S^A^Sg  =  1+  §2 


52*il)| 


1+  Si 


Eq.  (a) 


y^(Ai^-A2^)y  = 


y,  S2)-(y,Si^^Si,  §2 


Si 


2  |(SlM 


1+  Si 


Eq.  (b) 


The  complete  likelihood  ratio  test  then  becomes, 


(y’-2)~(y’^i)(^i’i2) 


-^n  L(y)  = 


iHis2ir- 


l|2  l(-1  ■  ^2)1 


+  MI  +  IIS2 


1+  Si 


1+  Si 


Eq.  (c) 


II  (l^  11  11^  |(§1»§2)l  II  |}2 

f^onlSill  SJ - jj — "2 —  =  |2lj|  where  a,  is  the  component  of  the 


small  target  SS-which  is  orthogonal  to  the  large  target  ^ 

That  is,  for  large  target  sufficiently  above  noise,  the  likelihood  statistic  results  in 


1+lhlr 


Eq.  (d) 


109 


Just  as  in  the  standard  detection  theory,  the  equation  shows  that  whitening 
followed  by  a  matched  filter  is  the  proper  detection  rule  .  The  unique  difference 
of  this  optimal  co-target  interference  rejection  rule  is  instead  of  the  filter  being 
matched  to  the  original  small  target  signal,  the  filteris  being  .matched  to  the 
small  signal  component  which  is  orthogonal  to  the  large  signal.  It  also  follows 
that  detection  and  false  alarm  probabilities  can  be  computed  from  the  available 
standard  curves  for  all  the  Swerling  target  models.  This  is  computed  by 
modifying  the  signal  energy  contained  only  in  the  orthogonal  component  in  the 
signal  to  noise  power  ratio.The  first  formulation  concludes  that  large  target 
suppression  is  optimal. 

A  second  formulation  capitalizes  on  the  capability  of  reprocessing  the 
data  given  the  large  target  is  present.  Given  a  detected  large  target,  via 
standard  methods,  one  would  like  to  search  for  small  targets  in  the  large  target 
vicinity.  Again,  only  two  hypothesis  are  considered.The  log  likelihood  statistic  is 


-tn  L(y)  =  (y  -  §,)  *  (aV  -  Ai’)(y  -  §,)  + 
where 

Ai  —  Aq  ^  A2  —  Aq  +  S2  §2 

Again  the  determinant  ratio  is 
given  as 


N 

IA2I 


1  +  Ao^  §2 


1-I 


-1 


Analogous  to  the  previous  derivation,  the  Matrix  Inversion  Lemma  yields, 

A”**  —  —2  _  ^0  §2  ~2  ^0^ 

1  +  S2Ao^S2  1  +  ||S2f 

resulting  in  the  following  log  likelihood  function 
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-^n  L(y)  = 


\{y-h^hf , 


1+  S 


-1 


The  above  equation  indicates  that  cancellation  algorithm  is  optimal.if  the 
presence  of  the  large  target  is  assumed  apriori. 

The  results  from  the  two  formulations  indicate  that  either  suppression  or 
cancellation  is  optimal  and  that  each  requires  a  signal  model.  Once  the 
orthogonal  or  residual  component  is  used,  one  applies  the  standard  detection 
techniques  and  the  Newman-Pearson  Lemma  lays  the  foundation  for  the  CFAR 
algorithms. 
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3.5.2  Description  and  Performance  Evaluation 


This  section  describes  and  evaluates  the  performance  of  several 
cancellation  algorithms  that  have  the  potential  to  reduce  co-target  interference. 
Four  effective  algorithms  have  been  identified.  They  are: 

1 .  Range  Response  Cancellation 

2.  Zero-Out  Blanking 

3.  Range  Sidelobe  Suppression 

4.  Doppler  Cancellation  in  Time  Domain 

The  first  three  techniques  are  in  the  range  domain  while  the  last  is  in  the 
doppler  domain.  Two  other  techniques,  doppler  cancellation  in  the  doppler 
domain  and  lattice  filter  cancellation  were  also  evaluated,  but  their  performance 
was  inadequate. 

Range  Domain  Issues.  Cancellation  in  the  range  domain  is  probably  the 
most  challenging  since  it  is  very  difficult  to  control  errors  that  affect  the  pulse 
compression  sidelobe  response.  The  main  error  factors  are  target  complex 
scattering  structure,  transmit  code  fidelity  video  filter  distortion  and  doppler 
mismatch  errors  (see  system  errors  section  3.5.3),  which  makes  difficult  the 
prediction  of  pulse  compression  sidelbobes.  Knowledge  of  the  sidelobe 
structure  is  required  if  the  techniques  are  to  be  fully  effective. 

The  first  approach  estimates  the  large  target  complex  amplitude,  range 
and  doppler,  predicts  the  sidelobe  contributions  to  other  bins  from  the  data  with 
the  response  model,  and  subtracts  these  directly  from  the  compressed  data. 
Alternatively,  to  minimize  processor  finite  word  length  distortions,  one  can 
subtract  the  large  target  from  the  uncompressed  data  and  then  recompressing 
the  neighboring  range  samples. 

3.5.2.1  Range  Response  Cancellation 

A  simple  suboptimal  technique  that  illustrates  range  domain  cancellation 
using  sidelobe  response  is  presented  below.  Consider  a  13:1  Barker  pulse 
compression  code  along  with  its  range  response  at  the  output  of  a  matched 
pulse  compressor  shown  in  Figures  3.5. 2.1-1  (a)  and  (b).  Given  the  large  target 
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range  estimate,  one  determines  the  range  samples  with  the  common  voltage 
values  at  the  output.  For  example,  range  samples  with  common  voltage  V2 
occur  at  times  t.12.  t^,  t^,  U>  and  ti2.  If  one  divides  the  peak  voltage  Vq  by  1 3 
and  subtract  the  values  at  these  sample  times,  the  large  target  is  cancelled. 
The  large  target  range  estimate  is  used  to  choose  the  set  of  range  bin  samples 
and  the  indexing.  Table  3.5.2. 1-1  lists  the  sidelobe  samples  with  common 
voltages  from  a  large  target  response  and  Barker  sidelobe  and  mainlobe 
relationships.  These  relationships  are  used  for  cancellation.  Figure  3.5.2.1-1 
(a)  also  shows  a  radar  return  pulse  compression  code  after  video  filtering.  The 
transitions  between  ±1 80®  phases  are  most  notably  affected  as  shown  by  the 
dashed  curve.  This  can  limit  performance  since  one  must  account  for  the 
smoothing  of  the  code.  Furthermore,  video  filters  are  typically  not  well  matched 
from  channel  to  channel  and  system  to  system,  and  their  responses  can  vary 
with  time.  Therefore,  periodic  calibration  may  bo  necessary  for  such  a 
cancellation  scheme  to  yield  the  required  performance. 


TABLE  3.5.2.1-1.  RANGE  SAMPLE  PROPERTIES  OF  13:1  BARKER  PULSE 
COMPRESSION  CODE  RESPONSE  FOR  A  TARGET  STRADDLED  IN  RANGE 


VOLTAGE 

CONSTRAINTS 

RANGE  SAMPLE  TIMES  WITH 

COMMON  VOLTAGE  AFTER  P.C. 

CODE 

Vi 

Li3>  Lg-  t-s-  tg-  V-  2^  ^i  i 

13  Vq  “  V2 

V2 

t-1 2'  *-8'  '-4*  U'  *8’  ^1 2 

13  V-1  *  V, 

•Vi 

t-l  1  •  t.7>  t.3>  *1  >  Ig 

and 

f-l  O'  f-6’  f-2’  '2'  *6'  ti  0 

lax 
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BEFORE  VIDEO  FILTER 
AFTER  VIDEO  FITER 
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(a)  13  T0 1  BARKER  BINARY  PHASE  CODE 


(b)  VOLTAGE  RESPONSE  OF  13:1  PULSE  COMPRESSION  (BARKER  CODE)  SAMPLE  TIMES  CORRESPONDS 

TO  A  STRADDLED  TARGET, 


Figure  3.5.2. 1-1  (a)  &  (b).  Barker  Pulse  Compression  Code  Properties. 

Figure  3. 5.2. 1-2  (a)  &  (b)  show  the  functional  block  diagrams  for  the 
normal  pulse  doppler  radar  system  and  the  range  ressponse  cancellation 
algortihm.  Pulse  oompression  is  done  right  after  doppler  filtering.  In  Figure 
3.5.2. 1-2  (b),  the  THRESHOLD  block  on  top  represents  the  detection  of  the 
large  target  peak  response  location  (mainlobe)  and  response.  Then  1/13  of  the 
peak  is  subtracted  from  the  adjacent  range  bins  which  predominantly  contain 
the  sidelobes  of  the  large  target.  This  is  done  in  the  CANCEL  block.  Although 
not  shown,  the  responses  of  the  large  range  bin  and  the  next  olosest  range  bin 
(before  or  after,  depending  upon  the  straddling)  are  set  to  zero.  Finally,  the 
residual  response  is  passed  through  a  second  THRESHOLD  for  detection  of 
small  targets. 
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(a)  PULSE  DOPPLER  PROCESSING:  NO  CANCELLATION 
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TARGETS 


(b)  CANCELLATION  USING  PC  RESPONSE 
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TARGETS 


SMALL 

TARGETS 


Figure  3.5.2. 1-2.  Functional  Block  Diagram  for  (a)  Standard  Pulse  Doppler 
Processing  (b)  Range  Domain  Cancellation  Using  the  Response 


The  performance  of  this  cancellation  technique  is  examined  next.  Figure 
3.5.2. 1-3  (a)  shows  the  case  where  only  the  large  target  is  present  at  range  bin 
15.  Only  the  doppler  filter  where  the  large  target  is  present  is  shown.  Also, 
system  errors  such  as  hardware  errors,  timing  control  errors,  etc.  are  not  taken 
into  account  in  all  the  following  performance  plots.  This  will  simplify  the  task  of 
ascertaining  the  effectiveness  of  the  different  cancellation  algorithms.  The 
expected  sidelobe  response  of  the  Barker  code  can  be  clearly  seen,  even 
though  a  small  amount  of  noise  is  present.  We  note  that  in  the  no-noise  case 
the  sidelobes  should  be  -22.3  dB  down  from  the  peak,  which  the  sidelobes  in 
the  figure  generally  are  about  that  value.  Figure  3.5.2.1-3  (b)  shows  the  same 
doppler  filter  after  cancellation  has  been  performed.  The  large  target  peak  at 
range  bin  15  has  been  effectively  cancelled,  while  the  sidelobes  have  also 
been  reduced  about  15  to  20  dB. 
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Figure  3.5.2.1-3  (a).  Pulse  compression  response  before  cancellation 

One  target  case. 
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Figure  3. 5.2. 1-3  (b).  Response  after  cancellation.  One  target  case. 


Figures  3.5.2. 1-4  (a)  &  (b)  show  the  most  important  case  of  a  large  target 
interfering  with  a  nearby  small  target.  The  small  target  at  range  bin  1 7  is  only  2 
range  bins  away  from  the  large  target.  Also,  its  power  is  20  dB  below  that  of  the 
large  target  peak.  The  two  targets  are  assumed  to  have  the  same  velocity  so 
they  will  be  in  the  same  doppler  filter  for  simplicity.  As  Figure  3.5.2.1-4  (a) 
shows,  the  s.mail  target  peak  (-20  dB)  and  the  large  target  sidelobe  at  range  bin 
17  (-22.3  dB)  combine  to  produce  the  response  in  that  range  bin  of  about  -16 
dB.  After  cancellation,  the  small  target  dearly  stands  out  in  range  bin  17,  and 
as  shown  in  Rgure  3.5.2. 1-4  (b)  its  response  is  indeed  20  dB  down  from  the 
large  target  peak,  as  expected.  As  this  peak  stands  out  so  clearly,  it  can  be 
easily  detected  and  a  small  target  is  declared.  Hence,  the  range  response 
cancellation  scheme  is  seen  to  be  very  effective  for  target  separation  as  dose 
as  2  range  bins. 
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LARGcTARGci  (a)  RB^15  SMALL  TARGET  @  RBifM 

Figure  3.5.2.1 -4  (b).  Response  after  cancellation.  Two  targets  case. 

To  snow  that  the  inter-target  separation  should  be  at  least  2  range  bins, 
we  lock  at  Figures  3.5.2. 1-5  (a)  &  (b),  which  depict  smaller  separations.  In  (a), 
the  separation  is  1.5  range  bin  for  a  small  target  20  dB  below  the  large.  After 
cancellation,  the  small  target  remains  at  range  bin  17  ( it's  been  shifted),  but  it  is 
not  much  above  the  sidelobe  residuals  of  the  large  target.  The  problem  here  is 
that  the  small  target  energy  in  range  bin  17  is  about  26  dB  below  the  peak, 
which  is  6  dB  worse  than  the  case  of  Figure  3.5.2. 1-4  (b).  Hence,  it  would  be 
diracult  to  detect  this  small  target.  Next,  in  Figure  3.5.2. 1-5  (b),  the  small  target 
is  only  1  range  bin  away.  As  can  be  seen,  after  cancellation,  it  was  cancelled 
out  along  with  the  large  target,  and  therefore  1  range  bin  separation  is  too 
small. 
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igure  3. 5.2.1 -5  (a).  Response  after  cancellation.  Two  targets  case. 
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Figure  3. 5.2.1 -5  (b).  Response  after  cancellation.  Two  targets  case. 
1  range  bin  separation. 
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Next,  we  present  the  case  where  the  small  target  power  is  40  dB  below 
that  of  the  large  target.  Figure  3.5.2. 1-6  shows  this  case.  The  small  target  is  2 
range  bins  away.  After  cancellation  the  small  target  actually  gains  some 
energy,  as  its  response  in  range  bin  17  is  37  dB  below  the  large  target  peak, 
while  we  expect  the  response  to  be  40  dB  below.  The  explanation  is  that  not  all 
the  large  target  sidelobe  energy  in  range  bin  17  was  cancelled. 


LARGE  TARGET  (3  RB#1 5  SMALL  TARGET  <§>  RB#!  7 

Rgure  3.5.2. 1-6.  Response  after  cancellation.  Two  targets  case. 

40  dB  dynamic  range. 

Next,  we  look  at  the  case  where  the  large  target  has  2  components:  one 
component  centered  at  range  bin  15  and  the  other  component  straddling  at  1/3 
range  bin  away  with  each  component  having  half  the  power.  The  small  target  is 
3  range  bins  away  and  is  20  dB  down.  Figure  3.5.2. 1-7  (a)  shows  the  response 
before  cancellation  and  in  (b)  we  have  the  response  after  cancellation.  As  the 
Figure  shows,  the  cancellation  is  quite  effective,  and  the  small  target  is  easily 
detected. 
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Figure  3.5.2. 1-7  (a).  Response  before  cancellation.  Two  targets  case. 
2-oomponents  large  target. 
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Figure  3.5.2.1-7  (b).  Response  after  cancellation.  Two  targets  case. 

2-compcnents  large  target. 
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3.5.2.2  Zero-Out  Blanking 


Blanking  techniques  to  remove  large  interfering  signals  have  been  used 
extensively  in  radars.  Algorithms  utilizing  spatial,  range  and  doppler 
observations  offer  a  simple  effective  method  of  maintaining  detection  sensitivity 
without  significantly  increasing  processing.  Blanking  effectiveness  decreases, 
however,  as  target  separation  decreases.  This  occurs  because  with  small 
separations,  masking  increases  and  blanking  removes  the  large  target  cells 
together  with  a  significant  portion  of  the  small  target.  The  classical  blanking 
techniques  maintains  detection  sensitivity  by  not  allowing  large  targets  to 
influence  the  threshold  estimates  in  a  constant  False  Alarm  Rate  (CFAR) 
system.  Given  a  large  target  detection  with  standard  processing  in  the  range- 
doppler  cell  domain,  the  technique  basically  sets  to  zero  the  sidelobe  and 
mainlobe  cells  containing  large  target  energy.  Then,  a  new  local  threshold  is 
estimated  to  which  the  neighboring  cells  are  compared.  With  very  low  doppler 
filter  sidelobes,  it  may  not  be  necessary  to  blank  sidelobe  doppler  cells  except 
for  the  two  cells  nearest  the  mainlobe  (depends  on  filter  broadening  due  to 
heavy  amplitude  weighting).  It  is  clear  that  that  blanking  is  undesirable  in  the 
doppler  domain  since  it  completely  desensitizes  the  full  PRF  for  the  large  target 
range  bin. 

A  variation  of  range  domain  blanking  has  been  developed  that 
progressively,  rather  than  abruptly,  degrades  detection  sensitivity  as  the  small 
target  is  closer  to  the  large  target.  After  large  target  detection,  the  set  of 
uncompressed  large  target  range  samples  in  the  target  doppler  filter  are  set  to 
zero  and  then  the  set  of  (2Npc  +  1)  of  range  samples  are  recompressed. 
Figure  3.5. 2. 2-1  shows  the  functional  block  diagram  for  the  zero-out  blanking 
algorithm.  As  the  diagram  shows,  after  the  large  target  is  detected,  its  location 
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Figure  3.5.2.2-1.  Zero-Out  Blanking  in  the  Range  Domain  Functional  Block 

Diagram. 

is  fed  into  the  BLAMKING  block.  Here,  the  uncomipressed  Npc  large  target 
range  samples  are  blanked  (set  to  zero).  Pulse  compression  is  done  once 
again  and  the  result  is  passed  through  thresr.clcing  for  smalt  target  detection. 

The  eftectiveness  of  the  zero-cut  blanking  is  next  determined.  Figure 
3.5.2.2-2  (a)  &  (b)  shew  the  case  where  the  small  target  is  5  range  bins  from  the 
large  target.  As  range  bin  20  is  in  one  of  the  nulls  of  the  large  target  response, 
when  it  is  combined  with  the  small  target  peak  there  (20  dE  down),  the  resultant 
response  is  19.6  cE  down  from  the  large  target  peak.  This  is  shewn  in  Figure 


LARGS  TARGST  -a  RESl 5  SMALL  TARGET  (g>  Re«20 

Figure  3. 5. 2.2-2  (a).  Pulse  compression  response  before  blanking. 
Two  target  case.  5  range  tin  separation. 
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3.5.2.2-2  (a).  After  the  blanking  algorithm  is  applied,  in  (b).  the  response  in 
range  bin  20  is  only  -29.2  dB,  or  a  loss  cf  about  9.2  dB  of  its  own  power. 
However,  the  small  target  energy  remaining  in  range  bin  20  is  sufficiently  above 
the  residual  in  the  surrounding  range  bins  to  allow  a  detection  of  the  small 
target. 


Figure  3.5.2.2-2  (b).  Response  after  blanking.  Two  target  case. 

5  range  bin  separaticn. 

In  Rgures  3.5.2.2-3  (a)  &  (b)  we  mcve  the  small  target  to  range  bin  19  so 
that  it  is  slightly  closer  to  the  large  target,  4  range  bins  away.  In  (a)  we  see  that 
since  range  bin  19  is  in  one  of  the  sidelcbe  peaks  of  the  large  target,  when  it  is 
combined  with  the  small  target  (20  dB  down)  the  resultant  response  is  15.1  dB. 
After  the  blanking  algorithm  is  applied,  in  (bj,  the  response  in  range  bin  19  is 
about  32.7  dB  down  from  the  large  target  peak.  Hence  the  small  target  has  lest 
about  12.7  dB  cf  its  power.  The  small  target  power  remaining  is  about  the  same 
as  the  residual  in  some  surrounding  range  tins.  It  would  be  difficult  to  declare  a 
detection  given  this  amount  cf  power.  Thus  given  a  4  range  bin  separaticn,  this 
is  tco  close  for  the  zero-out  blanking  algorithm  to  be  effective. 
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Figure  3.5.2.2-3  (a).  Pulse  compression  response  before  blanking. 
Two  target  case.  4  range  bin  separation. 
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Figure  3.5.2.2-3  (b).  Response  after  blanking.  Two  target  case. 

4  range  bin  separation. 
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At  5  range  bin  separation  and  40  dB  dynamic  range,  the  zero-out 
blanking  is  still  effective.  This  is  shown  in  Figure  3.5.2.2-4.  The  response  in 
range  bin  20  is  the  small  target  residual,  about  48  dB  down  from  the  large  target 
peak.  Compared  to  the  surrounding  range  bins,  however,  it's  still  enough  to 
declare  the  detection  of  a  small  target. 


Figure  3.5.2.2-4.  Response  after  blanking.  Two  target  case. 
5  range  bin  separation.  40  dB  dynamic  range. 
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3.5.2.3  Range  Sidelobe  Suppression 


In  radar  applications,  pulse  compression  waveforms  are  frequently  used 
to  increase  the  amount  of  energy  radiated,  and  a  pulse  compressor  is  required 
to  process  the  received  signal.  For  these  applications,  the  compressor  is  often 
an  implementation  of  a  matched  filter  (maximizes  detection  range)  whose  range 
sidelobes  can  be  relatively  high.  If  one  assumes  that  only  a  single  target  is 
present,  then  one  can  reject  those  detections  due  to  the  high  sidelobes,  and 
false  targets  are  eliminated.  However,  if  the  environment  contains  both  large 
and  small  targets,  all  the  sidelobes  from  the  large  target  will  interfere  with  the 
detection  of  potential  small  targets  in  the  range  region  where  the  sidelobes  are 
present.  Hence,  for  the  scenario  where  large  and  small  targets  are  present,  the 
control  of  the  sidelobes  is  extremely  important. 

Sidelobe  control  for  pulse  compression  waveforms  takes  on  many  forms. 
For  digital  codes,  this  is  achieved  by  forming  a  mismatched  filter,  as  opposed  to 
a  matched  filter  whose  filter  weights  specify  the  response.  There  are  algorithms 
to  minimize  integrated  sidelobe  ratio  (ISLR)  or  peak  sidelobe  ratio  (PSLR),  but 
both  approaches  extend  the  length  of  the  mismatch  filter  beyond  the  length  of 
the  matched  filter.  For  air-to-air  environment  where  sensitivity  is  important, 
extending  the  weight  length  extensively  will  induce  higher  SNR  losses  from  the 
weights.  An  acceptable  approach  is  one  which  uses  a  matched  filter  to  initially 
detect  targets,  then  re-processes  the  data  to  examine  regions  where  the 
sidelobes  from  the  initial  detections  hindered  sensitivity  using  sidelobe  control 
techniques.  The  following  discusses  such  an  approach  to  the  sidelobe  control 
problem  and  gives  some  examples. 

One  approach  to  control  sidelobes  is  to  re-process  the  data  for  a  set  of 
range  samples  with  a  weighting  sequence  that  provides  large  target 
suppression.  Figure  3.5.2.3-1  depicts  the  desired  properties  of  the  resultant 
response  using  the  weights.  The  response  has  very  low  sidelobes  in  the  region 
of  the  large  target  and  for  all  ranges  in  the  set  of  range  samples.  The  data  must 
be  re-processed  for  each  range  region  as  shown  in  the  figure.  It  is  assumed 
that  the  large  target  is  not  eclipsed  in  range,  since  one  can  always  find  a  PRF 
that  insures  that  the  large  target  will  not  be  eclipsed  once  the  range  is  known. 
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REGION  TO  EXAMINE 
POTENTIAL  TARGETS 


Figure  3.5.2. 3-1.  Placement  of  Low  Sidelobes 


The  weighting  sequence  that  gives  a  specific  low  sidelobe  response  can 
be  easily  constructed  by  using  the  Gram-Schmidt  procedure  to  orthogonalize  a 
set  of  vectors.  Let  c(n),  n=1 ,  be  the  sequence  which  specifieds  the  digital 
code  whose  length  is  Ng.  For  a  matched  filter,  the  weights  are  c(n).  For  the 
mismatched  filter,  begin  by  forming  the  vectors  given  by 


c(n-k) ,  n  =  k+1,...,  k+Nc 
0  ,  otherwise 


which  are  the  various  shifts  of  the  signal  where  jso  'S  fhe  assumed  shift  of  the 
signal  where  a  potential  target  is.  For  a  set  of  weights  given  by  a  weight  vector 
W.  the  response  r(k)  at  shift  k  is  given  by 

r(k)  =  Mt*Xk 


where  r(0)  is  the  mainlobe  response.  The  weight  vector  is  then  selected  such 
that  r(k)  is  zero  for  k  in  the  range  K,  and  K2  (one  of  the  range  regions),  the  sum 

of  the  weights  is  unity  and  r(0)  is  maximized.  This  criterion  will 

select  the  weight  vector  such  that  the  sidelobes  are  zero  in  the  desired  range 
region  and  the  signal-to-noise  loss  due  to  the  mismatched  weights  is 
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minimized.  With  these  definitions  the  Gram-Schmidt  procedure  can  be  used  to 
form  a  set  of  unit  vectors  Sj,  i=1  which  span  the  set  of  vectors  k=Ki,  K2. 
Then  the  weight  ^  is  found  by  extending  the  Gram-Schmidt  procedure  one 
stage  further  and  finding  the  unit  vector  ^  which  is  orthogonal  to  all  fij,  i=1 
and  minimizes  the  quantity  w*)^.  Hence  the  desired  weight  vector  is  given  by 


Xo 


i.1 


i.1 


Figure  3.5.2.3-2  shows  the  functional  block  diagram  for  the  range 
sidelobe  suppression  algorithm.  As  the  diagram  shows,  after  the  large  target  is 
detected,  its  location  is  utilized  in  the  weight  computation  procedure.  Then  the 
FFT  data  is  recompressed  and  multiplied  by  the  weight  sequence  to  null  out  the 
sidelobes.  Although  not  shown  in  the  figure,  the  range  bin  containing  the  large 
target  peak  is  zeroed  out.  Also,  in  the  implementation  of  the  suppression  algo- 
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Figure  3.5.2.3-2.  Range  Sidelobe  Suppression  Functional  Block  Diagram 


129 


rithm,  the  region  to  search  for  potential  small  targets  (between  the  K,  and  K2 
sidelobes)  covers  the  entire  sidelobe  region  of  the  large  target  on  both  sides  of 
the  large  target  peak.  Hence  all  12  sidelobes  of  the  large  target  response  will 
be  suppressed. 


The  performance  of  the  range  sidelobe  suppression  algorithm  is  similar 
to  the  performance  of  the  range  response  cancellation  looked  at  before.  Rgure 
3. 5. 2. 3-3  below  is  the  response  after  the  suppression  algorithm  has  been 
applied  for  the  case  where  only  1  target  is  present  at  range  bin  15.  The 
response  before  suppression  has  been  shown  before  in  Figure  3.5.2.1-3  (a). 
As  Figure  3. 5.2. 3*3  shows,  the  target  peak  at  range  bin  15  is  completely 
suppressed,  while  all  the  sidelobes  are  minimized  as  much  as  possible. 
Hence,  if  there  is  only  one  large  target  present,  range  suppression  will  not 
create  any  false  alarm. 


Figure  3.5.2.3-3.  Response  after  range  suppression.  One  target  case. 

The  two-target  case  is  shown  next  in  Rgure  3.5.2.3-4.  The  small  target  is 
20  dB  below  the  large  target  peak  and  is  2  range  bins  away  from  the  large. 
Range  suppression  clearly  reveals  the  small  target  at  range  bin  17  with  the 
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residual  response  there  slightly  less  than  20  dB  down.  It  would  be  easy  to 
detect  this  small  target  since  it  towers  over  the  residuals  in  surrounding  range 
bins.  The  small  target  has  lost  about  0.88  dB  of  its  own  power.  This  is  slightly 
more  than  the  signal  loss  after  range  response  cancellation  for  the  same 
scenario  (Figures  3.5.2. 1-4  (a)  &  (b)),  when  the  small  target  loss  was  0.34  dB. 
Thus,  the  performance  of  range  suppression  and  range  cancellation  are  quite 
similar  with  range  cancellation  slightly  more  effective. 
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Figure  3.5.2.3-4.  Response  after  range  suppression.  Two-  target  case. 

2  range  bins  separation. 

With  40  dB  dynamic  range,  the  range  suppression  algorithm  is  somewhat 
effective.  Rgure  3.5.2.3-5  shows  the  result  for  a  small  target  also  at  2  range 
bins  away  from  the  large  target  but  with  power  40  dB  down.  The  region 
between  range  bin  2  and  28  contains  the  sidelobe  response  of  the  large  target 
before  suppression.  After  suppression,  the  residual  response  in  range  bin  17  is 
of  the  small  target,  which  is  about  41  dB  down,  and  this  is  what  we  expect. 
However,  as  seen  in  the  Rgure,  the  residual  sidelobe  next  closest  to  the  large 
target  peak  at  range  bin  13  is  quite  significant,  and  can  create  a  false  alarm. 
This  residual  cannot  be  made  smaller,  and  therefore,  if  the  dynamic  range  is 
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greater  than  40  dB,  the  residual  will  dominate  after  suppression,  thus  creating  a 
false  alarm. 


RANGE  BIN 

L^RGE  TARGET  @  RB#1S  SMALL  TARGET  @  RB#17 

Figure  3.5.2.3-5.  Response  after  range  suppression.  Two-  target  case. 
2  range  bins  separation.  40  dB  dynamic  range. 


3.5.2.4  Doppler  Cancellation  in  the  Time  Domain 

There  are  interesting  possibilities  for  cancellation  in  the  doppler  domain. 
Digitally  formed  doppler  filters  exhibit  excellent  predictable  response  functions. 
That  is,  since  coherent  radars  maintain  spectrally  clean  transmit  and  reference 
signals,  doppler  filter  sidelobe  control  with  digital  processors  are  quite  feasible. 
For  example,  airborne  tactical  radars  utilize  90  dB  sidelobe  Dolph-Chebychev 
Filters  to  minimize  clutter  and  co-target  interference.  The  disadvantages  with 
such  filter  response  is  the  mismatch  (better  known  as  amplitude  weighting)  loss 
and  resolution  loss  since  the  filter  mainlobe  broadens  in  comparison  to  a 
uniformly  weighted  filter.  Since  detection  of  low  RCS  target  requires  a  taxing 


energy-aperture  product,  processing  losses  associated  with  heavy  amplitude 
weighting  may  not  be  acceptable.  Given  a  processing  chain  with  a  moderate 
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amplitude  weighted  filter  bank,  a  promising  technique  of  cancellation  in  the 
doppler  domain  is  described  next. 

Figure  3.5.2.4-1  shows  the  functional  block  diagram  of  the  doppler 
cancellation  in  the  time  domain  algorithm.  Basically,  a  three-point  quadratic  fit 
scheme  Is  used  to  estimate  the  doppler,  phase,  and  magnitude  of  the  large 
peak.  Then  this  estimated  time-domain  data  of  the  largo  target  is  subtracted 
from  the  l/Q  data  in  the  range  bins  where  the  largo  target  is  detected.  This 
resulting  data  is  then  passed  on  down  through  the  processing  chain  including 
doppler  filtering  and  pulse-compression.  The  entire  procedure  is  repeated  as 
many  times  as  desired,  with  each  time  taking  the  largest  peak  and  reducing  it 
further. 
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Rgure  3.5.2.4-1 .  Doppler  Cancellation  in  the  Time  Domain  Functional  Block 

Diagram 


The  performance  of  this  technique  is  quite  good.  It  has  the  ability  to 
detect  a  small  target  even  when  the  target  separation  is  only  one  doppler  filter 
apart.  This  is  shown  in  Rgure  3.5.2.4-2  (a)  &  (b).  In  (b),  the  small  target 
definitely  pops  out  clearly  while  the  largo  target  has  been  reduced  by  about  50 
dB.  The  performance  is  also  similarly  very  good  when  the  target  separation  is 
1 .5  filters,  as  shown  in  Figure  3.5.2.4-3  (a)  &  (b). 
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Figure  3.5.2.4-2  (b).  Response  after  doppler  cancellation.  Two  targets. 
1  doppler  filter  separation.  40  dB  dynamic  range 
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Figure  3.5.2.4-3  (b).  Response  after  doppler  cancellation.  Two  targets. 
1 .5  doppler  filters  separation.  40  dB  dynamic  range 


3.5.2.5  Other  techniques 

We  describe  two  other  cancellation  techniques  that  were  presented  in 
the  proposal.  However,  the  description  will  be  brief  since  their  performance 
was  found  to  be  inadequate.  Figure  3. 5.2. 5-1  shows  the  functional  block 
diagram  for  the  doppler  cancellation  in  the  doppler  domain  technique.  In  this 
technique,  the  large  target  signal  is  detected,  and  then  it  is  subtracted  out  from 
the  doppler  domain  data,  which  is  the  output  of  the  FFT  block. 


Figure  3.5.2.5-1 .  Doppler  Cancellation  in  the  Doppler  Domain  Functional  Block 

Diagram 


This  technique  was  found  to  be  quite  ineffective  at  removing  large-target 
in  the  doppler  domain.  Figure  3.5.2.5-2  is  a  plot  in  the  doppler  domain  of  a 
case  where  the  large  target  is  originally  at  doppler  filter  105  and  the  small  target 
is  at  filter  108.  As  the  Figure  shows,  the  large  target  cannot  be  completely 
removed:ln  fact,  it  still  has  a  very  large  magnitude.  Furthermore,  this 
cancellation  scheme  shifts  the  targets  to  incorrect  doppler  filters,  creating  false 
alarms  that  are  difficult  to  reject.  Because  of  these  characteristics,  this 
technique  was  not  pursued  further. 
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Figure  3.5.2.5-2.  Doppler  Cancellation  in  the  Doppler  Domain.  Large  Target  <§> 


filter  #1 05.  Small  Target  <§>  filter  #1 08. 


Another  technique  is  a  doppler  domain  suppression  technique,  which 
concatenates  a  low-order  filter  before  the  doppler  filters.  The  low-order  filter 
places  a  null  in  the  doppler  domain  corresponding  to  each  larget  target.  The 
parameters  that  set  the  low-order  filter  nulls  can  be  determined  by  estimating 
target  doppler  and  power.  A  Lattice  Filter  is  an  excellent  candidate  structure  for 
the  low-order  filter  since  it  is  a  time  and  order  recursive  fiiter.  This  means  that  it 
has  a  recursive  formulation  with  time,  and  also  as  the  order  of  the  filter  is 
increased  the  adaptively  determined  reflection  coefficients  do  not  have  to  be 
recalculated  except  for  the  added  stage  which  can  be  calculated  from  the 
existing  coefficients.  Rgure  3.5.2.5-3  shows  such  a  filter,  which  will  be  followed 
by  the  standard  doppler  filters,  when  a  strong  target  is  detected.  The 
coefficients  are  determined  adaptively  prior  to  doppler  filter  formation. 
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Figure  3.5.2.5-3.  Implementation  of  a  Lattice  Filter. 

Figure  3. 5.2. 5-4  shows  the  functional  block  diagram  of  the  doppler 
suppression  scheme  using  the  Lattice  filter.  As  the  Figure  shows,  the  large 
target  information  is  used  to  determine  filter  coefficients  that  would  place  nulls  in 
the  doppler  domain.  The  entire  data  is  then  passed  through  the  FFT  doppler 
filter  bank  and  the  usual  processing  sequence. 

SUPPRESSION  USING  LOW-ORDER  LATTICE  FILTER 


IQ 

DATA 


LARGE 

TARGETS 


SMALL 

TARGETS 


Rgure  3.5.2.5*4.  Doppler  Domain  Suppression  with  Lattice  Filter 
Functional  Block  Diagram 


The  lattice  filter  technique  was  also  found  to  be  inadequate  at  removing 
large-target  interference.  It  is  not  able  to  remove  the  large  target  and  thus  there 
is  no  way  to  detect  the  small  target.  Hence,  the  doppler  domain  suppression 
with  lattice  filter  technique  was  not  pursued  further. 
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8g^3ie-13  (03-13-89) 


3.5.3  Error  Analysis 


There  are  many  error  factors  which  must  be  addressed  in  the  design  of 
the  cancellation  or  suppression  mode.  In  this  section,  a  methodology  for 
incorporating  these  factors  is  presented.  It  is  based  on  the  utilization  of  budgets 
to  establish  a  structured  design  approach,  which  also  allows  the  control  of  the 
design. 

To  begin  with,  many  of  the  factors  can  bo  broken  down  into  contributions 
which  reduce  cancellation  or  suppression  ability,  increase  false  alarm  (FA)  rate 
or  decrease  angular,  range  and  doppler  accuracy.  Hence,  it  is  natural  to  show 
the  effects  by  establishing  the  following  budgets: 

1.  Cancellation/Suppression 

2.  False  alarm 

3.  Angular,  range  and  doppler  estimation  error 

Besides  these  budgets,  the  normal  budgets  of  losses  and  S/N  contributors 
complete  the  system  specification.  From  these  budgets,  probability  of  detection 
and  estimation  performance  can  be  directly  evaluated. 

The  potential  error  sources  in  the  system  are  shown  in  Table  3.5.3-1. 
Our  approach  in  doing  the  error  sensitivity  analysis  is  composed  of  three  steps: 

1 .  Determine  the  small  target  signal-to-interference  ratio  (SIR)  without 
cancellation  in  the  perfectly  errorless  system,  and  without  the  large 
target. 

2.  Introduce  a  single  source  of  error  from  the  Table  into  the  system. 

3.  Determine  the  small  target  SIR  after  cancellation. 

4.  Compute  the  loss:  The  difference  in  SIR  between  step  3  and  step  1 
is  the  loss. 

The  small  target  SIR  is  taken  to  be  the  ratio  of  the  small  target  signal  at  the  peak 
response  to  the  estimated  interference  (thermal  noise  included)  in  the 
surrounding  range-doppler  window.  The  estimation  process  is  similar  to  CFAR, 
where  a  window  is  formed  around  the  range-doppler  cell  of  the  small  target, 
and  the  average  interference  power  per  cell  is  estimated.  The  interference 


139 


consists  mainly,  of  the  large  target  sidelobes,  which  although  reduced  by 
cancellation,  still  have  residues  in  the  small  target  range-doppler  cell. 

TABLE  3.5.3-1 .  TABLE  OF  POTENTIAL  ERROR  SOURCES 
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Table  3.5.3-2  shows  the  loss  associated  with  the  three  range  domain 
techniques  when  there  are  no  errors  in  the  system.  The  two  targets  are 
separated  by  4  range  bins  and  there  is  a  20  dB  dynamic  range  difference 
between  them.  The  loss  shown  in  this  Figure  is  the  difference  between  small 
target  SIR  after  cancellation  and  the  small  target  SIR  computed  in  step  1  above. 

TABLE  3.5.3-2.  SMALL  TARGET  LOSS  WITH  NO  ERRORS 


Technique 

Loss  (dB) 

Range  Response  Cancel. 

0.099 

Zero-Out  Blanking 

9.07 

Range  Suppression 

1.49 

This  Table  shows  that  in  the  errorless  system,  the  range  response  cancellation 
technique  results  in  negligible  loss  while  the  zero-out  blanking  technique 
reduces  the  small  target  signal  greatly.  Hence,  zero-out  blanking  entails  a  large 
penalty  each  time  it  is  used,  and  is  not  highly  recommended.  Range 
suppression  gives  good  performance  with  a  small  amount  of  loss. 

The  first  error  source  to  be  looked  at  is  the  effect  of  the  video  filter 
bandwidth.  The  cancellation  loss  associated  with  the  three  range  domain 
techniques  are  given  in  Figure  3.5. 3-1.  The  range  response  cancellation 
results  in  the  least  amount  of  loss  while  blanking  gives  the  most  loss.  It  is  seen 
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Vld«o  Bandwidth 


- Q—  Range  can. 

“ — • —  Blanking 
- ■ —  Range  suppression 

4  RB  separation 
20  dB  dyn.  range 


Figure  3.5.3-1.  Cancellation  loss  vs.  Video  Bandwidth 


that  small  video  bandwidth  results  in  large  losses  for  all  three  techniques,  and 
so  the  video  filter  must  be  designed  carefully.  In  the  cancellation  simulation,  the 
errorless  case  corresponds  to  when  the  video  bandwidth  is  greater  than  6. 

Next,  Figure  3.5.3-2  shows  the  effect  of  A/D  quantization  bits.  In  an 
actual  radar  processing  system,  the  l/Q  data  is  sampled  by  the  A/D  converter 
into  digital  data,  and  naturally  this  will  reduce  the  accuracy  of  the  data  and  the 
cancellation  performance. 
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Figure  3.5.3-2.  Cancellation  loss  vs.  A/D  Quantization  bits 

As  the  plot  shows,  the  loss  is  essentially  independent  of  the  number  of 
A/D  bits,  so  the  quantization  does  not  have  any  real  effect  on  the  performance. 
The  losses  are  at  the  same  level  as  in  the  case  where  there  are  no  errors  in  the 
system. 

The  next  error  to  be  looked  at  is  l/Q  Amplitude  Imbalance  error.  This 
error  affects  the  l/Q  data  by  increasing  the  amplitude  of  the  Q  component  of  the 
data.  Figure  3.5.3-3  shows  the  effect  on  the  cancellation  performance. 
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Figure  3.5.3-3.  Cancellation  loss  vs.  I/Q  Amplitude  Imbalance 

When  amplitude  imbalance  is  in  the  system,  the  range  response 
cancellation  actually  results  in  signal  gain  for  the  small  target,  as  evidenced  by 
the  gain  of  about  0.2  dB.  Somehow,  the  amplitude  imbalance  has  reduced  the 
interference  and  noise  surrounding  the  small  target,  which  increases  the  small 
target  SIR.  As  for  range  suppression,  it  has  a  loss  of  about  -1.5  dB,  which  is  the 
same  level  of  performance  of  cancellation  as  in  the  ideal  case  of  Table  3.5.3-2. 

I/Q  Phase  Imbalance  also  affects  the  Q  component  of  the  data,  although 
in  a  different  manner.  It  changes  slightly  the  phase  of  the  Q  component.  Figure 
3.5.3-4  shows  the  effect  of  this  error  on  the  cancellation  performance. 
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Rgure  3.5.3-4.  Cancellation  loss  vs.  I/Q  Phase  Imbalance 

l/Q  Phase  Imbalance  has  about  the  same  effect  upon  the  performance  as 
the  l/Q  Amplitude  Imbalance  error.  Range  response  cancellation  has  a  slight 
gain  over  the  ideal  case  while  range  suppression  is  about  -1 .5  dB,  which  is  as 
good  as  the  ideal  case  of  Table  3.5.3-2. 

Power  Supply  Amplitude  Modulation  error  is  examined  next.  This  error  is 
due  to  the  system  power  supply  modulating  the  target  signal  with  a  small  60  Hz 
sinusoidal  ripple.  The  modulation  is  given  in  terms  of  percent  modulation.  The 
cancellation  loss  performance  is  given  in  Figure  3.5.3-5. 
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Figure  3.5.3-5. 


Cancellation  loss  vs.  Power  Supply  Amplitude  Modulation 


Power  Supply  Amplitude  Modulation  does  not  seem  to  affect  the  system. 
The  losses  are  the  same  level  as  the  ideal  case. 


The  next  error  is  timing  jitter.  This  is  due  to  the  system  clock  not  being 
precise  at  its  transitions.  Figure  3.5.3-6  shows  the  effect  of  this  error  on 
cancellation.  As  is  apparent  there,  the  losses  are  independent  of  timing  jitter. 
The  losses  are  the  same  level  as  the  ideal  case. 
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Timing  Jitter  (%) 


Figure  3.5.3-6.  Cancellation  loss  vs.  Timing  Jitter 

During  each  frame  the  radar  receives  many  pulses  back  from  the  target. 
The  range  to  the  target  is  changing  continuously  during  the  frame  time  due  to 
the  velocity  of  the  target.  Hence,  the  target  will  be  in  a  slightly  different  range 
bin  position  for  each  received  pulse.  This  is  called  range  walk  error.  Figure 
3.5.3-7  shows  the  result  of  this  error  on  the  cancellation  loss. 
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Figure  3.5.3-7.  Cancellation  loss  vs.  Range  Walk  Error 

There  is  no  discernible  pattern  to  the  lossed  created  by  range  walk  error. 
For  the  range  response  cancellation  and  zero-out  blanking,  the  velocity  of  Mach 
1  (300  m/sec)  seems  to  create  the  least  loss.  Incidentally,  it  must  be  noted  that 
in  the  ideal  case  of  Table  3.5.3-2.  range  walk  was  present,  at  the  velocity  of  300 
m/sec. 


The  final  error  to  be  looked  at  are  the  phase  transition  errors.  These  are 
the  errors  in  the  phase  transition  from  one  Barker  bit  to  the  next.  Two  possibile 
errors  sequences  of  length  13  were  looked  at  and  summarized  in  Table  3.5.3-3. 
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TABLE  3.5.3-3.  LOSS  WITH  2  PHASE  TRANSITION  ERROR  SEQUENCES 


PHASE  ERRORS 

SIGNAL  LOSS  (dB) 

RANGE  RESPONSE 
CANCELLATION 

ZERO-OUT 

BLANKING 

RANGE 

SUPRESSION 

r,r,3*,4«,3M3M4, 

4®,5“,-5”,5',0“,-10® 

-0.1 

-9.07 

-1.49 

5',6'’,5“,6“,7“,ir,18*, 

8*,9«,-1*,9»,19*,9* 

-0.1 

-9.07 

-1.49 

The  losses  are  exactly  the  same  as  in  the  ideal  case  of  Table  3.5.3-2.  Hence, 
phase  errors  are  not  important. 

Finally  in  Table  3.5.3-4  we  list  all  the  error  sources.  They  are  ranked  in 
terms  of  the  loss  of  the  range  response  cancellation  technique.  Several  things 
can  be  concluded  from  this  table.  First,  video  bandwidth  is  the  most  severe 
source  of  error.  A  large  video  bandwidth  of  at  least  6  is  crucial  to  good 
cancellation  performance.  Other  error  sources  are  negligible  by  comparison. 
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TABLE  3.5.3-4.  ERROR  SOURCES  RANKED  WITH  LOSS 


Error  Source  /Loss  In  dB 


Video  Bandwidth  (Alpha=3.) 


Video  Bandwidth  (Alpha=4.) 


Range  Walk  (125  m/sec) 


Video  Bandwidth  (AlphaaS.) 


Range  Walk  (758  m/sec) 


Video  Bandwidth  (Alphaa6.) 


Quantization  (10  bits) 


Quantization'  (12  bits) 


Quantization  (14  bits) 


Timing  Jitter  (0.001%) _ 


Timing  Jitter  (0.005%) 


Timing  Jitter  (0.01%) 


Range  Walk  (300  m/sec) 


Phase  Errors  (a) 


Phase  Errors  (b) 


Power  supply  Amp.  Mod.  (0.1  dB,  60  Hz) 


Power  supply  Amp.  Mod.  (0.2  dB,  60  Hz) 


Range 

can. 


Blanking 


-8.76 


-8.04 


-7.66 


-7.71 


-10.25 


-7.01 


-7.01 


-7.01 


-7.01 


-7.01 


-7.01 


-7.01 


-7.01 


-7.01 


-7.01 


-7.01 


-7.01 


Power  supply  Amp.  Mod.  (0.4  dB,  60  Hz) 

1.97 

-7.01 

0.58 

l/Q  Phase  Imbalance  (0.1  dB,  30®) 

2.01 

-6,98 

0.67 

l/Q  Amplitude  Imbalance  (0.1  dB.  10®) 

2.17 

-6.98 

0.67 

l/Q  Phase  Imbalance  (0.1  dB,  10®) 

2.17 

-6.98 

0.67 

l/Q  Phase  Imbalance  (0.1  dB.  20®) 

2.18 

-6.98 

0.68 

l/Q  Amplitude  Imbalance  (0.2  dB,  10®) 

2.24 

-6.95 

0.75 

l/Q  Amplitude  Imbalance  (0.4  dB,  10®) 

2.39 

-6.9 

0.92 

(a)  -  1®.2®.3®,4®.3M3®.14®.4®,5®. 
-5®,5°.0°,-10® 

(b)  =  5^6®.5®,6®,7®,17®,18®,8®,9®, 
-r,9°,19°,9° 

Assume  20  dB  dyn.  range  except  where 
noted. 

Small  target  is  4  RB  away. 

3.5.4  Processing  Raouirements 


Adaptive  processing  techniques  can  have  moderate  to  extreme  resource 
demands  depending  on  the  available  number  of  degrees  of  freedom,  N^.  For 

example,  ECM  noise  cancellation  throughput  requirements  are  proportional  to 

Ng.  Batch  processing  demands  memory  capacity  proportional  to  the  number  of 
channels.  It  is  recognized  that  with  judicious  application  of  large  target 
cancellation  algorithm,  one  will  have  a  minimal  demand  of  processing 
resources. 

Three  criteria  can  be  used  to  evaluate  the  processing  resource  demands 
for  each  candidate  technique.  Loading  or  throughput,  which  is  the  number  of 
computer  arithmetic  operations  required  by  each  algorithm,  is  the  most 
important  criterion.  Second  is  the  VAX  CPU  time  required  and  third,  as  a 
byproduct,  the  number  of  VAX  Fortran  77  instructions  required  for  executing  the 
algorithms. 

Corresponding  to  each  technique.  Figure  3.5.4-1  parametrically  provides 
the  computational  loading  formulas  and  Figure  3.5. 4-2  shows  an  example  of 
throughput  loading  for  a  set  of  waveform  parameters.  The  notation  for  the 
formulas  is  given  in  Table  3.5.4-1.  Finally,  in  Figure  3. 5.4-3,  the  Fortran  lines  of 
code  and  VAX  CPU  Secs,  required  for  each  algorithm  is  tabulated. 
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NO  CANCELLATION 

NFFrNRB*(5*log2(NFFT)  +  2*NPCR  +  NTH  +  6) 

CANCELLATION  with  PC  Response  " 

NFFT*NRB*(5*log2(NFFT)  +  2*NPCR  +  NTH  +  6) 

+  NFFT*NRB*(NTH+1) 

+  6*NTAR*(NPCR+2) 

BLANKING  (Zero  Out) 

NFFT*NRB*{5*log2(NFFT)  +  2*NPCR  +  NTH  +  6) 

+  NFFT*NRB*(2*NPCR  +  NTH+4) 

+  3*NTAR*(NPCR+1) 

RANGE  SIDELOBE  SUPPRESSION 

NFFT*NRB*{5*log2(NFFT)  +  2“NPCR  +  NTH  +  6) 

+  NFFT*NRB*(NTH+4)  +  NV*NFFr  (4*Nt-2) 

+  NV*NV*(2*Nt-1/2)  +  Nt'NV 

DOPPLER  CANCELLATION  IN  TIME  DOMAIN 

NFFrNRB*(5*log2{NFFT)  +  2*NPCR  +  NTH  +  6) 

+  NITER*tNPCR*NFFr(5*log2(NFFT) 

+  2*NPCR  +  NTH  +  6)  +  8*NFFT  ] 

Assume  complex  operation  =  10  real  operations 

complex  add  =  2  real  operations 

real  operation  real  add  or  real  multiply  or  clear 


Figure  3.5.4-1.  Signal  Processing  Loading  Computation  Formulas 
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#  of  Real  Operations 
per  Array  Time 

No  Cancellation 

540544 

103.00 

Range 

Response 

Cancellation 

602382 

114.78 

Blanking 

650494 

123.95 

Range 

Sidelobe 

Suppression 

718090 

136.83 

Doppler  Cancel, 
in  the 

Time  Domain 

885376 

168.71 

PRI 

= 

41  ms 

NFFT 

s 

128 

NRB 

3 

41 

NPCR 

3 

13 

NTAR 

3 

1 

NTH 

3 

24 

NK 

3 

15 

Array  Time 

AT 

-  5.248  ms 

Figure  3.5,4-2.  An  Example  of  Throughput  Calculations 


Lines  of 

FORTRAN  Code 

CPU  Secs. 

Range 

Response 

Cancellation 

145 

5.78 

Blanking 

122 

6.42 

Range 

Sidelobe 

Suppression 

172 

6.41 

Doppler  Cancel, 
in  the 

Time  Domain 

N/A 

N/A 

Figure  3.5.4-3.  Lines  of  code  and  CPU  Secs,  for  the  Algorithms 
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TABLE  3.5.4-1 .  NOTATION  FOR  PROCESSOR  LOADING  COMPUTATIONS 


NOTATION 

NFFT 

Number  of  doppler  filters 

NRB 

Number  of  range  bins 

NPCR 

Pulse  Compression  ratio 

NTH 

Number  of  cells  used  in  threshol  estimate 

PRI 

Pulse  repetition  interval 

NTAR 

Number  of  large  targets  detected 

NDET 

Number  of  (large  target)  detections 

(=3NTAR) 

NV 

The  number  of  sidelobes  suppressed 

Nt 

#of  digits  in  the  range  suppression  code 

NITER 

Number  of  iterations  desired  for  the 

doppler  can.  in  time  domain  technique 

CO 

complex  operation  (*  &  +  =10  ro) 

ca 

complex  add  (=co/2  other  use  co/5) 

ra 

read  add 

rm 

real  multiplier  (=co/4) 

The  results  of  these  Figures  can  be  summarized.  In  the  range  domain, 
the  range  response  cancellation  technique  requires  the  least  processing 
resources  overall  while  the  range  suppression  requires  the  most.  This  only 
reaffirms  the  conclusion  partially  reached  in  the  previous  Sections  that  the 
range  response  cancellation  technique  is  the  most  valuable,  it  is  effective  in 
cancellation,  not  subjected  to  as  much  errors,  and  also  the  simplest  to 
implement.  Hence,  this  technique  should  definitely  be  the  first  technique  to  be 
used  for  co-target  interference  rejection.  It  must  be  kept  in  mind,  however,  that 
the  range  suppression  technique  is  only  20%  more  resource  intensive  than 
range  response,  as  can  be  seen  in  Figure  3.5.4-2,  while  its  effectiveness  in 
revealing  the  small  target  is  almost  as  good.  Zero-out  blanking  requires  only 
8%  more  resources  than  range  response  cancellation,  but  its  performance,  as 
seen  in  the  previous  Sections  is  limited. 
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The  doppler  cancellation  in  time  domain  technique  is  the  most 
processing  intensive.  However,  it  is  a  very  effective  technique  and  should  be 
implemented  even  with  its  high  cost. 


3.5.5  Technique  Summary 

Table  3.5.5-1  summarizes  the  effectiveness  of  the  three  cancellation 
algortihms  that  are  effective  in  rejecting  co-target  interference.  Notice  that  the 
zero-out  blanking  technique  has  been  altogether  abandoned  since  it  requires  at 
least  5  range  bin  target  separation  to  be  effective.  The  two  remaining 
techniques  in  the  range  domain  show  very  similar  performance  except  for  the 
throughput  requirements,  where  range  suppression  requires  about  20%  more 
throughput  than  range  response  cancellation.  Airborne  radars  are  limited  in 
their  available  processing  resources,  and  this  could  be  the  deciding  factor  in 
favoring  range  response  cancellation  over  the  range  suppression  technique. 
Nevertheless,  it  is  important  to  have  two  available  techniques  in  the  range 
domain  to  expand  the  available  options.  In  the  doppler  domain,  even  though  it 
is  processing  intensive,  we  have  available  a  technique  that  is  effective  in 
discriminating  targets  as  close  as  1  doppler  filter  apart. 


TABLE  3.5.5-1.  EFFECTIVENESS  SUMMARY  FOR  ALGORITHMS 


TECHNIQUE 

REGION  OF 
EFFECTIVENESS 

(la  fcparatlon) 

LOSS 

(<») 

THROUGHPUT 

(Mrop) 

Range  domain 

RANGE  RESPONSE  CANCEL 

@  2  RB  and  up 

0.099 

114.78 

RANGE  SUPPRESSION 

@  2  RB  and  up 

1.49 

136.83 

PooBtef-domain 

DOPPLER  CANCEL 

@  1  filter  and  up 

0.60 

168.71 

IN  TIME  DOMAIN 
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3.6  Tracking  ImolementatiQn 


Given  that  co-target  interference  reduction  can  be  achieved,  the  system 
level  effectiveness  of  the  reduction  depends  on  the  use  of  the  resulting  track 
data.  For  example,  if  a  large  target  and  a  small  target  are  flying  in  formation  to 
deliberately  conceal  the  existence  of  the  small  target,  it  would  be  sufficient  to 
establish  the  presence  of  the  small  target  and  check  from  time  to  time  to  insure 
that  the  formation  remained  unchanged.  In  handing  off  to  local  fire  control 
assets,  the  message  would  be  that  one  (or  more)  small  targets  are  in  close 
proximity  to  the  designated  bright  target,  and  the  local  assets  could  then  acquire 
the  small  target  on  their  own.  This  would  essentially  be  a  sun/eillance  raid 
assessment  strategy.  On  the  other  hand,  the  local  assets  would  need  the  CTIR 
techniques  to  successfully  break  out  and  engage  the  small  target(s). 

On  the  other  hand,  if  target  (between  large  and  small)  separation  was 
sufficiently  large,  keeping  individual  tracks  of  all  targets  (which  would  involve 
repeatedly  applying  CTIR  techniques)  could  be  of  tactical  value  if  local  fire 
control  assets  could  benefit  from  such  information.  Thus,  the  system 
effectiveness  of  a  particular  tracking  strategy  using  CTIR  techniques  depends 
on  the  value  of  the  data  in  the  fire  control/engagement  (or  situation  assessment) 
portion  of  the  mission. 

In  this  Section  we  examine  the  value  of  using  CTIR  techniques  and  the 
resulting  track  data  in  a  comprehensive  Multiple  Target  Tracking  (MTT) 
simulation  developed  at  Hughes.  A  typical  MTT  system  consists  of  five  basic 
elements,  illustrated  in  Figure  3.6-1: 

•  Sensor  data  processing  and  measurement  formation, 

•  Data  correlation, 

•  Track  maintenance  logic  (initiation,  confirmation,  deletion), 

•  Estimate  filtering  and  prediction,  and 

•  Gating  computations. 
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Figure  3.6-1 .  ■  Elements  of  a  Typical  MTT  System. 

Sensor  data  processing  and  measurement  formation  involves  signal 
processing,  including  clutter  and  interference  rejection  processing,  range  bin 
and  doppler  filter  formation,  and  threshold  selection  and  testing  for  detection. 
Measurements  typically  consist  of  range,  doppler,  and  angle  discriminants. 
However,  when  large  dynamic  ranges  exist  between  closely  spaced  targets, 
there  are  likely  to  be  characteristics  of  the  measurements  that  differ  from  those 
obtained  from  "standard"  targets.  These  characteristics  are  called 
measurement  biases  and  must  be  understood  and  taken  into  account  from  the 
start,  since  they  may  have  a  significant  impact  on  tracking  system  design  and 
performance.  The  biases  are  examined  in  greater  details  later  on  in  this 
section. 

Data  correlation  involves  assigning  new  measurement  data  to  existing 
tracks.  This  process  can  be  very  difficult  in  closely  spaced  target  situations.  If 
targets  are  spaced  widely  enough  relative  to  track  and  measurement  accuracy, 
unambiguous  correlation  is  easy,  and  almost  any  technique  will  work 
(sequential  Nearest  Neighbor,  for  example).  If  targets  are  very  closely  spaced 
(relative  to  track  and  measurement  accuracy),  then  miscorrelation  will  occur 
frequently,  but  resulting  track(s)  will  be  stable.  In  this  case,  the  targets  are  too 
close  to  be  tracked  separately,  and  group  or  formation  tracking  is  required. 
Group  tracks  consist  of  a  group  centroid  estimate,  a  group  count,  and  potentially 
a  group  extent  estimate— the  centroid  estimate  is  based  on  measurements  of 
several  targets,  and  it  is  the  centroid  estimate  that  is  used  for  track  filtering  and 
prediction.  Group  tracking  is  appropriate  when  measurements  from  individual 
elements  of  the  different  targets  cannot  be  reliably  differentiated  or  resolved. 
Formation  tracking  differs  from  group  tracking  in  that  although  measurements  of 
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individual  targets  can  be  resolved,  the  targets  are  too  closely  spaced  to  allow 
maintenance  of  unambiguous  tracks.  A  formation  track  consists  of  several 
target  tracks,  each  propagated  by  an  average  formation  velocity.  Since  each 
element  is  tracked  individually,  no  estimate  of  group  count  or  extent  is  required. 
Errors  in  correlation  (assigning  the  wrong  measurement  to  a  track)  are 
expected,  and  occur  without  degrading  tracking  performance. 

When  targets  are  too  closely  spaced  to  be  tracked  unambiguously,  but 
are  too  widely  spaced  to  be  tracked  as  a  group  or  formation,  the  tracking 
process  Can  become  unstable.  Errors  in  correlation  can  lead  to  erroneous 
velocity  estimates,  which  in  turn  can  result  in  track  divergence,  track  loss,  and 
false  tracks.  This  type  of  behavior  can  occur  when  targets  are  close  enough 
that  the  probability  of  false  correlation  is  significant,  but  are  spaced  widely 
enough  that  correlation  errors  can  lead  to  substantial  kinematic  errors.  In  this 
environment,  special  approaches  are  required  to  minimize  the  probability  of 
false  correlation  (e.g.  Multiple  Hypothesis  Tracking)  and/or  minimize  the 
kinematic  errors  resulting  from  false  correlation  (such  as  Joint  Probabilistic  Data 
Association). 

Track  initiation,  confirmation,  and  deletion  are  straightforward  in 
unambiguous  tracking  regimes,  but  difficult  in  unstable  situations.  Track 
initiation  is  a  tentative  decision  that  certain  measurements  are  from  a  target  not 
already  being  tracked.  Track  confirmation  results  when  this  decision  is  made 
with  more  confidence.  Deletion  results  when  it  is  decided  that  a  track  is  no 
longer  being  updated  with  valid  target  measurements,  either  because  a  target  is 
no  longer  resolvable  (or  has  disappeared  for  any  other  reason),  or  because  the 
track  has  diverged  and  the  true  target  is  no  longer  within  the  track  gates. 

Filtering  and  prediction  of  estimates  of  target  kinematics  are  performed 
by  “smoothing"  in  some  sense  a  sequence  of  noisy  measurements.  When 
detection  probability  is  high,  and  targets  are  widely  spaced,  simple  filtering 
techniques  (such  as  fixed  coefficient  Alpha-Beta  or  Alpha-Beta-Gamma  filters) 
are  sufficient.  When  detection  probability  drops,  or  targets  are  more  closely 
spaced,  time-varying  approaches  such  as  Kalman  filters  become  more 
desirable  to  extract  as  much  information  as  possible  from  the  available  data. 
Although  it  is  possible  to  use  look-up  tables  for  fixed  coefficient  filters  to  adapt 
the  filter  gains  to  changes  in  target  behavior,  the  resulting  filters  are 
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complicated,  and  a  simple  (low  order)  Kalman  filter  is  straightforward  to  design 
and  implement.  Unless  large  numbers  of  targets  are  expected  or  very  high 
data-rates  are  involved,  the  processing  load  due  to  track  kinematic  filtering  is 
generally  much  lower  than  the  loading  due  to  signal  processing,  since  track 
computations  are  required  relatively  infrequently  (at  intervals  on  the  order  of  .5 
to  5  seconds).  The  LDR  Simulation  uses  a  9-state  Kalman  filter. 

In  order  to  decide  which  measurements  to  use  to  update  each  track,  a 
measure  of  the  distance  between  the  predicted  track  and  each  observation  is 
computed.  A  tracking  gate  is  a  threshold  applied  to  this  distance:  if  the  distance 
is  exceeded,  the  measurement  is  not  used  for  updating  a  track  (the 
measurement  is  not  within  the  track's  gates),  while  if  the  distance  is  less  than 
the  threshold  (the  measurement  is  within  the  track’s  gates),  the  measurement  is 
considered  for  use  in  updating  the  track  (but  not  necessarily  used).  Gating  is  a 
way  of  simplifying  the  correlation  process  that  immediately  rules  out  unlikely 
correlation  possibilities.  To  form  accurate  gates,  an  estimate  of  target  track  error 
is  required — the  larger  the  estimated  uncertainty,  the  larger  the  gates,  and 
conversely.  The  use  of  Kalman  filters  is  very  desirable  from  this  perspective, 
since  an  error  variance  is  explicitly  propagated  with  the  track  state  estimates 
that  takes  into  account  track  update  rate,  measurement  accuracy,  and 
anticipated  target  maneuver  capability. 

Designing  a  tracking  system  involves  selecting  one  or  more  techniques 
to  perform  each  of  these  functions.  It  is  desirable  to  avoid  unnecessary 
complexity  and  computational  load,  but  still  provide  acceptable  performance. 
Target  spacing,  detection  probability,  maneuver  capabilities,  and  the  effects  on 
measurements  resulting  from  the  specific  signal  processing  used  in  this  system 
must  be  taken  into  account,  and  weighed  against  available  processing 
resources,  number  of  targets  expected,  and  anticipated  antenna  loading. 

The  aforementioned  MTT  Simulation  available  at  Hughes  contains  all 
these  elements.  Specifically,  this  Simulation  uses  Multiple  Hypothesis  Tracking 
(MHT)  in  doing  data  correlation  to  existing  tracks  and  to  generate  track 
hypotheses.  The  LDR  Program  utilizes  the  information  about  co-target 
interference  and  combine  this  with  the  MHT  Simulation.  The  value  of  doing  this 
is  examined  in  a  later  section. 
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3.6.1  Filtering 


Of  the  5  basic  elements  of  an  MTT  system  described  in  Section  3.6,  only 
the  filtering  and  prediction  element  is  described  in  detail.  This  is  because  it  was 
necessary  to  modify  the  existing  MHT  Simulation  to  incorporate  a  filter  that  can 
accomodate  a  true  radar  measurement  model.  Filtering  and  prediction  are  the 
fundamental  elements  of  any  tracking  system.  They  are  used  to  estimate 
present  and  future  target  kinematic  quantities  such  as  position,  velocity,  and 
acceleration.  In  a  tracking  simulation  when  the  detection  probability  is  high, 
and  targets  are  widely  spaced,  simple  filtering  techniques  (such  as  fixed 
coefficient  a-p  filters)  are  sufficient.  When  detection  probability  drops,  or  targets 
are  more  closely  spaced,  time-varying  approaches  such  as  Kalman  filters 
become  more  desirable  to  extract  as  much  information  as  possible  from  the 
available  data.  This  latter  case  fits  more  closely  with  scenarios  where  the  LDR 
techniques  would  be  applied  to. 

The  choice  of  a  coordinate  system  is  an  important  element  of  a  tracking 
system.  In  most  MTT  systems,  it  is  important  to  express  tracks  in  a  common 
coordinate  system  to  facilitate  data  association  and  tactical  situation 
assessment.  A  pseudo-inertial  coordinate  system  such  as  NED  (North,  East, 
Down)  is  particularly  useful  for  airborne  systems,  but  it  is  also  applicable  for 
surface  tracking  systems.  This  coordinate  system  is  shown  in  Figure  3.6. 1-1. 
The  radar  system  is  at  the  origin. 
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Figure  3,6.1 -1.  Radar  centered  inertial  (NED)  coordinate  system 

The  north,  east,  and  down  coordinates  are  referred  to  from  now  on  as  the 
X,  y,  and  z  coordinates,  respectively.  The  Kalman  filter  implemented  in  the  LDR 
Tracking  Simulation  is  a  9-stato  filter,  tracking  the  position,  velocity,  and 
acceleration  states  for  each  of  the  three  coordintates.  Thus  the  state  vector 
containing  the  tracked  target  dynamics  can  be  written  as: 

■ 

xi(k) 

X2(k) 

X3(k) 

X4(k) 

x(k)  =  xsOc) 

X6(k) 

X7(k) 

X8(k) 

_  X9(k) 

where  the  dot  means  velocity  and  double  dot  means  acceleration. 
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The  target  dynamic  process  is  assumed  to  be  modeled  in  the  discrete 
Markov  form: 

X  (k+1 )  =  4>  X  (k)  +  q  (k)  +  f  (k+1 1  k), 

where  O  is  the  transittion  matrix;  q  (k)  is  the  zero-mean,  white,  Gaussian  noise 
process  with  known  covariance  Q;  and  f  (k+1|  k)  is  a  known  deterministic  input, 
such  as  the  relative  position  change  associated  with  own-ship  motion.  In  the 
LDR  Interaction  Simulation  f(k+1|  k)  is  zero  because  stationary  radars  are 
assumed.  The  discrete-time  Markov  process  can  be  defined  as  a  process  in 
which  the  statistical  representation  of  the  process  in  the  future  (scan  k+1)  is 
completely  determined  by  the  present  state  (scan  k). 

The  Singer  target  maneuver  model  is  used.  In  this  model,  the 
acceleration  is  modeled  as  a  first-order  Markov  porcess,  which  can  be  written 

a  (k+1 )  =  p^a  (k)  +  Vi  -p^  r  (k) 

The  maneuver  correlation  coefficient  is  defined  in  terms  of  the  maneuver 
time  constant  tm 

jr 

Pm=e 

and  Om  is  the  maneuver  standard  deviation.  r(k)  is  the  standard  normal  random 
variable  while  T  is  the  sampling  interval. 

With  the  Singer  model,  and  defining  3  to  be  l/x^.  the  transition  matrix  is 
given  by: 
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The  matrix  shown  is  for  the  3-state  case.  For  our  9-state  filter,  the  matrix  is 
duplicated  down  the  main  diagonal: 


■ 

* 

^3x3 

^3x3 

0  = 

^3x3 

'^>3x3 

^3x3 

^3x3 

^3x3 

^3x3 

where  03x3  is  the  3  x  3  zero  matrix. 

The  process  (or  maneuver)  excitation  covariance  matrix  is: 


Q  = 


Q 

0 

0 


*3x3 

^3x3 

^3x3 

3x3 

^3x3 

^3x3 

3x3 

^3x3 

^3x3 

where  the  matrix  Q3X3  is  given  in  Reference  3,  pg.  32. 

The  radar  measurements  give  the  target  range  (r),  range  rate  (r),  and 
angles  ("n  and  e).  These  measured  quantities  are  functions  of  the  state  variables 
in  the  state  vector,  corrupted  by  uncorrelated  noise.  Thus,  the  4-dimensional 
measurement  vector  is  modeled  as 

y  (k)  =  h  (x(k))  -h  v(k),  (Eqn.  3.6.1-1) 

where 


h  (x(k))  = 


r  (x(k)) 
r  (x(k)) 
(x(k)) 
£  (x(k))_ 


^/"l  2  2 

V  X  +  y  +  z 

XX  H-  yy  -I-  zz 

J  2  2  2 

V  X  -h  y  +  z 

tan*^  (y/x) 

-1,  ./""2  2, 

tan  (-Z/V  X  +  y  ) 
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v(k)  is  zero  mean,  white,  Gaussian  measurement  noise  with  covariance 
matrix  R  given  by 


The  measurement  accuracy  of  the  radar  is  specified  by  the  terms  in  the  R  matrix. 
For  example,  range  measurement  accuracy  is  given  by  Of. 

The  generic  Kalman  filter  equations  are  given  by; 

X  (k|  k)  =  X  (k|  k-1 )  +  K  (k)  [  y  (k)  -  H  x  (k|  k-1 )  ]  (Eqn.  3.6.1  -2) 

K(k)  =  P(k|  k-1)H^[H  P(k|k-1)H^+R]'’ 

P(k|  k)=  [l-K(k)H]P(k|k-1) 

X  (k+1 1  k)  =  <D  X  (k|  k) 

P(k+1|k)  =  <DP(k|  k)  (D^+Q, 

where  x  (k+1|  k)  is  the  1-step  predicted  estimate  of  the  state  vector  at  time  k, 

and  X  (kl  k)  is  the  smoothed  estimate  of  the  state  vector  at  time  k,  K  is  the  gain 
matrix,  H  is  the  measurement  matrix,  and  P  is  the  covariance  matrix  of  the 
estimation  error; 

P(k)  .  E{[x(k)-J(k)l[x(k)-i(k)l\ 

Note  that  in  the  generic  equation  3.6.1 -2  above,  the  measurement  matrix 
H  is  a  linear  matrix.  However,  our  measurement  model  in  Eqn.  3.6. 1-1  uses  a 
nonlinear  h  function.  In  these  cases,  it  is  necessary  to  linearize  the  h  function 
into  an  H  matrix  through  a  Taylor  expansion  method,  which  results  in  what  is 
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called  the  extended  Kalman  filteri3.  The  H  matrix  for  our  filter  is  a  4  x  9  matrix 
and  is  given  by: 


H(x(k|k-1))  = 


2  2 

% 


8h(x) 

dx 


X  =  X  (k|  k-l ) 


j  r  ,  . .  _  . . . 


where  r=vx^  +  y^+z^_ 

X,  y,  z,  X,  y,  z  are  formed  from  the  predicted  state  vector  x  (k|  k*1 ) . 


Also,  the  state  estimate  update  equation,  Eqn.  3.6. 1-2,  needs  to  be 
modified  for  the  extended  Kalman  filter.  It  is  given  by: 

X  (k|  k)  =  i  (k|  k-1)  +  K  (k)  [y  (k)  -  h(  x  (k|  k-1))  ]. 


3.6.2  Masking  Region  and  Measurement  Biases 

A  large  target  will  mask  a  smaller  target  from  detection  if  their  separation 
is  very  small.  Depending  upon  the  geometry,  the  masking  region  may  extend 
for  a  significant  amount  of  time.  Figure  3.6.2-1  (a),  (b),  and  (c)  show  the 
minimum  separation  required  in  range,  doppler,  and  angle  for  small  target 
detection  when  no  LDR  cancellation  techniques  is  used.  In  each  plot,  the 
power  ratio  refers  to  the  dynamic  range  between  large  and  small  target.  Using 
20  dB  dynamic  range  as  a  reasonable  level,  we  conclude  from  Figure  3.6. 1-2 
(a)  that  using  the  Barker  13:1  code,  a  separation  of  10  range  bins  in  range 
would  be  necessary  to  detect  the  small  target.  Of  course,  with  even  larger 
dynamic  range,  14  range  bin  separation  would  be  required  to  detect  the  small 
target.  For  the  cases  we  present,  however,  a  20  dB  dynamic  range  is  in  effect 
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SMALL  SIGNAL  -  20  DB 
I  3: 1  BARKER  CODE 


Figure  3.6.2-1  (a).  Minimum  range  separation  for  small  target  detection 

and  so  10  range  bins  is  used.  Similarly,  in  the  doppler  domain,  assuming  75 
dB  Dolph-Chebyshev  filter  amplitude  weighting,  we  see  that  a  2.75  doppler  filter 
separation  is  necessary,  and  in  angle  a  0.8  beamwidth  separation  is  necessary 
for  small  target  detection. 
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SMALL  SIGNAL  20  DB 

75  DB  DOLPH  CHEBYSHEV  WEIGHTING 


Figure  3.6.2-1  (b).  Minimum  filter  separation  for  small  target  detection 


Figure  3.6.2-1  (c).  Minimum  angular  separation  for  small  target  detection 
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Based  on  these  results  we  can  define  a  masking  region  as  a  region 
where  only  the  large  target  can  be  observed.  This  region  is  illustrated  in  Figure 
3.6.2-2. 


Target  1  Track 


Masking  Region 


Figure  3.6.2-2.  Masking  region  geometry 


The  conditions  for  the  two  targets  to  be  in  the  masking  region  are; 

Rdiff<  10  Range  bins 
and 

Rdiff  <  3  Doppler  filters 
and 

A2<jiff  <  0.8  Beamwidth 
and 

Eljjjff  <  0.8  Beamwidth 


With  or  without  using  LDR  techniques  to  remove  large  target  interference, 
whenever  the  two  targets  are  in  the  masking  region,  the  measurements  will  be 
modified  by  a  bias  term.  The  bias  is  due  to  interference  by  each  target  on  one 
another,  arxJ  is  most  severe  when  the  dynamic  range  is  small.  Figures  3.6.2-3 
(a),  (b),  and  (c)  show  the  large  target  measurement  biases  due  to  small  target 
interference  in  the  range,  doppler,  and  angle  domains. 


168 


SCf’ARATION  (RBI 


13 


Figure  3.6.2“3  (a).  Measurement  bias  in  range 
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Figure  3.6.2-3  (b).  Measurement  bias  in  doppler 
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Figure  3.6. 2-3  (c).  Measurement  bias  in  angle 


In  all  three  Figures,  the  bias  is  plotted  with  the  dynamic  range  in  dB  as 
the  parameter.  The  bias  is  worse  when  the  dynamic  range  is  small.  Some  of 
the  biases  can  be  very  substantial,  depending  upon  the  separation  and 
dynamic  range.  In  doppler,  the  bias  is  as  large  as  0.8  of  a  filter  when  the 
dynamic  range  is  5  dB  and  target  separation  is  2  filters.  In  range,  the  bias  is  as 
large  as  0.3  of  a  range  bin  when  the  dynamic  range  is  5  dB  and  separation  of 
about  1  range  bin.  The  maximum  bias  for  angle  is  about  0.14  beamwidth. 

Now  that  we  have  specified  the  measurement  biases,  we  finally  come  to 
the  model  of  target  track  observations  that  will  be  fed  into  the  MHT  Simulation  to 
evaluate  track  effectiveness.  There  are  two  scenarios,  with  and  without  LDR 
techniques,  and  these  are  summarized  in  Figure  3.6.2-4. 
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1.  Without  LDR  Techniques 
r^ase  a.  Targets  in  ROM 

Measurement  is  from  large  target  only  -  with  bias 

Case  b.  Targets  outside  ROM 

Measurements  are  from  both  targets  -  no  biases 

fL=*L 

2.  With  LDR  Techniques 
Case  a.  Targets  in  ROM 

Measurements  are  from  both  targets  •  with  biases 

fs-is  +-  +Ss  (j!s*l!LOfSNR«,=  0(B) 

Case  b.  Targets  outside  ROM 

Measurements  are  from  both  targets  -  no  biases 

2s  =  -s  +- 

t 

where  v  is  gaussian  noise  with  mean  0  and  covariance  R. 


Rgure  3.6.2-4.  Track  measurement  model  with  and  without  LDR  techniques 

The  covariance  matrix  R  is  the  same  measurement  variance  matrix  R  defined  in 
Section  3.6.1  for  the  Kalman  filter. 
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3.6.3  Results  of  Tracking 


We  present  color  RAMTEK  charts  that  show  the  target  tracks  as  well  as 
the  likelihood  probability  of  the  tracks.  In  all  these  charts,  the  X  and  Y  scale  are 
labeled  in  pixel  units,  where  each  pixel  is  1/10  nautical  miles.  The  radar  is 
always  at  the  origin  (0,0).  The  velocities  of  the  targets  are  about  Mach  1 . 

The  tracks  are  color-coded  in  terms  of  probability  of  tracks.  Red  is  the 
highest  probability  track  while  blue  and  violet  are  the  lowest  probability  tracks. 
The  yellow  crosses  in  the  charts  are  some  of  the  observations  that  were  not 
erased  from  the  screen.  Faint  violet  circles  are  the  covariance  gates  that  the 
MHT  algorithm  puts  around  an  observation. 

In  Figure  3.6.3-1,  we  have  the  launching  of  a  missile  when  LDR 
techniques  are  not  used.  In  this  chart  we  simulate  the  launching  of  a  missile 
from  an  original  target.  It  can  be  a  bomber  launching  a  cruise  missile.  The 
bomber  is  easily  tracked  from  about  (X=460,Y=180)  all  the  way  to  the  end  of 
the  run  (-25,180).  This  track  is  in  red.  However,  at  some  point  in  time  a  missile 
was  launched,  as  can  be  seen  by  the  second  red  track  above  the  bomber  track. 
We  cannot  detect  the  origin  of  this  launch  since  the  missile  is  masked  by  the 
bomber.  All  we  see  is  the  sudden  appearance  of  a  second  red  track  at 
(160,220).  But  by  this  time  many  seconds  have  gone  by  since  the  launch  of  the 
missile. 

In  Figure  3.6. 3-2,  we  have  the  launching  of  a  missile  when  LDR 
techniques  are  used.  In  this  chart  we  have  the  same  scenario  as  in  the  previous 
Figure  except  that  here  LDR  techniques  are  used.  The  bomber  is  easily  tracked 
from  about  (X=460,Y=180)  to  the  middle  of  the  chart  (225,180).  This  track  is  in 
red.  Then  it  launches  a  missile  flying  upward.  Initially  the  missile  is  very  close 
to  the  bomber  and  thus  Is  masked.  However,  by  using  LDR  techniques  it  is 
possible  to  obtain  observations  from  the  missile  within  the  masking  region. 
Because  of  this  the  MHT  algorithm  is  able  to  hypothesize  that  there  are  two  low 
probability  tracks  close  to  each  other  (green  and  blue).  It  continues  to  track  the 
two  targets  until  they  separate  far  enough  to  be  tracked  as  red  tracks.  These 
results  suggest  that  by  using  LDR  techniques  it  is  possible  to  detect  a  cruise 
missile  launch. 
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In  Figure  3.6.3-3,  we  have  the  case  of  two  parallel  targets  when  LDR 
techniques  are  not  used.  In  this  chart  the  two  targets  are  flying  closely  together 
in  parallel  at  only  1 5  m  apart.  With  this  separation,  the  small  target  is  masked 
by  the  large  target  for  the  length  of  the  run.  As  a  result,  the  MHT  algorithm  only 
hypothesizes  one  red  track.  Therefore,  the  small  target  escapes  detection. 

In  Figure  3.6.3-4  we  have  the  same  scenario  here  as  in  Figure  3.6.3-3. 
Two  targets  are  flying  closely  together  in  parallel  at  only  1 5  m  apart.  With  this 
separation  the  small  target  is  masked  by  the  large  target.  However,  with  LDR 
techniques,  the  MHT  algorithm  is  able  to  observe  the  small  target  in  the 
masking  region.  With  two  closely  spaced  observations  at  each  update  time,  the 
MHT  algorithm  hypothesizes  two  or  more  blue  tracks.  The  extreme  closeness  of 
the  observations  mean  that  they  will  fall  within  each  target's  covariance  gates, 
thus  making  it  difficult  for  MHT  to  hypothesize  two  definite  red  tracks. 

Figure  3.6. 3-5  shows  the  case  of  10°  target  crossing  when  LDR 
techniques  are  not  used.  In  this  chart  we  have  two  targets  crossing  each  other 
at  a  1 0°  angle.  Initially  they  start  out  at  about  (530,-20)  and  (530,0)  and  move 
to  the  left  of  the  chart.  At  first  the  targets  are  far  enough  apart  so  that  we  get 
observations  from  each  target.  However,  they  are  also  close  enough  so  that 
their  observations  fall  within  each  other's  covariance  gates,  making  it  more 
difficult  to  obtain  two  red  tracks.  Therefore  the  tracks  are  blue  or  green  in  color. 
At  about  (380,10),  the  targets  cross  and  are  very  close  together.  This  is  the  start 
of  the  masking  region  where  we  don't  observe  the  small  target.  As  a  result, 
MHT  only  picks  up  the  large  target  as  a  red  track.  At  about  (240,  20)  the  targets 
again  become  separated  enough  for  us  to  have  2  observations.  Here  MHT 
creates  several  hypothesis  tracks  in  blue  until  about  X=160  when  there  are  two 
tracks  in  light  green  until  the  end  of  the  run.  This  chart  shows  that  without  LDR 
techniques,  the  small  target  disappears  in  the  masking  region. 

In  Rgure  3.6.3-6  we  have  the  same  scenario  as  in  Figure  3.6.3-5  but  with 
LDR  techniqes  used.  The  difference  here  is  that  MHT  has  2  observations  in  the 
masking  region  due  to  LDR.  So  in  the  masking  region,  instead  of  tracking  only 
a  single  red  track,  MHT  hypothesizes  several  blue  tracks.  With  more 
information,  this  can  suggest  that  there  are  two  targets  flying  very  closely 
together. 
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In  Figure  3.6. 3-7  we  have  two  targets  crossing  at  10°  but  starting  from 
opposite  ends  of  the  chart.  The  targets  are  easily  tracked  even  when  they  are 
very  close  together.  The  targets  never  enter  the  masking  region  because  their 
dopplers  are  so  different  (head-on  approach).  Thus  LDR  techniques  are  not 
needed. 

These  results  show  that  having  the  extra  observation  from  the  small 
target  in  the  masking  region  using  LDR  techniques,  the  tracking  system  can 
hypothesize  about  the  existence  of  another  target.  In  some  cases,  as  in  Figure 
3.6. 3-2,  the  existence  of  the  cruise  missile  is  immediately  realized.  In  other 
cases,  while  the  MHT  algorithm  is  not  able  to  keep  clear  red  tracks  of  the  two 
targets  when  LDR  techniques  are  used,  it  is  able  to  hypothesize  weak  tracks. 
This  is  already  a  significant  improvement  over  the  situation  of  not  receiving  this 
second  target  observation. 
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Figure  3.6. 3-1 .  Lauching  of  a  missile  when  LDR  techniques  are  not  used 
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Figure  3.6.3-3.  Parallel  targets  when  LDR  techniques  are  not  used 
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4.0  CONCLUSIONS  AND  FUTURE  WORK 


The  results  of  this  program  indicate  potential  for  greatly  improving  small 
target  detection  and  track  in  the  presence  of  much  larger  targets.  The 
phenomenology  of  LDR  must  be  studied  using  physical  systems  to  determine  if 
the  nature  and  extent  of  the  phenomenology  indicated  is,  in  fact,  exploitable.  A 
follow-on  study  which  physically  models  realistic  targets,  measures  radar  cross 
section  of  the  individual  targets  for  comparison  with  the  LDR  prediction 
methods,  and  examines  co-target  data  collected  over  a  wide  range  of 
frequencies,  waveforms,  and  relative  position  is  a  logical  first  step  in  confirming 
the  results  of  the  LDR  study.  Adjustments  to  the  assumptions  used  for  the 
tracking  simulation  could  then  be  made  and  the  LDR  improvement  over 
conventional  tracking  methods  could  be  determined. 

If  the  LDR  results  are  confirmed  by  the  above  approach,  a  flight  test 
program  using  full  sized  targets  and  operational  systems  should  be  conducted 
to  collect  and  analyze  data  in  the  full-scale  environment  under  a  variety  of 
conditions. 
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6.0  APPENDIX 


BCSJVlcKle-IJioiLtatm^ 

First  of  all,  the  LDR  interaction  mode!  may  be  optimistic  in  predicting  the 
strength  of  the  interaction  signal.  Less  energy  may  be  scattered  from  the 
missile  than  predieted.  For  a  stand-alone  cylinder,  the  scattering  will  be  found 
in  a  cone  around  the  cylinder  in  the  specular  direction.  However,  the  cruise 
missile  cylinder  has  many  components  attached  to  it,  and  so  this  result  may  not 
be  true.  Scattering  will  probably  be  limited  to  only  certain  angles  in  the 
specular  direction.  Also,  since  the  specular  point  approach  was  not  used, 
shadowing  is  not  taken  into  account,  and  it  is  possible  that  for  many  aspect 
angles,  scattered  energy  from  the  cylinder  may  be  blocked  by  the  wings  of  the 
missile  or  other  components.  Thus,  the  LDR  interaction  mode!  may  be 
predicting  too  much  interaction  between  the  two  targets. 

Secondly,  the  model  relies  heavily  on  the  monostatic/bistatic  equivalent 
theorem  for  scattering  from  the  airplane.  As  stated  in  Section  3. 1.2. 1.6,  the 
bistatic  theorem  breaks  down  for  large  bistatic  angles. 
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MISSION 

OF 

ROMELABORA  TORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this. 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportabiiity; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro^magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Bectromagnetic  Technology, 
Photonics  and  Reliablirty  Sciences. 


